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ABSTRACT 


This  report  reviews  the  basic  prooerties  and  terrestrial  surface 
responses  of  microwave  sensors,  both  active  (radar)  and  passive 
(microwave  radiometer).  Appropriate  surface  responses  (bistatic 
coefficient,  emissivity  etc.)  are  defined  and  used  to  provide 
general  formulas  for  the  system  response  (receiver  power,  antenna 
temperature,  doppler  spectrum)  in  terms  of  system  geometry  and  configur¬ 
ation.  The  dielectric  properties  of  terrestrial  surfaces  are  reviewed, 
and  data  on  the  complex  dielectric  constant  of  rocks,  soils  and 
vegetation  are  tabulated.  Theoretical  models  for  surface  response 
are  provided  for  a  variety  of  surface  classes,  including  the  surface 
of  uniform  layers,  the  surface  with  large  scale  unduiation,  the 
slightly  rough  surface,  the  surface  of  individual  scatterers  (Lommel- 
Seeliger  surface)  and  the  empirical  models  (Lambert  surfaces). 

Finally,  the  report  includes  a  large  number  of  typical  surface  responses 
at  microwave  frequencies.  Measured  bistatic  scattering  coefficients, 
radar  return,  and  surface  brightness  temperature  for  many  terrestrial 
surface  classes  are  presented  and  interpreted  in  terms  of  the  influence 
of  such  parameters  as  surface  roughness,  dielectric  constant,  polarization, 
frequency,  angle  of  incidence  etc.  on  the  response. 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1.  Geometry  of  the  scattering  problem.  25 

2.  Bistatic  scattering  geometry.  26 

3.  Range  ring.  27 

4.  Beam  and  pulse  limited  geometries.  28 

5.  Isodop  contours  and  doppler  spectrum.  29 

6.  Geometry  of  a  radiometer  system.  30 

7.  The  doppler  spectrum.  31 

8.  The  centroid  concept  for  small  targets.  32 

9.  Relative  complex  dielectric  constant  of  water. 

(Typical)  45 

10.  Dielectric  constant  of  soils  vs  moisture  content.  46 

11.  Relative  dielectric  constant  of  vegetation  at  8.5  GHz.  47 

12.  Loss  tangent  of  fir  wood  as  a  function  of  frequency.  48 

13.  Dielectric  constant  of  fir  wood  as  a  function  of 

frequency.  49 

14.  (a)  "Slightly  rough"  surface  h  <<  x.  Scattered  beam- 

width  determined  by  transmitter  beanwidth;  -  predomin¬ 
antly  coherent.  84 

(b)  Undulating  "rough"  surface  h  »  x,  slopes  «  1. 
Scattered  beamwidth  -  r.m.s.  surface  slope;  specular 

(not  coherent).  84 

(c)  The  Lambert  type  of  "rough"  surface  with  large  re¬ 
entrant  cavities.  84 

(d)  The  Lommel-Seeliger  type  of  "rough"  surface  with 

many  independent  scatterers.  84 

15.  Geometry  of  the  scattering  problem.  85 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Fi gure  Page 

16.  Geometry  of  layered  surface.  86 

17.  The  small  scale  component  has  Gaussian  surface  height 
correlation  coefficient  with  p(r)  =  exp(-r2/ £2 )  (i.  is 
small  scale  correlation  length,  h  is  the  mean  square 
height  of  the  small  scale  process.  Solid  curve  k$=l , 
kh=0.2;  dotted  curve  ka=l,  kh=0.1.  The  complex 
dielectric  constant  55(1- j. 55)  corresponds  roughly  to 

sea  water  at  X-banci.  87 

18.  Backscattering  from  cylindrical  vegetation  model.  88 

19.  Lomnel-Seeliger  cross  section  vs.  frequency.  89 

20.  Optical  depth  vs.  frequency.  90 

21.  Photo  of  loam  with  stubble  and  rough  sand.  119 

22.  Bi static  scattering  cross  sections,  in  the  plane 

of  incidence,  for  smooth  sand  and  loam  120 

23.  Bistatic  patterns  for  smooth  sand,  dry  grass,  and 

rough  sand.  121 

24.  Bistatic  scattering  behavior,  in  the  azimuth  cone, 

for  smooth  sand,  loam,  and  soybean  foliage.  122 

25.  Experimental  and  theoretical  reflection  coefficients; 
solid  curves  are  theoretical  values  for  smooth 

surfaces  -  Reference  1.  123 

26.  Reflection  coefficients  measured  by  Sherwood  and 

Ginzton  for  a  tidal  flat  with  some  organic  material; 

solid  curves  are  theoretical  for  e  =  10  -  Reference  9.124 

r 

27.  Effects  of  polarization  on  the  azimuthal  bistatic 

scattering  pattern  for  smooth  sand.  125 

28.  Effects  of  polarization  on  the  azimuthal  bistatic 

scattering  pattern  for  soybean  foil  age.  125 

29.  Effects  of  polarization  on  the  bi static  cross  section, 

in  the  plane  of  incidence,  for  smooth  sand.  126 

30.  Effects  of  polarization  on  the  bistatic  cross  section, 

in  the  plane  of  incidence,  for  dry  grass.  126 


v 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure  Page 

31.  Radar  backscattering  cross  sections  for  three  smooth 

surfaces  (limestone,  asphalt,  and  concrete)  at  X-  127 
band. 

32.  Radar  backscattering  cross  sections  for  three  rubble 
type  surfaces  (sand,  pea  gravel,  and  crushed  stone) 

at  X-band.  (C.P.  -  Circular  Polarization)  127 

33.  Radar  backscattering  cross  section  for  large  blocks 

of  pumice  at  Mono  Crater,  California.  128 


34.  Pumice  blocks  at  Mono  Crater.  The  photograph  shows 
the  large  angular  blocks  of  pumice  that  underlie 
the  rugged  topography  near  the  summit  of  one  of  the 
Mono  Craters.  Note  the  large  size  of  the  vesicles. 


The  scale  is  in  inches.  129 

35.  Radar  backscatter  for  green  soybeans  at  X-band.  130 

36.  Radar  backscatter  for  green  oats  at  X-band  - 

Reference  6.  130 

37.  Measured  and  computed  difference  between  vertical 
and  horizontal  backscatter  for  "smooth"  limestone 

at  X-band.  1 31 

38.  Measured  and  computed  difference  between  vertical 
and  horizontal  backscatter  for  "smooth"  asphalt 

at  X-band.  131 

39.  Angular  dependence  of  the  radar  backscatter  for 

slightly  rough  gravel;  average  diameter  of  rounded 
gravel  =  1  cm.  -  Reference  4.  132 

40.  Radar  return  vs.  normalized  roughness  for  crushed 

stone,  gravel,  and  sand  surfaces  at  X-,  Ku-,  and 
Ka-bands  -  Reference  4.  133 

41.  Seasonal  changes  of  grass  at  X-band  -  Reference  3.  134 

42.  Seasonal  changes  of  grass  at  Ku-band  -  Reference  3.  135 

43.  Seasonal  changes  of  grass  at  Ka-band  -  Reference  3.  136 


vi 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure 

Page 

44. 

Effects  cf  rain  on  a  smooth  asphalt  road  at  Ka-band  - 
Reference  3. 

137 

45. 

Effects  of  rain  on  two-inch  grass  at  Ka-band  - 
Reference  3. 

137 

46. 

Frequency  dependence  of  radar  backscattering  for  a 
smooth  asphalt  surface. 

138 

47. 

Frequency  dependence  of  radar  backscattering  for  a 
smooth  sand  surface  -  Reference  4. 

139 

48. 

Frequency  dependence  of  radar  backscattering  for  a 
slightly  rough  gravel  surface;  average  diameter  of 
rounded  gravel  =  1  cm  -  Reference  4. 

140 

49. 

Frequency  dependence  of  radar  backscattering  for  a 
rough  crushed  stone  surface;  average  diameter  of 
stone  =  3  to  5  cm  -  Reference  4. 

141 

50. 

Frequency  dependence  of  radar  backscattering  for  a 
rough  grey  punice  surface. 

142 

51. 

Frequency  dependence  of  radar  backscattering  for 
a  vegetated  surface  (oats). 

143 

52. 

Polarization  dependence  of  the  radar  backscattering 
for  a  "smooth"  limestone  surface  at  X-band  - 
Reference  5. 

144 

53. 

Polarization  dependence  of  the  radar  backscattering 
for  a  smooth  sand  surface  at  X-band  -  Reference  4. 

145 

54. 

Polarization  dependence  of  the  radar  backscattering 
for  a  slightly  rough  gravel  surface  at  X-band  - 
Reference  4. 

145 

55. 

Polarization  dependence  of  the  radar  backscattering 
for  a  rough  stone  surface  at  X-band  -  Reference  4. 

146 

56. 

Polarization  dependence  of  the  radar  backscattering 
for  a  vegetated  surface  (oats)  at  X-band  - 
Reference  6. 

147 

57. 

Brightness  temperature  of  smooth  asphalt  at  10  GHz. 

148 

vi  i 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure  Page 

58.  Brightness  temperature  of  smooth  asphalt  at  35  GHz.  1A9 

59.  Brightness  temperatures  for  different  size  fractions 

of  limestone  at  13.5  GHzs  -  Reference  11.  150 

60.  Brightness  temperature  for  dense,  broken,  grey 

pimice  at  10  GHz.  151 

61.  Brightness  temperature  for  light  weight,  dark 

pimice  at  10  GHz.  151 

62.  Brightness  temperature  for  large  blocks  of  obsidian 

at  10  GHz.  152 

63.  Photographs  of  volcanic  materials  at  Mono  Craters, 

California  -  Reference  7.  153 

64.  Brightness  temperature  of  alfalfa  at  10  GHz.  154 

65.  Brightness  temperature  of  green  wheat,  in  head, 

at  10  GHz.  154 

66.  Brightness  temperature  of  green  oats,  in  head, 

at  10  GHz.  155 

67.  Computed  and  measured  brightness  temperatures 
for  coal  (near  Cadiz,  Ohio),  limestone  (near 
Marblehead,  Ohio),  and  asphalt  surfaces  at  10  GHz, 
vertical  polarization  (Data  normalized  to  300°K 

at  the  Brewster  Angle).  156 

68.  Brightness  temperatures  of  playa  sediments,  with 

variable  moisture,  at  13.4  GHz  -  Reference  11.  157 

69.  Brightness  temperature  of  three  pumices  (Mono 

Crater,  California)  t  X-band.  158 

70.  Brightness  temperatures  for  irrigated  and  non- 

irri gated  soybeans  at  X-band.  159 

71.  Brightness  temperature  of  smooth  asphalt  at  10  GHz 

and  35  GHz.  160 

72.  Brightness  temperature  of  lapilli  (volcanic  ash)  at 

10  GHz  and  35  GHz  (Mono  Crater,  California).  161 

vlii 


LIST  OF  ILLUSTRATIONS  (Contd. ) 


Figure  Page 

73.  Close-up  of  Tapi  Hi  surface.  Note  the  concentration 


of  coarse  pumice  lapilli  in  botches  and  bands,  and 
the  sparse  ground  vegetation  and  scattered  pine 
needles.  Scale  is  in  inches. 

162 

74. 

Brightness  temperature  of  "rough"  dark  pumice  at  10 
GHz  and  35  GHz  (Mono  Crater,  California). 

163 

75. 

Bistatic  scattering  cross  section,  in  the  plane  of 
incidence,  for  smooth  sand  at  X-band. 

164 

76. 

Bi static  scattering  cross  section,  in  the  azimuth 
cone,  for  smooth  sand  at  X-band. 

166 

77. 

Bi static  scattering  cross  section,  in  the  plane  of 
incidence,  for  smooth  loam  at  X-band. 

168 

78. 

Bi static  scattering  cross  section,  in  the 
azimuth  cone,  for  smooth  loam  at  X-band. 

170 

79. 

Bi static  scattering  cross  section,  in  the 
plane  of  incidence,  for  rough  sand  at  X-band. 

172 

80. 

Bistatic  scattering  cross  section,  in  the 
plane  of  incidence,  for  stubble  at  X-band. 

(4>  =  0  or  180°  except  where  marked). 

174 

81. 

Bistatic  scattering  cross  section,  in  the  plane 
of  incidence,  for  dry  grass  at  X-band. 

178 

82. 

Bi static  scattering  cross  section,  in  the  plane 
of  incidence,  for  soybean  foliage  at  X-band. 

180 

83. 

Bistatic  scattering  cross  section,  in  the  azimuth 
cone,  for  soybean  foliage  at  X-band. 

181 

84. 

Bi static  scattering  cross  section  of  the  sea  (sea 
state  3)  at  C-band,  vertical  polarization  - 
Reference  13. 

182 

85. 

Bi static  scattering  cross  section  of  the  sea  (sea 
state  1)  at  C-band,  vertical  polarization  - 
Reference  13.  (Depression  Angle  =  Grazing  Angle  = 
90°  -  e,). 

183 

ix 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure  Page 

86.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  2)  at  C-band,  vertical  polarization  - 

Reference  13.  184 

87.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  3)  at  C-band,  vertically-polarized  trans¬ 
mitter  and  horizontally-polarized  receiver  - 

Reference  13.  185 

88.  Bi static  scattering  cross  section  of  the  sea  (sea 
state  1)  at  C-band,  vertically-polarized  trans¬ 
mitter  and  horizontally-polarized  receiver  - 

Reference  13.  186 

89.  Bi static  scattering  cross  section  of  the  sea  (sea 
s*ate  2)  at  C-band,  vertically-polarized  trans¬ 
mitter  and  horizontally-polarized  receiver  - 

Reference  13.  187 

90.  Radar  backscattering  cross  sections  for  several 

slightly  rough  surfaces  -  Reference  3.  188 

91.  Radar  backscattering  cross  section  for  normal 

desert  at  X-band.  189 

92.  Radar  backscattering  cross  section  for  Pisgah 

Lava  Flow  at  X-band  and  Ku-band.  189 

93.  Radar  backscattering  cross  section  for  Bristol 

(Dry  Salt)  Lake  at  X-band  -  Reference  16.  190 

94.  Radar  backscattering  cross  section  for  Lavic 

Dry  Lake  playa  at  X-band  and  ku-band.  190 


95.  Lavic  Lake  Playa.  The  surface  of  the  playa  is 

almost  perfectly  flat  except  for  small  scale  rough¬ 
ness  generated  by  contraction  (dessi cation)  cracks. 

These  cracks  break  the  playa  surface  into  polygons 
of  three  distinct  sizes:  A  -  measured  in  terms  of 
1 0 ’ s  of  feet;  B  -  measured  in  terms  of  a  few  feet; 
and  C  -  measured  in  terms  of  a  few  inches.  In  the 
intermediate  distance  is  the  radar  truck  being  used 
to  make  calibration  measurements  and  in  the  far  dis¬ 
tance  to  the  left  are  the  black  basalt  lava  flows  from 
Sunshine  Crater.  The  view  is  to  the  northwest.  191 


x 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure 

Page 

96. 

Radar  backscattering  cross  sections  for  the  lunar 
and  planetary  surfaces. 

192 

97. 

Frequency  dependence  of  the  radar  backscattering 
cross  section  for  the  Moon  -  Reference  21. 

192 

98. 

Polarization  dependence  of  the  radar  backscattering 
cross  section  for  the  Moon  at  23  cm  -  Reference  22. 

193 

99. 

Polarization  dependence  of  the  radar  backscattering 
cross  section  for  the  Moon  at  68  cm  -  Reference  22. 

193 

100. 

Radar  cross  section  of  the  lunar  surface  vs.  angle 
of  incidence  for  Surveyor  V  (13.3  GHz). 

194 

101. 

Radar  cross  section  of  the  lunar  surface  vs.  angle  of 
incidence  for  Surveyor  VI  (13.3  GHz). 

194 

102. 

Radar  cross  section  for  the  lunar  surface  vs. 
angle  of  incidence  for  Surveyor  VII  (13.3  GHz). 

194 

103. 

Radar  backscattering  cross  section  for  wooded 
terrain  at  X-band. 

195 

104. 

Radar  backscattering  cross  section  for  Arizona 
farmland  at  X-band  -  Reference  16. 

196 

105. 

Radar  backscattering  cross  section  for  24"-36" 
green  soybeans  at  S-,  X-,  and  Ka-bands.  Shaded 
area  is  envelope  of  measured  curves. 

197 

106. 

Radar  backscattering  cross  section  for  !Q"-36" 
oats  at  S-,  X- ,  and  Ka-bandr>. 

198 

107. 

Radar  backscattering  cross  section  for  wheat  at 

297  MHz,  5.87  GHz,  and  9.375  GHz  -  Reference  25. 

199 

108. 

Radar  backscattering  cross  section  for  cities 
at  220  MHz  and  428  MHz. 

200 

109. 

Radar  backscattering  cross  section  for  cities 
at  X-band. 

200 

no. 

Radar  backscattering  cross  section  for  the  sea 
at  428  MHz,  vertical  polarization  -  Reference  27. 

201 

xi 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure 

Page 

m. 

Radar  backscattering  cross  section  for  the  sea 
at  428  MHz,  horizontal  polarization  -  Reference  27. 

201 

112. 

Radar  backscattering  cross  section  for  the  sea 
at  1.228  GHz,  vertical  polarization  -  Reference  27. 

202 

113. 

Radar  backscattering  cross  section  for  the  sea 
at  1.228  GHz,  horizontal  polarization  -  Reference  27. 

202 

114. 

Radar  backscattering  cross  section  for  the  sea 
at  4.455  GHz,  vertical  polarization  -  Reference  27. 

203 

115. 

Radar  backscattering  cross  section  for  the  sea 
at  4.455  GHz,  horizontal  polarization  -  Reference  27. 

203 

116. 

Radar  backscattering  cross  section  for  the  sea 
at  8.91  GHz,  vertical  polarization  -  Reference  27. 

204 

117. 

Radar  backscattering  cross  section  for  the  sea 
at  8.91  GHz,  horizontal  polarization  - 
Reference  27. 

204 

118. 

Radar  backscattering  cross  section  for  the  sea 
with  oil  slicks  at  8.91  GHz,  vertical  polarization  - 
Reference  27. 

205 

119. 

Radar  backscattering  cross  section  for  the  sea 
with  oil  slicks  at  8.91  GHz,  horizontal  polarization  - 
Reference  27. 

205 

120. 

Radar  backscattering  cross  section  for  the  sea 
at  13.3  GHz. 

206 

121. 

Radar  backscat\.r:  iq  cross  section  1  ■>'  the  sea 
at  3.8  GHz,  (a vet ... ^e  of  vertical  horizontal 

polarization)  -  Reference  29. 

206 

122. 

Radar  backscattering  cross  section  for  the  sea 
at  9.4  GHz,  vertical  polarization  -  Reference  24. 

207 

123. 

Radar  backscattering  cross  section  for  the  sea 
at  24  GHz,  vertical  polarization  -  Reference  24 

207 

xii 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure  Page 

124.  Radar  backscattering  cross  section  for  the  sea 

at  35  GHz,  vertical  polarization  -  Reference  24.  208 

125.  Radar  backscattering  cross  section  for  the  sea 

at  9.6  GHz,  vertical  polarization  -  Reference  30.  208 

126.  Radar  backscattering  cross  section  for  the  sea 

at  24  GHz,  vertical  polarization  -  Reference  30.  209 

127.  Radar  backscattering  cross  section  for  the  sea 

at  35  GHz,  vertical  polarization  -  Reference  30.  209 

128.  Radar  backscattering  cross  section  for  the  sea 

at  48.7  GHz,  vertical  polarization  -  Reference  30.  210 

129.  Effects  of  polarization  on  the  radar  backscattering 

for  the  sea  at  24  GHz  -  Reference  30.  210 

130.  Near-grazing  radar  backscatter  for  trees  (X-band 

and  S-band)  and  cotton  seedlings  (X-band).  211 

131.  Near-grazing  radar  backscatter  for  cultivated 
terrain  at  10  GHz  -  Reference  32.  Parameter  is 

per  cent  of  returns  exceeding  given  ordinate.  212 

132.  Near-grazing  radar  backscatter  for  wooded  terrain  at 
10  GHz  -  Reference  32.  Parameter  is  per  cent 

of  returns  exceeding  given  ordinate.  212 

133.  Seasonal  variation  in  the  radar  backscatter 
distribution  for  cultivated  terrain,  at  10 

GHz,  for  a  1.25°  depression  angle  -  Reference  32.  213 

134.  Seasonal  variation  in  the  radar  backscatter 
distribution  for  cultivated  terrain,  at  10 
GHz,  for  a  2.5°  depression  angle  - 

Reference  32.  214 

135.  Seasonal  variation  in  the  raday'  backscatter 
distribution  for  cultivated  terrain,  at  10 

GHz,  for  a  5°  depression  angle  -  Reference  32.  215 


xiii 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figiu-e 

Page 

136. 

Seasonal  variation  in  the  radar  backscatter 
distribution  for  wooded  terrain,  at  10  GHz, 
for  a  0.7°  depression  angle  -  Reference  32. 

216 

137. 

Distribution  of  average  values  of  a0  for  different 
types  of  terrain,  at  X-band,  horizontal  polarization, 
for  a  1°  depression  angle  -  Reference  33. 

217 

138. 

Distribution  of  average  values  of  a0  for  different 
types  of  terrain,  at  X-band,  horizontal  polarization, 
for  a  5°  depression  angle  -  Reference  33. 

218 

139. 

Distribution  of  average  values  of  a0  for  different 
types  of  terrain,  at  X-band,  horizontal  polarization, 
for  a  30°  depression  angle  -  Reference  33. 

219 

140. 

Brightness  temperature  for  asphalt  at  10  GHz  and 

35  GHz. 

220 

141. 

Brightness  temperature  for  glacially  polished 
limestone  at  10  GHz. 

221 

142. 

Brightness  temperature  for  Pittsburg  #8  Coal  at 

10  GHz. 

22? 

143. 

Deviation  envelope  and  mean  value  curve  for 
the  radiometric  temperature  of  Palo  Alto 

Marine  mud  at  13.5  GHz  -  Reference  34. 

223 

144. 

Deviation  envelope  and  mear  value  curve  for 
the  radiometer  temperature  of  Palo  Alto 

Marine  mud  at  37  GHz  -  Reference  34. 

224 

145. 

Deviation  envelope  and  mean  value  curve  for 
the  radiometric  temperature  of  LaJoila  Beach 
sand  at  13.5  GHz  -  Reference  34. 

225 

146. 

Deviation  envelope  and  mean  va’ue  curve  for 
the  radiometric  temperature  of  LaJo'ila  Beach 
sand  at  37  GHz  -  Reference  34. 

226 

147. 

Brightness  temperature  for  lapilli  (volcanic 
ash)  at  10  GHz  and  35  GHz. 

227 

xiv 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure 

Page 

148. 

Brightness  temperature  for  disced,  bare  soil 
(2"-3"  ridges)  at  10  GHz. 

228 

149. 

Radiometric  temperature  for  silt  loam  soil 
(Davis,  California)  at  13.4  GHz  and  37  GHz  - 
Reference  35. 

229 

150. 

Radiometric  temperature  for  sandy  loam  soil  with 

10%  vegetation  cover  (Rocky  Mtn.  Arsenal)  at 

13.4  GHz  and  37  GHz  -  Reference  35. 

229 

151. 

Brightness  temperature  for  dense,  broken  grey 
pumice  at  10  GHz  and  35  GHz. 

230 

152. 

Brightness  temperature  for  light  weight,  dark 
pumice  (block)  at  10  GHz  and  35  GHz. 

230 

153. 

Brightness  temperature  for  large  blocks  of 
obsidian  at  10  GHz. 

231 

154. 

Radiometric  measurements  of  playa  sediments, 
lava  flow,  and  lava  flow  beneath  playa  sedi¬ 
ments  at  13.5  GHz  -  Reference  34. 

232 

155. 

Radiometric  measurements  of  playa  sediments, 
lava  flow,  and  lava  flow  beneath  playa  sedi¬ 
ments  at  37  GHz  -  Reference  34. 

233 

156. 

Brightness  temperature  of  green  oats,  in  head, 
at  10  GHz. 

234 

157. 

Brightness  temperature  of  ripe  oats  at  10  GHz. 

234 

158. 

Brightness  temperature  of  green  wheat  at  10  GHz. 

235 

159. 

Brightness  temperature  of  green  sorghum,  in 
tassel,  at  10  GHz. 

235 

160. 

Brightness  temperature  of  green  alfalfa  at  10  GHz. 

236 

161. 

Bn'ghtness  temperature  of  irrigated,  green  soy¬ 
beans  at  1C  GHz. 

236 

162 . 

Brightness  temperature  of  irrigated,  green  soy¬ 
beans  at  10  GHz. 

237 

xv 


LIST  OF  ILLUSTRATIONS  (Contd.) 


I 


Figure 

Page 

163. 

Brightness  temperature  of  irrigated  soybeans 
at  10  GHz. 

237 

164. 

Brightness  temperature  of  irrigated,  ripe  soybeans 
at  10  GHz. 

238 

165. 

Deviation  envelope  and  mean  value  curve  for  the 
radiometric  temperature  of  Palo  Alto  marsh 
foliage  at  13.4  GHz  -  Reference  34. 

239 

166. 

Deviation  envelope  and  mean  value  curve  for  the 
radiometric  temperature  of  Palo  Alto  marsh 
foliage  at  37  GHz  -  Reference  34. 

240 

167. 

Radiometric  temperature  for  dense  marsh  vegetation 
(Seal  Beach  Tidal  Marsh)  at  13.4  and  94  GHz  - 
Reference  35. 

241 

168. 

Radiometric  temperature  for  pasture  gras  (Thousand 
Oaks)  at  13.4  GHz  and  94  GHz  -  Reference  35. 

241 

169. 

Radiometric  temperature  for  flower  stocks  (Oxnard 
at  13.4,  37,  and  94  GHz  -  Reference  35. 

241 

170. 

Radiometric  temperature  for  celery  crops 
(Oxnard)  at  13.4,  37,  and  94  GHz  - 
Reference  35. 

242 

171a. 

Brightness  temperature  of  the  Sal  ton  Sea 
at  1.55  cm  -  Reference  36. 

243 

171b. 

Brightness  temperature  of  the  ocean  at  low  and 
medium  wind  speeds  for  1.41,  8.31  and  19.34  GHz. 
Reference  39. 

244 

17?. 

Deviation  envelope  and  mean  value  curve 
for  the  radiometric  temperature  of  sea 
water  (Palo  Alto)  at  13.5  GHz  - 
Reference  34. 

245 

173. 

Deviation  envelope  and  mean  value  curve 
for  the  radiometric  temperature  of  sea 
water  (Palo  Alto)  at  37  GHz  - 
Reference  34. 

246 

xv  i 


l 


LIST  OF  ILLUSTRATIONS  (Contd.) 


Figure 

Page 

174. 

Comparison  of  mean  radiometric  temperatures 
at  13.5  GHz  -  Reference  34. 

247 

175. 

Compari son 
at  37  GHz  - 

of  mean  radiometric  temperatures 
■  Reference  34. 

248 

176. 

Comparison 
at  35  GHz, 

of  measured  radiometric  temperatures 
horizontal  polarization  -  Reference  38. 

249 

xvli 


TABLE  OF  CONTENTS 


Chapter  Page 

I  INTRODUCTION  1 

II  THE  BASIC  SENSOR  RESPONSE  3 

A.  The  Radar  Response  3 

B.  Temporal  Fluctuations.  The  Doppler  Spectrum  9 

C.  The  Microwave  Radiometer  12 

D.  Large  Extended  Targets  -  Bias  Errors  16 

E.  Small  Extended  Targets.  The  Centroid  Concept.  19 

F.  System  Limitations  Due  to  Tluctuatiohs~~in 

Ar.tenna  Pointing  23 


REFERENCES  33 

III  THE  CONSTITUTIVE  PARAMETERS  35 

A.  Introduction  35 

B.  The  Dielectric  Constant  of  Water  and  Ice  38 

C.  Rocks  and  Rowders  39 

D  s0iTs  "  1 

E.  Vegetation  41 


REFERENCES  50 

IV  THEORETICAL  MODELS  FOR  SCATTERING  AND  EMISSION  52 

A.  Introduction  52 

B .  The  Flat  Surface  54 

C.  The  Slightly  Rough  Surface  65 

D.  The  Rouch  Continuous  Surface  -  Physical  Optics 

Models  71 

E.  The  Composite  Surface  74 

F.  The  Lambert  Model  “  78 

G.  Scattering  by  Individual  Elements 

TLomrnelT-Seeliger  Scattering)-  79 


REFERENCES  91 


xvi  i  i 


TABLE  OF  CONTENTS  (Contd.) 


Chapter  Page 

V  EXPERIMENTAL  DATA  93 

A.  Introduction  93 

D.  Bistatic  Scattering  Data  95 

C.  Backscattering  Data  98 

D.  Brightness  Temperature  Data  103 

E.  Special  Surface  Categories  107 

F.  Backscattering  Dati>  108 

G.  Brightness"  Temperature  Data  115 


REFERENCES  251 


THE  RESPONSE  OF  TERRESTRIAL  SURFACES  AT  MICROWAVE  FREQUENCIES 


I.  INTRODUCTION 

The  purpose  of  this  report  is  to  summarize  the  electromagnetic 
response  of  terrestrial  surfaces  at  microwave  frequencies.  Although 
a  large  number  of  data  on  both  the  scattering  properties  (radar  return) 
and  the  thermal  emission  (radiometric  or  brightness  temperature)  of 
such  surfaces  has  been  collected  during  the  last  decade,  there  has  not 
been  available  a  comprehensive  review  which  could  serve  the  newer 
applications  of  remote  sensing  at  microwave  frequencies. 

A  significant  fraction  of  this  data  is,  of  course,  available  in 
the  several  "data  banks"  that  are  now  in  operation,  and  these  banks 
have  been  valuable  in  providing  rather  complete  records  of  surface 
response,  and  in  compelling  a  more  uniform  format  for  data  reporting. 

They  have  not,  however,  addressed  themselves  to  the  problem  of  providing 
a  critical  understanding  of  the  data  they  contain.  While  this  report 
cannot  aspire  to  duplicate  the  completeness  of  the  data  banks,  it  is 
our  hope  to  provide  an  appropriate  theoretical  and  conceptual  frame¬ 
work  for  understanding  the  significance  of  the  microwave  data  collections. 

This  kind  of  understanding  is  nov;  particularly  important  because 
of  the  widespread  interest  in  exploiting  microwave  censors,  for  a  variety 
of  purposes,  by  investigators  in  other  (i.e.,  non-engineering)  disciplines, 
such  as  geology,  agronomy,  etc.  Many  of  these  newer  users  of  the  data 
are  unfamiliar  with  the  nature  of  the  relationships  between  the  surface 
parameters  (roughness,  dielectric  constant,  etc.)  and  the  microwave 
response.  Even  such  specialists  as  radar  or  radiometer  system  designers 
may  also  benefit  from  a  more  systematic  discussion  of  the  role  of  the 
various  factors  which  control  the  instrument  response. 
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In  the  succeeding  sections,  we  shall  discuss  first  the  basic 
properties  of  microwave  sensors,  particularly  the  manner  in  which  the 
geometry  and  other  system  parameters  affect  the  system  output.  Then 
the  various  kinds  of  surface  models  are  reviewed  and  the  relation  be¬ 
tween  the  predictions  of  the  model,  the  properties  of  the  surface,  and 
the  measured  response  are  compared.  Finally  a  representative  selection 
of  data,  for  both  active  and  passive  sensors,  is  used  to  discuss  in  a 
more  empirical  way  the  range  of  surface  response  to  be  expected  from 
terrestrial  surfaces.  The  principal  emphasis  in  the  report  is  on  the 
average  properties  of  the  response,  although  some  discussion  of  the 
statistical  properties  is  also  given. 


2 


II.  THE  BASIC  SENSOR  RESPONSE 


Before  discussing  the  nature  of  the  interaction  between  sensor  and 
surface,  it  is  desirable  to  review  a  number  of  system  concepts  in  order 
to  introduce  appropriate  definitions  of  the  parameters  used  to  describe 
the  surface.  Because  a  number  of  textbooks  on  radar  systems ,[i ,2,3,4] 
and  radar[5]  and  radio  astronomy[6]  are  now  available,  no  extensive 
discussion  of  system  operation  will  be  given  here.  We  have,  however, 
attempted  to  mention  some  of  the  limitations  of  the  conventional  treatmen 
and  to  provide  operationally  useful  definitions  from  which  (e.g.,  by 
computer  simulation  if  need  be)  the  system  performance  iiiay  be  estimated 
under  any  conditions.  A  large  number  of  convenient,  approximate  formulas 
for  various  specialized  geometries  are  given  in  References  1-6. 

A.  The  Radar  Response 

The  most  useful  descriptor  of  the  active  microwave  response  of  a 
terrestrial  surface  is  the  average  bistatic  scattering  cross-section  per 
unit  area,  a°s(e..  ,es , 4>s )  for  angle  of  incidence  (with  respect  to  the 
mean  surface  normal);  scattering  angles  es,  <j>s;  incident  polarization 
state  specified  by  first  subscript  (i)  and  polarization  component  of 
scattered  radiation  specified  by  second  subscript  (s)  (see  Fig.  1). 

This  parameter  is  defined  as  follows.  Consider  a  plane  wave  of  in- 

A 

tensity  I  watts/meter  incident  with  unit  propagation  vector  n\  in 
the  xz  plane  and  polarization  state  .  That  is,  the  incident  electric 
field  is  E  p.  exp[-jkn.  •  r  +  ju,t]  where  f  =  u>/2ti  is  the  frequency; 

X  =  2-rr/k  is  the  wave  length;  IQ  =  | Eq | ZQ  and  ZQ  is  the  impedance  of 
free  space.  Let  that  component  of  the  field  scattered  by  an  element 
of  area  dA  (dA  is  the  projection  of  the  actual  surface  onto  the  xy 
plane)  which  has  polarization  state  p$  (where  ps  •  p*  =  p^  •  p*  =  1) 
have  an  intensity  dl  at  distance  r  in  the  scattering  direction. 

Then  the  bi static  cross-section  per  unit  area  is  by  definition 


where  the  symbol  <  >  indicates  that  the  measurement  is  repeated  many 
ti i«ies  with  different  members  of  an  ensemble  of  similar  surfaces,  and 
the  average  scattered  intensity  taken.  It  is  a  consequence  of  the 
reciprocity  theorem  that 


( a )  o*  ( 0 .  <t>  • ,  0  d>  )  =  a  .  ( 0  ,6  ,  0  •  d)  • ) 
v  '  isv  vi’  svs'  siv  s’vs  vi' 

where  the  incident  plane  wave  can  now  have  any  azimuth  direction  <j>^ . 

Since  the  scattering  cross  section  is  defined  in  terms  of  statistical 

average  power  densities,  there  are  certain  limitations  on  its  applicability. 

For  the  definition  to  be  consistent,  it  is  necessary  that  the  receiver  be 

in  the  far  field  of  the  scattering  pattern  of  the  element  dA  in  order  that 

2  2 

dls  be  proportional  to  l/r$.  This  implies  that  rg  >  2i  /\  where  £  is  a 
coherence  length  for  the  surface  structure.  Secondly,  the  surface  element 
must  be  sufficiently  large  and  irregular  that  the  (ensemble)  average 
scattered  power  is  proportional  to  dA.  A  smooth  sea  surface  near  normal 
incidence,  for  example,  fails  to  satisfy  these  conditions;  the  scattering 
contains  a  coherent  part  (i.e.,  a  part  which  remains  identical  in  ampli¬ 
tude  and  phase  when  one  member  of  an  ensemble  is  substituted  for 
another),  and  great  difficulty  is  encountered  in  trying  to  interpret  the 
results  of  experiments  over  such  a  surface  in  terms  of  the  parameter  oQ. 
Finally,  it  must  be  remembered  that  in  many  system  applications,  it  is 
important  to  know  the  statistical  properties  of  the  return  (i.e.,  the 
first  probability  distribution  of  aQ,  the  correlation  between  returns 
measured  at  two  different  adjacent  frequencies,  etc.).  For  most  surfaces 
it  turns  out  that  many  simple  statistical  properties  can  be  inferred 
from  a  knowledge  of  aQ  alone.  However,  determination  of  the  statistics 
for  some  surfaces,  particularly  for  higher  order  statistical  properties, 
requires  a  detailed  scattering  model,  or  actual  measurements. 
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To  determine  the  average  power  collected  by  a  real  radar  system 
viewing  a  surface  target,  (see  Fig.  2)  consider'  a  radar  antenna  at 
T  in  the  xz  plane,  a  distance  n  from  the  origin  0,  which  radiates 
power  PQ( t )  with  a  one  way  power  pattern  f^il^.  .  For  convenience 
the  origin  of  coordinates  is  chosen  to  lie  on  the  axis  TO  of  the  antenna 
pattern.  The  angle  is  between  TO  and  the  line  from  the  antenna  to 
the  surface  element  dA.  The  angle  g.  is  an  azimuthal  pattern  coordinate, 
conveniently  taken  as  the  angle  -Y-O-dA  projected  onto  a  plane  perpendicular 
to  TO.  The  receiving  antenna  is  at  R,  specified  by  r^,e  , «j>  with  pattern 
f r(i|»  ,£  ).  (Because  of  these  choices  of  angle,  which  are  appropriate 
for  high  gain  antennas  aimed  at  the  origin,  the  magnitude  of  the  pattern 
fr(0,£;r)  does  not  necessarily  attain  its  maximum  value  along  the  direction 
RO.)  The  ensemble  average  power  P  (t)  collected  by  the  receiver  at  time 
t  is  then 


(3)  P, 


A  G. 
.(r)  =  emr  1 


(4ir) 


1 


p0( t_ ( ri+rs )/ c){fi  U-j  »Ci)fr(4«r.?r) 


r 


1 

s 


where  the  integral  is  taken  over  the  entire  xy  plane,  and 

Aemr  =  maximum  collecting  aperture  of  receiving  antenna 
Gt  =  maximum  gain  of  transmitting  antenna 
ei,9s’^s  =  ^3Ca^  incidence  and  scattering  angles  at  dA. 
c  =  velocity  of  propagation 

The  polarization  specifications  in  this  equation  must  be  handled 
with  care.  The  equation  represents  the  actual  received  power  due  to 
an  actual  transmitting  antenna.  The  patterns  f  (i|;  )  and  f.(ij.-^)  do  not^ 
in  general  represent  pattern  measurements  made  with  a  fixed  polarization 
test  antenna  such  as  a  dipole.  They  represent  the  total  power  density 
(whatever  its  specific  polarization  state  may  be)  transmitted  in  a 
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given  direction.  By  the  same  token  the  subscripts  on  a.  cannot 
represent  arbitrary  polarization  states  but  must  correspond  to  the 
actual  polarization  states  of  the  two  antenna  patterns.  The  more 
general  case  of  arbitrary  transmitting  and  receiving  antenna  polari¬ 
zation  requires  *n  extensive  formal  development  to  handle  correctly, 
since  the  complex  scattering  amplitudes  for  the  surface  must  be  in¬ 
troduced.  Simplification  of  the  resulting  equation  depends  on  making 
plausible  assumptions  about  the  correlation  between  the  various  terms 
of  the  scattering  matrix.  For  the  purposes  of  this  report,  it  is 
sufficient  to  point  out  that  the  situation  remains  simple  if  the  more 
conventional  choice  is  made  of  specifying  two  orthogonal  polarization 
states  (e.g.,  horizontal,  p*j,  p^  in  the  xy  plane  perpendicular  to  n^ , 
n  respectively;  and  vertical  p.  =  n.  x  p^1  and  =  n  x 

Thus  if  horizontal  polarization  is  transmitted,  then  the  term 
If i f r°is 1  in  Ecl‘  (3)  must  &e  replaced  by 


(4)  4<*r«r>«hh  +  fifr° 


hv 


where,  for  example,  (see  Fig.  2)  fY(ij>.  ,5^)  is  the  power  pattern  of  the 
transmitting  antenna  as  measured  with  a  "vertical"  dipole  oriented 
along  pV.  The  pattern  factors  are  normalized  so  that  the  maximum 
value  of  f^U.j)  +  is  equal  to  unity  and  it  is  assumed  that 

the  four  components  of  the  scattering  matrix  a-  •  are  uncorrelated. 

In  most  cases,  well  designed  high  gain  antennas  have  almost  constant 
polarization  over  the  main  beam  and  first  few  side  lobes,  so  that 
only  one  term  in  the  sum  of  Eq.  (4)  will  be  important. 

It  is  often  convenient,  for  high  gain  antennas,  to  have  estimates 
of  the  value  of  the  integral  of  Eq.  (3).  If  the  antenna  patterns  can 
be  approximated,  at  least  over  the  main  beam  region,  by  f ^  )  = 

exp[-A(l-cos^^ )]  and  f  ( ip  )  =  exp[-B(l-cosi|>r)3,  one  can  show  that 
for  a  cw  system, 
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<5> 


„  2  4  2  2  4  2  2  2,  2  2  2  2  2 

C  =  [A  r  cos  8.+B  r.cos  e +A  Br-r  (cos  e.+cos  e +sin  e.sin  e„cos  <j>) 
Ls  11  s  is  1  s  i  s  s 


or  for  back  scattering,  r.  =  rg. 


(6) 


P  .  Po  Aemr  Gt  /  ’"’o 
r~  (4nr.  \Acos6i 


Equations  (5)  and  (6)  represent  the  received  power  levels  for  the  cw 
case.  When  pulse  radars  are  used,  it  is  of  interest  to  know  the  time 
variation  (i.e.,  pulse  shape)  of  the  ensemble  average  received  pulse. 

This  can  most  readily  be  computed  by  converting  dA  in  Eq.  (3)  to  an 
infinitesimal  area  contributing  to  the  return  at  a  specified  time, 
that  is  (see  Fig.  3)  a  range  ring  of  radius  p  and  vi dth  dp.  One  can 
now  utilize  the  relationship  between  "receiver  time"  t1  (defined  so 
that  t1  is  zero  when  the  leading  edge  of  the  pulse  reflected  from 
the  subradar  point  reaches  the  receiver)  and  the  angle  \p  between  the  z 
axis  and  the  line  of  sight  to  the  range  ring.  If  the  radar  is  at  altitude 
h,  then  Eq.  (3)  becomes 


6.  A 

(7)  P(t.)=_L_emr 
(4nh ) 


2tt  t' 

Po(f  -  &)  f2(*,*)o0(*,*)  cos3^(^J  dcf> 

O  0 


where  the  dummy  variable  of  integration  is  T,  and  iji  and  r  are  functions 
of  T,  viz 


(8)  cos  *(T)  =  (1  +  cT/2h)-1 
r(T)  =  h/cos  ijj(T). 
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When  a  constant  amplitude  pulse  of  power  level  P  and  pulse  length  t 
seconds  is  used,  and  both  the  antenna  patterns  and  oQ  are  independent  of 
i|>»  Eq.  (7)  becomes 


(9)  Pr(t') 


Po  Gt  flem 


8irh‘ 


P2(^)  0Q(ip)  COS3\l> 


(t'-T)U(t'-r) 


where  U(x)  is  the  unit  step  function  and  ip  is  again  a  function  of  T. 

In  this  form,  the  time  response  of  the  pulse  may  be  used  to  illustrate 
the  difference  between  the  so-called  "beam-limited"  and  the  "pulse- 
length-limited"  geometries.  If  t*  <  t  (see  Fig.  4a)  then  the  lower 
limit  of  the  integral  of  Eq.  (8)  is  zero,  and  the  area  of  ground 
contributing  to  the  return  power  at  receive  time  t'  is  a  circle. 

The  build  up  of  energy  in  the  receiver  is  independent  of  pulse  length 
but  may  be  limited  by  the  antenna  pattern  if  the  antenna  beanwidth  is 
smaller  than  the  angle  subtended  by  the  illuminated  circle.  This 
situation  is  referred  to  as  a  "beamwidth-limi ted"  geometry. 

On  the  other  hand,  when  t'  >  t,  the  effective  illuminated  area 

and  thus  the  magnitude  of  the  received  power  is  an  expanding  annulus 

whose  width  is  controlled  by  t;  this  is  the  "pulse-length  limited" 

geometry.  It  is  clear  from  Eq.  (8)  that  at  receiver  time  t1  the 

area  actually  contributing  to  receiver  power  is  that  subtended  by  the 

cones  ^(t1)  and  ( t ' -  t)  (see  Fig.  4b).  Thus  the  "effective"  area  is 

2  2 

Aeff  ■  tt (h  sec  '{i(t'))  -  jr(h  sec  ^(t'-x))  and  this  is  approximately, 

for  t'>>  x,  equal  to  =  2nh  tan  i^(t')  [(cx/2)/sin  <j>(t')].  This 
will  be  recognized  as  the  annulus  whose  outer  radius  is  fixed  by  the 
leading  edge  of  the  pulse,  and  whose  width  is  the  projection  of  a 
pulse  of  length  ct/2  onto  the  xy  plane.  Thus  it  is  often  said  that 
the  "effective"  pulse  length  is  t/2. 
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B. 


Temporal  Fluctuations.  The  Doppler  Spectrum 

The  instantaneous  power  P^( t)  returned  from  a  natural  surface  is 
not  the  constant  value  indicated  by  the  ensemble  average  F  (t)  of 
Eqs.  (3)  or  (7),  but  varies  in  time,  or  from  pulse  to  pulse  (unless 
both  the  surface  and  the  transmitter  are  immobile,  that  is,  unless 
both  the  transmitter  position  and  the  surface  are  "frozen").  The 
statistical  properties  of  Pr(t)  depend  on  the  nature  of  the  trans¬ 
mitted  signal,  the  velocity  of  the  sensor  platform  relative  to  the 
mean  surface,  and  the  local  motion  of  the  surface  itself  (e.g.,  the 
movement  of  leaves  and  branches,  or  wave  motion  over  the  ocean). 

The  actual  computation  of  statistics  of  the  return  signal  cannot  be 
carried  out  without  a  detailed  model  of  the  surface  scattering 
mechanism.  However  certain  statistics,  such  as  the  power  spectral 
density  of  the  scattered  signal  due  to  transmitter  motion  (in  the 
case  where  the  surface  is  frozen)  can  be  calculated  from  geometrical 
considerations  alone. 

Consider,  for  example,  a  transmitter  moving  with  velocity  v  with 
respect  to  a  fixed  motionless  surface  (see  Fig.  5),  and  emitting  a 
signal  of  frequency  fQ(t).  Then  the  frequency  of  the  signal  returned 
to  the  receiver  at  time  t  from  the  area  dA  in  direction  n^  is 

(10)  f(t)  =  fQ(t  -  2r/c)  (1  +  2  v  •  n./c) 

where  r  is  the  range  to  the  actual  surface  point  at  dA.  This  equation 
is  correct  to  first  order  in  v/c.  The  second  order  terms  (of  order 
(v/c)  )  must  be  computed  according  to  the  special  theory  of  relativity; 
expressions  of  the  type  (1  +  v  •  n. )/(l  -  v  •  n^)  f  /c  are  based  on  a 
false  analogy  with  acoustical  doppler  shifts  and  give  incorrect  second 
order  terms. 


The  lines  on  the  surface  for  which  f ( t )  is  constant  are  called 
isodops  and  in  the  case  that  fQ( t )  (the  transmitter  frequency)  is 
constant  (excluding  f.m.  doppler  and  short  pulse  systems,  for  example) 
they  are  found  as  the  intersection  of  the  cone  v  •  n.  =  constant  with 
the  surface.  Note  that  the  doppler  shift  in  no  way  depends  on  the 
slope  of  the  surface  at  dA,  but  only  on  the  direction  of  v  and  the 
line  of  sight  n^ . 

For  a  plane  surface  the  isodops  are  conic  sections,  namely  the 
intersection  of  the  xy  plane  with  the  cone  obtained  by  rotating  n^ 
around  v,  i.e.,  the  cone  n.  •  v  =  constant.  These  isodop  curves  are 
ellipses  if  the  cone  does  not  intersect  the  horizon,  and  hyperbolas 
if  any  part  of  the  cone  does  rise  above  the  horizon. 


The  average  doppler  spectrum,  i.e,  the  average  power  per  unit 
positive  frequency  interval  may  be  found  by  integrating  over  the  area 
between  two  isodop  lines.  Thus  if  the  isodops  chosen  correspond  to 
frequencies  f-j  and  f^,  with  interval  &f  =  |f2  -  f -j  | ,  then  the  spectrum 
5(f)  is.  for  the  backscattering  cw  case 


(11)  S(f)  =  ff  -  fl»-GtP°(f°) 


(4t0  Af 

lim  Af+0  Isodop' strip 

between  frequencies 
f|  and  f2 


f2(VM  ao(eP 

r 


In  general,  this  integral  is  impossible  to  evaluate,  although 
numerical  integration  for  certain  special  cases  are  available. 

Even  for  high  gMn  antennas  the  peak  of  the  doppler  spectrum 
does  not  correspond  to  the  doppler  shift  from  the  ground  point  il¬ 
luminated  by  the  axis  of  the  main  beam.  This  effect  is  largest  over 
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surfaces  such  as  the  ocean,  for  which  o(e)  varies  rapidly  with  e, 
and  is  responsible  for  the  land-sea  bias  effect  in  doppler  navigation 
systems . 

The  power  spectral  density  of  the  return  signal,  as  given  by 
Eq.  (11),  combined  with  the  assumption  that  the  return  is  a  Gaussian 
random  noise  signal,  provides  sufficient  statistical  detail  for  most 
radar  design  purposes.  The  bandwidth  of  the  signal  is  approximately 
equal  to  the  difference  between  the  highest  and  lowest  doppler  shifts 
occurring  in  the  area  illuminated  by  the  beam.  However,  certain  special 
purpose  systems  (e.g.,  those  depending  on  frequency  agility  to  detect 
ship  targets  in  the  ocean)  may  require  a  more  detailed  knowledge  of  the 
fluctuation  of  the  return  signal. 

If  the  surface  is  also  in  motion,  the  doppler  spectrum  is  broadened 
in  a  manner  that  must  be  computed  from  the  kinematic  behavior  of  the 
surface  scattering  elements.  The  bandwidth  of  this  "natural  line  width" 
effect  is  approximately  equal  to  the  doppler  shift  associated  with  the 
rms  velocity  of  the  random  motion  of  the  surface.  Let  the  surface 
element  at  x,y  have  a  "natural"  doppler  broadening  spectrum  Sn(f-fQ) 

(see  Fig.  5b)  with 


|  Sn(f)df  =  1 
o 

when  the  surface  is  illuminated  with  a  plane  wave  of  frequency  f 
coming  from  the  transmitter  position  at  T  with  the  transmitter  motion¬ 
less  (v  =  0).  Then  the  total,  broadened  spectrum  St( f )  when  the 
transmitter  is  also  moving  with  velocity  v  becomes,  for  a  cw  signal, 


(13)  St(f)  = 


Aemr  Gt  W 


(4tt)‘ 


f2(V5i>c0le!)  Wf*Voftl  " 

r 


xy  plane 
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This  expression  gives  the  doppler  spectrum  when  both  transmitter 
and  parts  of  the  surface  are  in  motion.  In  the  most  general  case,  Sn(f) 
may  also  depend  on  position,  incidence  angle,  etc.,  along  with  aQ(e.j). 
When  the  surface  is  motionless,  Sn(f-f  )  is  a  delta  function  6(f-fQ) 
and  Eq.  (13)  reduces  to  Eq.  (12),  since  the  integrand  is  zero  except 
over  the  isodop  specified  by  v  •  n^ . 

For  the  bi static  case,  the  doppler  spectrum  at  the  receiver 
can  be  computed  from  the  above  formulas  by  replacing  the  doppler 
shift  fQ(2v  •  n./c)  by  f  (v  •  rT  -  v  •  ns)/c  where  ns  is  the  unit 
vector  from  dA  to  the  receiver.  The  isodops  then  become  the  loci 
v.  •  (n.j  -  ns)  =  constant.  When  the  transmitted  signal  is  not  a  cw 
signal  but  a  short  pulse,  the  return  pulse  will  have  a  spectrum  given 
by  the  convolution  of  the  pulse  spectrum  with  the  doppler  spectrum 
given  above. 

C.  The  Microwave  Radiometer 

The  microwave  radiometer  (see  Fig.  6)  estimates  the  brightness 

temperature  of  the  thermal  radiation  incident  upon  it  radiated  by 

the  ground  surface  and  the  surrounding  atmosphere.  If  the  radiometer 

1/ 

antenna  is  illuminated  by  radiation  of  brightness  temperature  Td(*j>.  ,£.) 

Dll 

(where,  as  in  Fig.  6,  ^  and  £•  are  polar  and  azimuth  angles  for  the 

'  '  k 

antenna  power  pattern  f^iJsS)  and  Tg  represents  the  brightness  tem¬ 
perature  of  the  radiation  in  state  k  incident  on  the  antenna)  then  the 
corresponding  antenna  temperature  is 

2 

l  /fk(*,e)  TgU,ii)da 

04)  tJA>=  }tL-i - 

l  Jfku>5)<ta 

k=l  K 

where  e  is  the  angle  between  the  antenna  axis  and  the  surface  normal, 
the  superscript  j  identifies  the  nominal  polarization  of  the  receiving 
antenna,  and  the  subscripts  k  indicate  any  two  orthogonal  polarization 
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states,  one  of  which  is  customarily  taken  as  identical  to  state  j  of 
the  antenna.  In  well  designed  antennas,  only  the  design  polarization 
is  significant  over  the  main  beam,  but  the  orthogonal  polarized  terms 
may  contribute  a  non-negligible  fraction  of  the  total  received  power 
(e.g.,  2%)  in  the  side  lobes  and  back  lobes. 


The  antenna  temperature  Ta(ec)  is  the  sensor  output,  i.e.,  the 
quantity  measured  by  the  microwave  radiometer  system.  However,  the 
quantity  desired  is  usually  (vp=0)  =  T^(6Q),  that  is,  the  brightness 
temperature  of  the  radiation  incident  along  the  main  beam  direction 
of  the  antenna.  Thus  Eq.  (14)  must  be  regarded  as  an  integral  equation 
relating  the  desired  brightness  temperature  and  the  measured  antenna 
temperature.  Except  in  radio-astronomical  measurements,  the  inversion 
process  is  not  usually  carried  out.  It  can  be  performed  in  a  number 
of  ways  including:  (i)  the  "bootstrap"  method  in  which  T  ( e ,4>)  is 

d 

measured  over  4tt  steradians.  A  trial  brightness  T^( © ,4>)  is  assumed 
equal  to  the  measured  antenna  temperature  (T‘  =  T  (©,<)>))  and  this  is 

D  d 

inserted  into  Eq.  (14)  to  obtain  a  desired  value  of  r(e,<j>),  (the 
value  the  antenna  temperature  would  have  if  the  brightness  were 
T _ ( © ,4>) ) .  A  new  trial  brightness  temperature  T)^(e)  is  estimated 


T'(e)  and  the  process  is  repeated  until  sue- 

C 


from  Tj^(e)  =  2Tfl(e) 
cessive  trials  converge;  (ii)  the  Fourier  transform  method  in  which 
the  integral  in  Eq.  (14)  is  recognized  as  the  convolution  of  the 
pattern  function  and  the  brightness  distribution.  By  converting  ij >,£  to 
"rectangular"  co-ordinate?  x'y'  (e.g.,  the  coordinates  of  the  point 
where  the  line  of  sight  specified  by  intersects  a  plane  perpendicular 
to  the  axis  RO  of  the  main  beams)  and  essuming  pattern  multiplication 
holds,  Eq.  (14)  may  be  put  in  the  form  of  a  two  dimensional  convolution. 
Ir.  that  case  the  Fourier  transform  of  T^x'y')  is  the  ratio  of  the 
transforms  of  T  (x'y')  and  f(x'y').  Only  rather  smooth,  well  behaved 
antenna  patterns  are  suited  to  this  technique;  (iii)  a  side  lobe  cor¬ 
rection  procedure  in  which  estimates  are  made  of  the  beam  efficiency, 
and  the  power  entering  the  antenna  via  the  side  lobes  and  back  lobes. 
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This  is  the  simples l  method  ai  d  gives  reasonably  accurate  estimates 
if  the  beam  efficiency  is  h i oh . 


The  brightness  temperature  T^(0Q)  incident  upon  the  antenna  is 
closely  related  to  the  bistatic  scattering  coefficients  a°.(  .,6  .a). 

jK  1  S  S 

Consider,  see  Fig.  6,  an  antenna  viewing  a  surface  at  an  angle  eQ  from 
nadir.  The  brightness  temperature  along  the  direction  OR  is  the  sum 
of  two  contributions;  the  first  is  the  emission  from  the  ground,  equal 
to  a  brightness  temperature  ej(0o)Tg  where  ej(0o)  is  the  emissivity 
(for  radiation  in  polarization  state  j)  of  the  ground,  and  T^  is  the 
physical  or  "thermometer"  temperature  of  the  ground.  This  is  assumed 
to  be  uniform  to  a  distance  of  several  skin  depths  below  the  surface. 

The  second  contribution  is  the  sky  noise  (which  is  due  primarily  to 
atmospheric  absorption  at  microwave  frequencies)  reflected  by  the  ground 
towards  the  antenna.  If  the  sky  noise  has  a  brightness  temperature 
distribution  Ts ( es »4>s )  then  the  toUl  brightness  temperature  of  the 
incident  radiation  coming  from  direction  eQ  in  the  j  polarization 
state  is[l 5,16] 


ground  emission 


attenuation 


reflected  sky  radiation 

_ H _ 


[a^  .(e  ,9  ,4  )  +  a0.,  (e  ,0  ,d>  )] 
L  jj  o’  s  ys'  jkv  o  s  Ts' J 


dfls  1 
4ttCOS0q  j 


r 

+  |  Tm(p)  e(p)  [a( r)/a(p ) ]  dp 

0  atmospheric  path 

emission 
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where 


:j(9o>  ■  1  -  II  +  0jk<9o’6s'*s)] 


dns 

4ttcosF 


=  emissivity  =  1  -  albedo 


a(p)  =  exp[-  a(p')  dp'] 
o 

-  attenuation  of  the  atmosphere  in  nepers  from  the 
surface  to  a  point  a  distance  p  from  the  surface 
along  the  line  of  sight,  (i.e.,  p  =  Z  sec0Q) 


(16)  a(p )  =  power  attenuation  coefficient  in  nepers/meter  of 

the  atmosphere  at  a  distance  p  from  the  surface 
along  the  line  of  sight, 
r 

a(r)/a(P)  =  exp[-  j  a(p')  dp'] 
p 

Tm(p)  =  physical  (thermometer)  temperature  of  the  attenuating 
medium. 


e(p)  =  emission  coefficient  (per  meter)  of  the  atmosphere 

This  equation  includes  a  factor  a(r)  to  account  for  the  path  loss 
between  ground  and  antenna,  and  a  term 

r 

Tm(p)  •"dp 

0 

to  account  for  the  emission  corresponding  to  this  warm  absorbing  path. 
The  actual  surface  brightness  consists  of  two  terms,  T^  giving  the 
ground  emission,  and 
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f  E  °jk  '  Ts  •  dos 


giving  the  sky  radiation  reflected  towards  the  receiver. 


If  losses  (i.e.,  attenuation)  due  to  scattering  by  the  intervening 
medium  between  surface  and  sensor  is  negligible  compared  to  losses  due 
to  absorption,  then  one  may  assume  e(p)  =  a(p).  Otherwise  e(p)  is  equal 
to  that  part  of  the  total  attenuation  due  to  absorption.  The  above 
equation  assumes  that  the  sky  radiation  ( es » )  incident  on  the  ground 
is  unpolarized,  and  that  the  atmospheric  attenuation  and  absorption 
coefficients  a  and  e  are  the  same  for  either  of  two  orthogonal 
polarization  states. 

It  is  often  tempting  to  write  the  terrain  dependent  part  of 
Eq.  (15)  in  the  form 


£j(9o>  Tg  +  Ts  with  V  ’  ■  ej(eo> 


%  ^ 

where  p  and  Tg  represent  some  appropriate  average  values  for  the 
surface  reflectivity  and  illumination.  This  is  possible  for  a  perfectly 
flat  surface,  but  it  is  clear  that  for  the  general  surface,  no 
independent  definitions  of  *YS  and  p  can  be  made  consistent  with  the 
choice  p  +  e  =  1 . 


D.  Large  Extended  Targets  -  Bias  Errors 

There  are  a  number  of  configurations  in  which  a  relatively  high 
gain  (narrow  pencil  beam)  antenna  is  used  to  estimate  the  properties 
of  an  extended  surface.  Because  the  system  response  is  a  convolution 
of  the  antenna  pattern  with  the  target  response  (see,  e.g.,  Eq.  (3)  for 
the  back-scattering  response,  Eq.  (11)  for  the  doppler  response, 

Eq.  (14)  for  the  brightness  temperature  response)  a  number  of  "bias" 
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errors  can  occur  when  the  target  properties  change  rapidly  with  look 
angle.  It  is  often  convenient  to  make  first  order  corrections  for 
these  eff/1cts,  in  order  to  avoid  the  need  for  inverting  the  convolution 
integral . 

We  summarize  here  these  first  order  corrections  when  the  antenna 
pattern  is  assumed  to  be  of  the  form 

f(*)  =  e-A(,-cos*> 

since  this  approximates  a  Gaussian  beam  when  A  is  large.  The  half 
power  beam  width  (between  3  dB  points)  of  such  a  pattern  is 


eB  =  2  (2*n  2/A)*5 

radians,  and  the  directivity  D  is  approximately  D  ^2A  •  (A  »  1). 
In  this  case  if  the  beam  axis  makes  an  angle  eQ  with  respect  to  the 
mean  surface  normal,  we  find  the  radar  response  (c.w.)  to  be 


P  G.  A 

P  /n  \  -  o  t  emr 

r  0  (W 


TT0°(eo) 

AcoseQ 


where 


Si  "  2  Sec  eo 


2  a0(®0)  2 
7  sin^e0  -  2  -  (coteo/4)  -  y  (7  sin^eQ  -  1) 


“o"(ecl 


o'  o 
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Thus,  the  correction  due  to  the  biasing  effect  of  the  surface  and  the 

geometry  of  the  experiment  is  of  order  1/D,  and  does  not  go  to  zero 

even  for  an  isotropic  surface  (o'  =  o'"  =  0);  the  fixed  error  is  due 

4 

predominantly  to  the  1/r  weighting  of  the  scattering  contribution 
from  dA. 

For  a  radiometer  antenna  viewing  a  surface  for  which  the  total 
brightness  temperature  incident  on  the  antenna  is  Tb(eQ) ,  the  antenna 
temperature  is  given  by 


(17)  Veo> *  W  +  AT/D 

aT  =  Tg( 6q)  -  Tg(eQ)  (9  sin2e  -1)  seceocoseceo. 

This  assumes  that  the  beam  efficiency  of  the  antenna  is  100%.  That  is, 
it  does  not  take  account  of  energy  entering  the  antenna  via  side  and 
back  lobes.  If  the  actual  antenna  has  a  beam  efficiency  of  100F%, 
(i.e.,  if  the  actual  antenna  pattern  car.  be  written  in  the  form 

fU)  =  exp(-A(l-cosv))  +  fs(i|>) 

where  f$ ( 4>)  is  the  side  lobe  contribution,  and  Jf  dft//fdft  =  (1-F), 
then  the  bias  error  correction  is  just  F  times  as  large  as  that  given 
in  Eq.  (17).  For  the  radiometer  measurement  a  "grey  body"  surface 
(T^( p)  =  T£(e)  =  0)  gives  no  bias  effect,  since  (see  Eq.  (14))  Tb(e)  is 
weighted  only  by  dn. 

The  most  significant  of  the  bias  errors  occurs  in  estimating 
the  doppler  spectrum  of  a  radar  signal.  Because  of  the  difficulty  in 
computing  the  spectrum  by  straightforward  integrations  over  the  isodop 
contours,  it  is  simpler,  in  this  case,  to  estimate  the  bias  effect  by 
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computing  the  moments  of  the  doppler  spectrum  S(f).  We  define  the 

a.  i_ 

ntn  moment  of  the  spectrum  with  respect  to  the  transmitter  frequency 
fo  by 


Anf 


s(f)  (f  -  fQ)n  df/ 


S(f)  df 


If  now  Af  is  the  doppler  shift  associated  with  the  line  of  sight 
along  the  axis  TO  of  the  beam  (i.e.,  Af  =  2f  v/c  sine„  cos(n)  where  v 

OQ  0 

is  taken,  as  in  Fig.  7,  in  the  xy  plane  with  v-ty  =  v  cos(n)),  then 
relations  can  be  found  between  Af  (the  "nominal"  doppler  along  the  main 
beam  axis  line  of  sight)  and  the  moments.  The  first  two  moments  give 
the  relations 


Afo  =  A  f  -  2(v/c)  fQ  5 


(18) 


(Af0)2  =  2 1  A1  f 


0  “7  4vf  f 

At  I  +  (2vf0/c)Z  5^  -  6 


Af 


where 


a '  (e  ) 

6  =  (1/D)  2  — 7-t— r-  ccs 0  -  12  sine  ^cosn 

°o^°o  0  0 


^  =  (1/D)  (cos2eQ  cos2n  +  sin2n) 


Thus,  by  measuring  the  first  moment  of  the  doppler  spectrum  a  f,  and 
using  the  first  of  Eq.  (18),  the  nominal  doppler  Af  may  be  estimated. 

E .  Small  Extended  Targets.  The  Centroid  Concept. 

In  certain  situations,  (as  for  example  in  determining  the  position 
of  an  island  in  the  ocean)  one  may  wish  to  estimate  the  sensor  response 
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as  the  beam  passes  over  a  small  extended  target.  (That  is,  the  "target" 
is  small  in  the  sense  that  it  subtends  an  angle  comparable  to  or  smaller 
than  a  beamwidth).  If  aQ(o)  and  T^e)  are  known  for  all  elements  dA 
of  the  target  area,  (and  this  implies  that  both  aQ  and  may  be  rapidly 
varying  functions  of  x,y,  the  position  coordinates  of  dA) ,  then  the 
sensor  response  may  be  computed  for  each  point  of  the  antenna  trajectory 
as  it  passes  over  the  target  area.  It  is  often  more  convenient  to 
estimate  sensor  response  in  terms  of  parameters  which  can  be  computed 
independently  for  the  target  and  for  the  experimental  geometry.  This 
is  possible  by  introducing  the  centroid  concept.  If  we  consider  only 
the  c-oi  radar  and  microwave  radiometer  systems,  the  contribution  of  a 
target  area  (see  Fig.  8)  to  the  total  received  power  dP^  or  to  the 
antenna  temperature,  dT  ,  can  be  written 

a 


dTa<%> 


(19)  a. 


dPr(o0) 


aj1  |  f(x,y)  T^U.y,  0 1 ( x ,y ) ]  dxdy 

/fU)  da  -  antenna  beamwidth 


We. 


o0  [x.y,0'(x,y)]  f2(x,y)  ^ 
A  r 


e'(x,y)  =  local  angle  of  incidence  at  x,y, 


where  the  integration  is  over  the  target  area  A,  and  the  antenna  pattern 
f(ip,C)  is  written  as  a  function  of  x  and  y.  Now  these  equations  are  of 
the  form 

dR(e0)  =  C  j  S(x,y)  g(xy)  dx  dy 
A 

where  dR(eQ)  =  sensor  response  at  look  angle  eQ 
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xc  =  /  X  s(x,y,e0)  dA//  s(x,y,eo)  dA 

(20) 

yc  =  /  y  s  dA//  s  dA 

where  all  integrals  are  over  the  area  of  the  target.  Then 
dPr  =  g(xcyc)  [  /s(x,y,eo)da]  +  second  order  terms 
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and  the  sensor  response  is  written  as  the  product  of  a  target  term 
independent  of  transmitter  pattern  and  geometry,  and  a  geometrical 
term  g(xcy(.)  dependent  only  on  the  experimental  geometry.  Thus  for  the 
radiometer, 


(21) 


■ 


cose. 


f<Vc> 


[ /Tb ( e o )  dA]  +  second  order 


where  r  ,  oc  are  the  range  and  angle  of  incidence  for  the  centroid 
point,  and  the  centroid  is  defined  by 


x 

y 


(x,y,e0)  dx  dy/ /  T(x,y)  dx  dy 


Similarly,  the  radar  response  is 


(22)  dPr(e0) 


Wen.  f  <V*c>  [  ,  ..  r  , 

"wo7  0 


dA  +  second  order 


where  the  radar  centroid  is  determined  by 


°o^x,y,eo^  dA//  o0(x*y)  dA  • 

The  magnitude  of  the  second  order  terms  can  only  be  found  by  comparing 
the  exact  integration  (Eq.  (19))  with  Eqs.  (21)  and  (22).  The  corresponding 
pointing  errors,  (i.e.,  the  errors  incurred  by  assuming  the  defined 
centroid  lies  on  the  axis  of  the  main  beam  when  the  measured  system 
response  dR  is  a  maximum)  are  quite  small,  and  do  not  impose  a  serious 
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limitation  on  the  centroid  concept  under  most  circumstances.  Never¬ 
theless  it  is  usually  desirable  to  check  system  performance  by 
numerical  integration  in  systems  requiring  high  pointing  accuracy, 
particularly  if  the  target  has  a  very  irregular  distribution  of 
surface  properties  o°(x y)  or  T^x.y). 

F.  System  Limitations  Due  to  Fluctuations 

in  Antenna  Pointing 

It  is  not  always  appreciated  that  in  many  operational  microwave 
radiometer  systems,  the  limiting  factor  in  system  performance  may  be 
determined  not  by  the  inherent  precision  of  the  instrument,  but  by 
"noise"  introduced  by  fluctuations  in  look  angle,  or  by  the  fluctuation 
due  to  surface  inhomogeneities  or  to  atmospheric  path  loss  variations. 
Suppose,  for  example,  that  a  small  target  contributes  brightness 
temperature  Tt(e)  to  the  antenna  temperature,  and  the  antenna  under¬ 
goes  r.m.s.  fluctuations  of  (de)2  in  look  angle  due  to  vehicle 
oscillation  or  the  antenna  steering  control  system.  Then  trie  total 
fluctuations  in  antenna  output  are  approximately 


p 

Tt(e)l 

2  -l 

p 

AT  ^ 
d 

UT0)2  ♦ 

L 

30 

(de)2 

where  AT  is  the  inherent  sensitivity  of  the  radiometer, 
o  J 

sTo s  K 

K  =  constant  of  order  unity 
B  =  RF  bandwidth  of  radiometer 
t  =  radiometer  time  constant 
Tg^  =  effective  receiver  noise  temperature. 
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It  may  often  happen  that  the  second  term  is  of  the  same  order  or  larger 
than  the  first  and  the  supposed  precision  of  measurement  A Tq  cannot  be 
achieved  in  practice.  Similar  effects  occur  when  a  moving  antenna 
views  a  fluctuating  background  brightness,  or  when  the  attenuation  of 
the  atmosphere  fluctuates.  For  example  if  the  surface  has  a  brightness 
temperature  and  the  atmosphere  has  physical  temperature  Tatm  and 
attenuation  a((ka«l)  between  surface  and  antenna,  then  the  antenna 
temperature  is  approximately 


Thus,  fluctuations  in  T  due  to  fluctuation  in  a  are  approximately 

a 

<‘V2  -  <Tb  -  w2  ■ 

Since  (Tb  -  T  )  may  be  of  the  order  of  100°K  over  reflective  surfaces, 
fluctuations  in  a  of  the  order  of  0.04  dB  may  introduce  fluctuations 
of  order  1°K  in  T  .  Very  small  amounts  of  condensed  water,  as  in 
clouds,  can  produce  attenuations  of  this  order  at  the  higher  microwave 
frequencies. [14] 

Finally,  systematic  errors  can  occur  in  target-seeking  split- 
lobe  or  double  beam  systems,  in  which  the  background  temperature  differs 
for  the  two  beams  (e.g.,  the  temperature  due  to  the  sea  background  for 
an  island  seeking  system).  This  shift  in  null  temperature  is  of 
order  [ aTg (0)/ 39]  A6  where  AS  is  the  angular  separation  of  the  two 
beams.  Again,  the  error  incurred  by  assuming  the  target  is  located  in 
the  null  direction  may  be  larger  than  the  nominal  pointing  accuracy  of 
the  system  over  a  homogeneous  background. 
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PROJECTION  OF  OD 


Fig.  2.  Bistatic  scattering  geometry. 
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Fig.  3.  Range  ring. 
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(b)  t'  >  T  PULSE  LIMITED 


Fig.  4.  Beam  and  pulse  limited  geometries. 
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Fig.  6.  Geometry  of  a  radiometer  system. 
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Fig.  8.  The  centroid  concept  for  small  targets. 
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III. 


THE  CONSTITUTIVE  PARAMETERS 


A.  Introduction 

The  scattering  from  natural  surfaces  is  controlled  by  both  the 
structure,  or  roughness  of  the  surface,  and  the  electrical  properties, 
that  is  the  complex  electric  permittivity  and  magnetic  permeability. 
The  magnetic  properties  play  a  rather  small  role  at  microwave  fre¬ 
quencies  since,  except  for  a  few  minerals,  the  permeability  may  be 
taken  as  unity. 


To  define  the  various  parameters  used  to  describe  the  electrical 
properties  of  matter,  it  is  convenient  to  start  with  Maxwell's  equation 
for  a  source  free  region 


where  H  is  the  magnetic  field,  D  the  electric  displacement  and  J£  the 

total  conduction  current.  In  order  to  solve  for  the  fields  within  a 

body,  it  is  necessary  to  provide  a  connection  between  D,  J  and  the 

electric  field  E.  This  is  done  by  introducing  the  constitutive 

-12 

parameters,  e_  (the  permittivity  of  vacuum,  e  ^8.85  x  10"  farad- 
_  1  0  0 

meter  )  (the  real  relative  dielectric  constant  of  the  medium)  and 

a  (the  real  conductivity  of  the  medium  in  mhos-meter"^ ) .  When  the  fields 

vary  as  eju)t,  then  Maxwell's  equation  becomes 


(1 )  VxH  =  [jojegGp  +  a]  E  =  ;i-xe  ec  E. 

The  complex  proportionality  factor  ec  may  be  written  in  a  number  of  ways. 
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(2) 


ec  =  (e„  +  e / j uie Q )  (the  complex  dielectric  constant) 
ec  =  er  -  je^m  (the  real  and  imaginary  parts  of  ec) 
e  =  e  (1  -  j  tan  6)  (tan  6  is  the  "loss  tangent") 

C 


where?  tan  6  =  a/ cue  e  =  e.  /e  etc.  On  the  other  hand,  in  the  theory 
o  r  im  r 

of  metals  it  is  customary  to  emphasize  the  complex  conductivity  ac 
by  the  definitions 


(3) 


VxH  =  [a  +  j(aieoer)]  E  =  ac  E  =  (a  +  j  aim)  E 
°im  =  U)eoer  irna9inafy  Part  conductivity) 


In  fact,  at  any  frequency  one  can  measure  only  the  complex  number  e  or 
ac.  How  one  assigns  significance  to  the  real  and  imaginary  parts 
depends  on  how  detailed  a  model  is  used  to  describe  the  interaction  of 
radiation  and  mutter  at  the  molecular  level.  For  example,  the  relax¬ 
ation  model  for  sea  water  provides  different  mechanisms  for  the  losses 
associated  with  the  polar  molecules  and  with  the  ion  transport.  Thus 
for  such  models  one  often  writes 

ec  =  er  '  je"  '  j(adc/u)Fo) 

where  now  e.  is  the  sum  of  two  terms  e.  =  e"  +  a  .  /we„.  H?-e  a  .  is 
lm  lm  dc  o  dc 

the  d.c.  conductivity  of  the  substance.  At  any  other  frequency  only 

the  total  loss  e-  can  be  measured,  and  the  difference  between  the 
ini 

measured  loss  and  the  term  a .  /me  is  ascribed  to  the  polar  relaxation 

dc  o  r 

mechanism.  Naturally  for  sea  water  the  terms  e"  and  a.  /me,,  have 

dc  o 

quite  different  dependences  on  temperature,  frequency  etc.  so  that 
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in  this  case  the  separation  is  quite  plausible.  Nevertheless,  the 

interaction  at  any  frequency  depends  only  on  the  single  complex 

number  e  (k). 
c 

In  the  time  domain,  a  corresponding  connection  obtains  between 
E  and  D,  and  it  is  not  hard  to  show,  by  considering  Eq.  (3)  as  the 
Fourier  transform  of  the  time  dependent  equations  that 


D(t)  =  e  E(t)  +  h(-r)  E(t-r)  dx 

00 

0 


(4) 


,  t  +ju)T 

h(T)  =  u  j  C-c(-)  '  E^e  d 


e  =  lim  e  (u) 
00  C 


that  is,  D  is  the  convolution  of  E  with  the  "impulse  response"  of  the 

dielectric  where  the  latter  is  the  Fourier  transform  of  (e  (u)  -  e  ). 

c  00 

(The  infinite  frequency  part  is  subtracted  off  to  avoid  the 
appearance  of  5  functions  in  Eq.  (4)).  Because  a  dielectric  may  be 
regarded  as  a  causal  linear  system,  (h(t)  =  0  for  t  <  0),  er(oj)  and 
e.jm(u))  constitute  a  pair  of  Hilbert  transforms,  a  connection  sometimes 
referred  tu  as  the  causality  condition,  the  Kramers -Kronig  relation, 
or  the  dispersion  relation,  viz 


1  ?  Eim(a),)  da)' 


1  E-SmV^ 

r  /  \  -i  _  •  im 

L.£vA -  £  ] - - 

L  r  ®J  TT  U)-U) 


(5) 


£j  (a))  ~  ~  ~ 
im'  it 


(e  (a)')  -  £  ) 


OJ-O) 


du' 


These  relations  are  often  helpful  in  determining  the  behaviour  of 

or  e.  over  short  frequency  intervals, 
r  im 
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There  ^s  a  second  group  of  parameters,  originating  in  optics, 
which  are  also  used  to  describe  the  electrical  properties  of  matter. 

If  one  considers  a  plane  wave  propagating  in  the  x  direction  in  a 
medium  with  parameters  eoec,  then  the  spatial  dependence  of  the 
fields  may  be  written  as  exp[  -  yx  +  jut]  where  y  is  a  complex  number 
referred  to  as  the  complex  propagation  constant.  The  real  and  imaginary 
parts  are 

y  =  a'+jg  =  juijy0ecec  =  {2*/\q)  =  3 (2tt/Xq)  n£ 

a'  =  attenuation  constant  (nepers/meter) 

3  =  wave  number  (radians/meter) 

XQ  =  vacuum  wavelength  =  (2ir) 

complex  index  of  refraction. 

The  wavelength  in  the  medium  is  X  =  2-rr/e ;  tne  skin  depth  6  or  penetration 
depth  in  the  medium  is  6  =  1/a'  meters.  Note  that  the  power  attenuation 
coefficient  used  in  Eq.  (16)  of  Chapter  I  is  a(p)  =  2a\ 

B.  The  Dielectric  Constant  of  Water  and  Ice 

Because  of  the  overwhelming  influence  of  water  in  modifying  the 
dielectric  constant  of  many  terrestrial  surfaces  at  microwave  frequencies, 
it  is  desirable  to  have  simple  estimates  for  ec.  These  may  be  obtained 
from  the  formula[l] 

e  j-e 

.  d _ .  p 

Gc  "  ri+jf/'f0)  "  J 

which  is  applicable  to  both  water  and  ice.  Here 
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a  =  ionic  conductivity.  This  is  of  order  4  mhos/meter  for  sea 
water,  and  is  proportional  to  the  salinity,  a  is  negligible 
for  distilled  water, 
e  =5.5  for  water  and  3.2  for  ice 

oo 

ed  =  87.7  -  0.4(T-273)  for  water  (Tin  °K) 

£(J  -  90  +  T/550  for  ice  °C) 

f  =  9.1  GHz  for  water  at  273  °K 

f  =  12.6  GHz  for  water  at  283  °K 
o 

f  =  17.2  GHz  for  water  at  293  °K 
o 

f  =21.6  GHz  for  water  at  303  °K 
o 

f  =  10  KHz  for  ice  at  273  °K 
f  =  3  KHz  for  ice  at  263  °K 
fQ  =  1  KHz  for  ice  at  253  °K. 

The  dissipative  properties  of  ice  are  not  well  represented  by 
this  formula  at  microwave  frequencies,  particularly  for  sea  ice,  the 
properties  of  which  can  vary  widely  depending  on  the  conditions  under 
which  it  was  formed.  For  this  reason,  measured  values[2l  of  loss  tangent 
should  be  used  for  ice  at  microwave  frequencies.  The  general  features 
of  the  elestrical  properties  of  water  (and  ice)  are  shown  in  Fig.  9„ 
in  which  is  plotted  against  e^,  with  frequency  as  the  parameter. 

C.  Rocks  and  Powders 


Rocks  and  minerals  in  the  solid  form  have  dielectric  constants 
in  the  range  of  er  =  2  to  10  at  microwave  frequencies,  and  show  a 
tendency  for  the  dielectric  constant  to  be  proportional  to  density  and 
almost  independent  of  frequency.  The  loss  tangert  for  such  materials, 
however,  can  range  over  several  orders  of  magnitude,  and  may  vary  widely 
with  frequency,  particularly  if  the  material  contains  water  inclusions. 
Penetration  depths  lie  in  the  range  of  0.5  to  10  meters  for  frequencies 
between  0.5  GHz  to  35  GHz, [5]  for  dry  materials.  Tables  I  and  II 
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summarize  a  number  of  measurements  reported  in  References  1,  4  and  5. 
Reference  5  also  includes  calculated  penetration  depths  for  the  samples 
listed  in  iable  II. 

When  rocks  are  ground  into  powder,  the  properties  of  the 
powder  car.  be  obtained,  approximately  from  the  powder  density  and  the 
bulk  dielectric  constants  of  the  solid  material.  A  convenient  general 
mixing  formula  has  the  form  of 

V£i  .  v  £b~£l 

£m+u  ‘  cb+u 

where 

e  is  the  complex  dielectric  constant  of  the  mixture 

ei  is  the  complex  dielectric  constant  of  one  component  (typically 

air,  e.j  =  1  for  powders) 

eb  is  the  bulk  complex  dielectric  constant  of  the  other  component. 

V  is  the  volume  fraction  of  the  total  volume  occupied  by  com¬ 

ponent  eK. 

u  is  a  form  number,  0  <  u  <  ®. 

The  form  number  has  a  particular  value  for  mixtures  of  a  particular  type, 
depending  principally  on  the  shape  of  the  particles.  A  value  u  =  2  has 
been  found  to  give  reasonable  results  for  powders  of  low  density[5] 

(-  1  gm/cm  )  whereas  u  =  2eb  appears  to  give  better  agreement  for  higher 
density  "mixtures"  such  as  sand.  It  may  be  of  interest  to  note  that 
the  dielectric  constant  of  powdered  rocks[5]  seems  to  depend  almost 
entirely  on  the  powder  density;  some  20  of  the  samples  from  Table  II 
exhibit  an  in  the  range  1.9  to  2.1  when  powdered  to  a  uniform  density 
of  I  gm/cm”3. 
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Soils 


The  dielectric  constant  of  soils  at  microwave  frequencies  depends 
predominantly  on  their  moisture  content.  Fig.  10  shows  measurements  of 
er  vs  moisture  content  for  two  common  soils. [3]  The  imaginary  part 
eim  1S  a1so  s^own  f°r  t*1e  0-3  SHz  case.  It  is  clear  from  these  figures, 
and  other  similar  data, [2, 6]  that  both  the  real  and  the  imaginary  parts 
of  the  dielectric  constant  for  typical  soils  are  much  higher  than  would 
be  expected  as  the  basis  of  laboratory  measurements  on  dry  samples,  and 
at  any  given  frequency  are  approximately  proportional  to  water  content, 
and  independent  of  soil  type.  One  consequence  of  this  fact  is  clearly 
illustrated  in  Fig.  68  (Chapter  V)  which  shows  the  systematic  changes  in 
brightness  temperature  of  a  soil  surface  as  moisture  content  is  varied. 

At  frequencies  below  100  me  the  electrical  properties  of  rocks 
and  soils[8]  are  exceedingly  complex,  and  no  attempt  should  be  made  to 
extrapolate  the  results  giver:  here  to  lower  frequencies. 

E.  Vegetation 

The  dielectric  constant  of  vegetation  is  also  controlled,  at 
microwave  frequencies,  by  the  water  content.  For  example,  Fig.  11 
shows  the  real  and  imaginary  parts  of  the  dielectric  constant  of 
typical  leaf  material  at  8.5  GHz.  At  frequencies  near  X-Band  the  results 
can  be  roughly  represented  by  a  formula  of  the  type 

‘vegetation  -<F/2>Re<‘„>  ' 

where  F  is  the  fraction  of  water  by  weight  in  the  vegetation  and  e  is 
the  cowplex  dielectric  constant  of  water.  At  lower  frequencies  (say 
1  GHz  and  below)  the  ionic  conductivity  of  the  dissolved  salts  in  the 
leaf  must  also  be  taken  into  account  in  estimating  the  dielectric 
constant. [9] 
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A  similar  situation  obtains  for  the  woody  parts  of  vegetation. 

For  example, [10]  oak  and  wych  elm  wood  at  70  -  75%  moisture  content  have 
er  =  45  at  0.1  GHz;  Figs.  12  and  13  show  the  complex  dielectric  constant 
and  loss  tangent  of  fir  wood  for  a  number  of  moisture  contents. 
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Freq . 


TABLE  I 

(Data  from  Refs.  2  and  4) 


n  GHz 

er 

tan  6 

Density  Material 

35 

2.3 

0.08 

Red  Granite 

35 

2.4 

0.06 

White  Granite  Crushed 

35 

1.7 

0.012 

White  Pumice  Crushed 

3C 

1.6 

0.02 

Black  Pumice  Crushed 

L.G 

35 

0.175 

3.01 

Chondritic  Metorite 

14 

4.6 

0.006 

1.81 

Halite 

14 

3.8 

0.012 

1.56 

Halite 

14 

3.3 

n  nnn 

u  .uuy 

1.42 

Halite 

10 

4 

0.1 

Limoni te  (Coarse) 

10 

4 

0.01 

Limoni te  (Fine) 

3.5 

3.3 

0.0055 

1.11 

Magnesite  Hard  Packed 

B.5 

2.45 

0.002 

1.22 

Quartz  Powder 

3.1 

10 

0.02 

Basalt  (Hawaii)  Oven  Dry 

3.1 

10 

0.08 

Basalt  (Hawaii)  0.36%  Water 

3.1 

7  to  9 

0.1 

Granite  (Quincy) 

3.1 

8.4 

0.006  to  .018 

2.65 

Limestone  (Lucerne  Valley) 

3.1 

5.5 

0.001 

Rhyolite 

10 

4.8 

0.005 

2.45 

Basalt  (Vesicular) 

10 

4.4 

0.013 

2.63 

Biotite  Granite 

xO 

5.1 

0.081 

2.35 

Obsidian 

10 

4.8 

0.009 

2.74 

Olivine  Basalt 

10 

5.4 

0.086 

2.68 

Serpentine 

10 

5.0 

0.027 

1.62 

Volcanic  Ash 

10 

4.7 

0.017 

2.27 

Altered  Tuff 

10 

5.5 

0.016 

2.03 

Tuff 

10 

4.7 

0.01 

2.3 

Horn  Blende 

10 

3.0 

0.012 

0.78 

Mono  Pumice 

14 

2.9 

0.011 

1.63 

Desert  Sand 

14 

8.2  to 

8.6 

0.004  to  .02 

2.65 

Limestone  (Lucerne  Valley) 

14 

4.7  to 

6 

0.01  to  0.1 

2.35 

Asphalt 

14 

4.5  to 

5.2 

0.02  to  .06 

2.1 

Concrete 

10 

2.7 

0.008 

1.9 

Pumice  (Mono) 

10 

1.7 

0.016 

0.45 

Pumice  (Mono) 
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TABLE  II 
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All  samples  vacuum  dried. 


0  20  40  60  60 


Fig.  9.  Relative  complex  dielectric  constant 
of  water,  (typical) 
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MOISTURE  CONTENT  (PERCENT) 
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Fig.  10.  Dielectric  content  of  soils  vs. 
moisture  content  (Ref,  3) 
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11.  Relative  dielectric  constant  of 
vegetation  at  8.5  GHz.  (Ref.  7) 
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LOSS  TANGENT  TAN&= 


Dielectric 
(Ref.  10) 
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IV. 


THEORETICAL  MODELS  FOR  SCATTERING  AND  EMISSION 


A.  Introduction 

In  this  section,  a  number  of  models  for  the  scattering  and  emissive 
properties  of  natural  surfaces  are  reviewed.  Such  models  provide  an 
insight  into  the  mechanisms  which  control  the  scattering  from  a  variety 
of  surface  types.  Thus,  the  models  can  be  used  to  identify  the  par¬ 
ameters  which  are  likely  to  affect  the  microwave  behavior  of  the  surface 
in  a  significant  way.  As  in  section  1,  emphasis  will  be  placed  on 
theoretical  descriptions  which  are  suitable  for  computer  modelling  of 
the  sensor  response.  Thus,  general  formulas,  or  chains  of  relatively 
simple  formulas,  will  be  given  rather  than  a  large  number  of  special 
forms  derived  by  algebraic  manipulation.  Wherever  possible,  theoretical 
predictions  will  be  compared  with  measured  responses,  in  order  to 
illustrate  the  extent  to  which  the  models  may  be  useful. 

An  excellent  survey  of  the  theoretical  foundations  of  rough  surface 
scattering  is  given  in  the  recent  review  by  Barrick,  et.  al.,[l]  which 
supplements  the  older  survey  by  Beckman  and  Spizzichino.[2]  Many  "special 
case"  formulas  may  be  found  in  these  two  references. 

Before  discussing  scattering  laws  in  detail,  it  may  be  desirable 
to  review  briefly  some  general  ideas  about  the  classification  of  rough 
surfaces,  and  the  nature  of  the  scattering  from  the  several  classes. 

In  the  first  place,  a  surface  may  be  perfectly  flat,  in  which 
case  an  incident  plane  wave  is  reflected  specularly  as  a  plane  wave., 
with  a  reduced,  polarization  dependent  amplitude  given  by  the  wel  known 
Fresnel  reflection  coefficient.  If  such  a  plane  surface  has  superim¬ 
posed  on  it  some  small  irregularities  with  their  rms  deviation  from  the 
plane  (rms  height  -  h)  much  less  than  a  wavelength,  then  the  principal 
scattering  is  still  a  coherent,  specularly  reflected  plane  wave,  but  the 
power  density  is  reduced  by  a  factor 
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| R | 2  exp[-2(kh  cos  e^)^3 

p 

where  |R|  is  the  appropriate  Fresnel  coefficient  and  the  factor 
exp[-2(kh  cos  e.)2] 

represents  the  effect  of  roughness  in  diminishing  the  forward  scattered 
energy.  The  remaining  power,  (of  order  of  magnitude  |R|  2(kh  cos  e^)  ) 
is  scattered  incoherently  into  all  directions,  i.e.,  may  be  considered 
a  diffusely  scattered  contribution.  (See  Fig.  14a).  There  are  two  points 
to  be  made.  The  first  is  a  distinction,  not  always  observed,  between 
"specular"  and  "coherent".  The  specular  power  leaves  the  mean  surface 
in  the  same  direction,  approximately,  as  the  specularly  reflected 
power  from  a  flat  surface.  It  may  or  may  not  be  phase  coherent  with 
the  incident  wave.  The  coherent  part  of  the  specular  power  is  phase 
coherent  with  the  incident  wave,  that  is,  if  the  surface  is  removed  and 
replaced  by  another  member  of  an  ensemble  of  such  surfaces,  the  co¬ 
herently  reflected  power  will  have  the  same  amplitude  and  phase  at  any 
point  as  it  did  before.  The  incoherent  power  will  not,  in  general. 

The  slightly  rough  surfaces  are  usually  treated  by  perturbation  methods. 

The  second  point  to  be  made  is  that  the  factor 

-2(kh  cos  ej2 
e  1 

illustrates,  in  more  quantitative  form,  Rayleigh's  criterion  for  dis¬ 
tinguishing  between  a  smooth  and  a  rjugh  surface.  We  normally  say  a 
surface  behaves  as  a  smooth  (or  slightly  rough  surface),  and  thus  retains 
an  appreciable  fraction  of  coherent  power  in  the  scattered  beam,  if 

2 

exp(-2(kh  cos  e^)  )  =  1,  or  kh  cos  <  1/4. 

This  transitional  behaviour  is  illustrated  in  Fig.  40. 
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The  number  1/4  is  quite  arbitrary;  various  authors  have  chosen 
different  numerical  criteria  depending  on  their  particular  application. 

It  will  be  observed  that  near  grazing,  cos  0  «  1 ,  a  surface  acts  more 
nearly  like  a  smooth  surface  than  that  at  normal  incidence.  If,  on  the 
other  hand  (kh  cos  e  >>12),  one  i?.ay  call  the  surface  rough.  Then  two 
limiting  situations  occur.  In  the  first  the  surface  may  be  gently 
undulating,  i.e.,  with  radius  of  curvature  much  less  than  x.  In  this 
case  (see  Fig.  14b)  each  part  of  the  surface  (say  over  local  areas  of 
one  or  two  Fresnel  zones)  is  locally  smooth,  and,  from  the  point  of  view 
of  geometrical  optics,  reflects  the  incident  radiation  as  if  it  were  a 
tilted  flat  surface.  The  total  scattered  field  is  then  the  sum  of  the 
‘reflections  from  these  locally  flat  areas,  and  it  seems  reasonable  (and 
is  correct)  to  estimate  that  the  scattered  power  forms  a  specular  beam 
with  a  beamwidth  approximately  equal  to  the  rms  slope  of  the  surface. 
Furthermore,  since  the  scattered  field  is  the  sum  of  many  local  contri¬ 
butions,  it  will  not  be  phase  coherent  with  the  incident  beam.  The  total 

p 

power  in  the  specular  beam  is,  again,  approximately  |R|  times  the  power 
in  the  incident  beam.  Surfaces  of  this  type  are  usually,  in  practice, 
treated  by  physical  optics  methods. 

Finally  a  surface  may  have  many  areas  of  large  slope,  vertical 
or  re-entrant  walls,  etc.,  as  in  a  layer  of  crushed  rock,  or  may  be 
made  up  of  many  independent  scatterers,  such  as  the  leaves  of  vegetation. 

In  this  case  there  is  seldom  observed  any  forward  lobe  in  the  scattering 
pattern.  The  scattering  is  more  or  less  isotropic,  and  is  referred  to  as 
diffuse  scattering.  This  class  of  surface  is  most  often  treated  in  terms 
of  empirical  scattering  laws,  such  as  Lambert's  law  and  its  generalizations, 
and  the  Lommel-Seeliger  scattering  type  laws.  (See  Figs.  14c  and  14d.) 

B.  The  Flat  Surface 

The  simplest  "surface"  is  the  perfectly  flat  surface,  i.e.,  a  hcif- 

space  with  complex  constitutive  parameters  e  and  y  .  The  scattering 

c  c 
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properties  of  such  a  surface  depend  solely  on  e  ,  ;>  and  the  illumination 
pattern  of  the  antenna.  Simple  as  it  may  appear,  even  the  primitive  case 
of  dipole  illumination  (often  referred  to  as  the  "Sommerfeld  Problem") 
presents  formidable  difficulties,  and  has  generated  a  vast  literature. 
Little  of  this  is  relevant  to  microwave  scattering,  but  there  are  some 
important  applications  to  back  scattering  from  tne  ocean  surface  at  HF 
and  lower  frequencies. [3] 

For  high  gain  antennas,  and  geometry  not  too  near  grazing,  one 

may  consider  the  flat  surface  to  produce  a  reflected  beam  identical  to 

the  beam  from  an  antenna  at  the  image  position  but  diminished  in 

amplitude  by  the  appropriate  Fresnel  coefficient.  The  reflected  field 

is  coherent  with  the  incident  field,  since  replacing  a  flat  surface  by 

another  identical  one  does  not  change  the  boundary.  (Near  grazing,  the 

total  field  is  the  sum  of  the  direct  and  reflected  fields;  this  produces 

a  lobe  structure  in  the  fields  generated,  for  example,  by  ground  based 

radars.  Such  effects  are  not  considered  here. [2])  Although  the  reflected 

2 

field  intensity  falls  off  as  l/rs  ,  it  is  not  proportional  to  the 
illuminated  area,  and  thus,  one  cannot  define  a0  for  the  flat  surface. 
Instead,  it  is  customary  to  describe  the  scattered  field  in  terms  of  an 
effective  reflection  coefficient. 

While  the  flat  surface  does  not  possess  any  scattering  properties 
of  interest,  its  radiometric  properties  constitute  one  of  the  standard 
geometries,  and  provide  an  insight  into  the  behaviour  of  many  terrestrial 
surfaces.  The  flat  surface  also  affords  an  opportunity  to  discuss  a 
number  of  details  concerning  polarization  states  and  the  geometry  of 
specular  reflection.  Consider  a  flat  surface  with  unit  normal  n, 
specified  by  the  polar  angles  e,  41,  illuminated  by  a  beam  along  n\  (see 
Fig.  15).  Then  the  direction  of  the  reflected  ray  is  along  the  unit 
vector  ng ,  the  so  called  "specular"  ■'irection.  The  relations  between 
these  three  vectors  are  given  by 
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(1) 


n$  =  n.  +  2  cos  a  n 
n  =  (ns~nJ/(2  cos  a) 

where  a  is  the  local  angle  of  incidence 

(2)  cos  a  =  n-n^  =  n*n$ 
cos  2«  =  n. -ns 

These  relations  allow  the  specular  scattering  direction  n$  to  be  found 
given  n,  ru  etc.  If  explicit  values  of  0$,  $s  etc.  are  required  it 
is  only  necessary  to  write  out  the  unit  vectors  in  polar  form, 

n .  =  -  t  cos  0.  +  T,  sin  e. 

1  i  1  » \  ! 

(3)  n  =  i  cos  e  +  i  sin  0  sin  $  +  i  sin  0  cos  . 

J  -  b  J  J  J  A  J  J 

ft  =  i  cos  0  +  i  sin  0  sin  <(■  +  i  sin  0  cos  <t> 

Z  y  X 

Thus,  one  may  easily  find 

(4)  cos  es  ~  (ft..  +  2  cos  a  ft)  •  iz 

sin  <t>  =  (ft.  +  2  cos  a  ft)  •  i  /  sin  0 
s  i  y  s 

etc,  with  similar  formulas  for  cos  0,  sin  $  if  the  term  (n^  +  2  cos  a  ft) 
is  replaced  by  (ns-n^)/2  cos  a. 

Many  special  forms  of  these  relations  may  be  found  in  References 
1  and  4.  Here  we  note  only  that  when  ft  is  specified  and  cos  0$  =  1 
(so  that  sin  eg  is  zero)  then  $s  is  indeterminate;  it  is  desirable  to 
set  (j>_  =  7T  in  this  case.  Similarly  when  fts  is  specified  and  cos  e  =  1 
then  <p  is  indeterminate;  in  this  case  it  is  desirable  to  choose  <j>  =  0. 
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There  r>  one  other  feature  of  plane  surface  reflection  of  interest 
in  connection  v/ith  the  randomly  tilted  surface.  If  the  normal  n  is 
allowed  to  lie  in  a  small  solid  angle  da,  then  the  reflected  ray  will 
lie  in  a  small  solid  angle  da,.  The  connection  between  the  two  solid 
angles,  i.e.,  the  Jacobian  of  the  transformation  between  e,  $  and 

V  V  is 

(5)  da  -  4  cos  a  da. 
s 

To  calculate  the  field  strength  of  the  reflected  ray,  one  must 
know  the  polarization  state  of  the  incident  ray  with  respect  to  the 
plane  of  incidence  (i.e.,  the  plane  of  in,  n  ,  and  n^).  However,  the 
polarization  state  is  often  specified  in  terms  of  "horizontal"  and 
"vertical"  unit  vectors  (see  Fig.  15).  These  are  defined  such  that 
h.  =  -  i  and  V  =  n  x  n . .  The  "vertical"  and  "horizontal"  unit  vectors 

*  ^  ^  -f 

for  the  reflected  ray  are  giver,  by  h  ■  n  x  i  /sin  e  (with  h  =  -  l 

sszs  sy 

if  es  »  0)  and  v$  *  h$  x  n$.  Now  any  arbitrary  incident  polarization 
may  be  written  in  the  form 


(6) 


Pi  =  a,  1., 


ta2  vi 


where  a^  and  a^  are  complex  numbers  such  that  P^’P-j*  =  I.  (For  example, 

for  linear  polarization  with  the  electric  field  making  an  angle  6  with 

the  "horizontal,"  a-j  =  cos  6,  b^  =  sin  6;  for  circular  polarization, 

a1  =  1/  a2  =  +  j/  The  state  of  polarization  orthogonal  to 

is  found  from  p^ -p^*  =0).  To  calculate  the  reflected  wave,  however, 

the  incident  polarization  p^  must  be  resolved  along  unit  vectors 

perpendicular  (s.  =  n .  x  n/sin  a  with  s.  =  -  i  when  a  *  0)  and  para lie 
^  •  I  I  y 

(ti  =  s^  x  n^ )  to  the  plane  of  incidence.  If  now  the  incident  electric 
field  has  the  form 


(7) 


-jk(n.  -r)+ju)t 
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then  the  reflected  wave  will  have  the  form 


-v  .  -jk(n  -r)+jait 

<8>  Erofl  ■  Eo  Psp  e 

where  the  polarization  vector  pgp  of  the  specular  reflected  wave  is 

(9)  p$p  =  {(iys.)  Rx(a)}  s$  +  {(fyt.)  Rj, (a) }  t$. 

Note  that  p$p  is  not  a  unit  vector,  but  contains  the  reflection  coef¬ 
ficients;  here  sg  and  t  =  s$  x  ng  are  the  "perpendicular"  and  "parallel" 
unit  vectors  for  the  plane  of  incidence  defined  by  the  specular  direction. 
The  Fresnel  coefficients  R^a)  and  R„(a)  for  perpendicular  and  parallel 

polarizations  are 

> 

2 

(y  COS  a  -  v  ye  -  Sin  n) 

(10)  ra(„)  =  -£ - -  -c-c-  -- 

(»c  cos  a  +  |wcsc  -  sinEo) 

(  rir 

f; _  COS  a  -  Jy,E„  -  Sin  a 

m«)  -  -  - -  = c 

.  2 

er  COS  a  +  'jycec  ”  Sin  a 

The  angles  n ■ ,  n.  betv zen  the  "vertical"  and  "parallel"  states  for  the 

I  ~4 

incident  and  reflected  fields  respectively  are  conveniently  found  from 
sin  n j  =  and  sin  ns  = 

If  the  amount  of  powe*  in  a  particular  polarization  state  is 
desired,  one  may  describe  the  state  by  p  =  b,  h  +  b0  v  where  b-,  and  b0 

5  I  S  l  S  \  c 

r8 


are  two  complex  numbers  (analogous  to  a-|  and  &2  for  the  incident 

polarization  state),  with  p  *p  *  =  1»  to  find  the  strength  of  the 

-*  s  s 

component  ps  of  the  specular  reflected  field  when  the  incident  field 
is  polarized  in  state  p^ ,  one  must  take  the  dot  product  PCp‘Ps**  Thus, 
one  may  regard 

<">  R1S<“>  ■  V’s* 

as  the  general  Fresnel  coefficient  for  reflection  from  state  p.  into 
state  ps<  For  example  the  strength  of  the  vertically  polarized  com¬ 
ponent,  when  the  incident  field  is  vertically  polarized,  is 

(12a)  Eyv  =  sin  sin  ng  Rx (ot)  +  cos  ni  cos  ns  R,,(a). 

Similarly,  one  can  show  the  left  circular  polarized  (left  cp)  component 
of  the  reflected  wave  when  right  cp  is  incident  is  given  by 

J  in.-  +  n.) 

02b)  Erl  =  Eq  e  1  s  (R^(o)  -  R,,{cx))/2. 

The  relations  above  permit  one  to  compute  the  amplitude  and 
phase  of  any  component  of  the  reflected  field  due  to  any  arbitrary 
incident  polarization,  if  the  surface  normal  is  specified.  Many  con¬ 
venient  special  cases  are  given  in  References  1  and  4,  but  for  computer 
evaluation  of,  s^y,  polarization  or  angle  of  incidence  effects,  a 
direct  sequential  computation  based  on  the  above  chains  of  formulas 
and  definitions,  will  often  prove  more  convenient  than  the  use  of  a 
single  complicated  formula. 

The  emissive  properties  of  a  flat  surface  art  as  simple  as  the 
reflective  properties,  as  long  as  the  incident  sky  illumination  is 
unpolarized.  In  this  case  let  the  sky  radiation  incident  on  the 
surface  from  direction  -n  (see  Fig.  15)  have  brightness  temperature 
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T  ,(e  ,<}>  ),  where  n  is  the  specular  direction  appropriate  to  the 
sxy  s  s  ^  s 

surface  normal  n  and  the  viewing  direction  n. .  Then  the  brightness 
temperature  of  the  radiation  traveling  along  -n.  is 


(13a)  T°  (0i)  =  Ts  (es  *s)  Rx(a) 


+  T  (1 
9 


Rx(oi) 


for  the  component  perpendicular  to  the  local  plane  of  incidence,  and 


2  2 
(13b)  T*  (0.)  =T$ky  (es  *s)  j R„(ot)  |  +Tg  (1  -  | R,,  (a)  |  ) 


for  the  component  parallel  to  the  plane  of  incidence,  where  the  direction 
0S  <f>s  ^efers  to  the  specular  direction  defined  by  n.  and  n.  In  general, 
one  wishes  to  know  the  brightness  temperature  TD(o.)  that  would  be 

D  1 

received  by  an  ideal  antenna  of  arbitrary  polarization  state  p^  (i.e., 
the  antenna  as  a  transmitter  would  produce  an  electric  field  of  tre  form 
of  Eq.  (7)).  This  is  found  from 


(14)  TB(e.)  = 


V*i 


-  I2  -rB, 


This  relatively  simple  representation  of  the  brightness  temperature 
"received"  by  an  antenna  with  polarization  state  p^  is  only  possible 
because  of  the  assumption  that  the  incident  sky  radiation  is  unpolarized 
(and  this  is  a  good  assumption  at  microwave  frequencies  where  the  sky 
radiation  is  due  solely  to  atmospheric  absorption).  In  the  case  of 
partially  polarized  sky  radiation  the  more  general  techniques  of  radio 
astronomy , [11 ]  utilizing  either  a  Stokes  parameter  or  a  coherency  matrix 
description  must  be  used. 

The  discussion  just  given  for  the  brightness  temperature  o^  a 
flat  surface  assumes  that  the  Kinetic  temperature  Tg  of  the  ground  is 


60 


uniform  down  to  a  distance  of  several  skin  depths.  If  the  Kinetic 
temperature  is  not  uniform,  but  is  given  by  T(z)  (z  <  0)  as  a  function 
of  depth  below  the  surface,  and  if  the  power  attenuation  coefficient  a(z) 
(i.e.,  a(z)  =  -  2  Im  [w  yQ  eQ  ec] 2)  is  also  a  function  of  depth,  then 
the  Kinetic  temperature  above  must  be  replaced  by 


o 

r 


aU)  ds/cos  0 


(15) 


eff 


Tg(z)  a(z)  e 


dz/cos  6, 


where  0r  is  the  angle  of  refraction,  i.e.,  the  angle  such  that  sin  ©r  = 
sin  O./ReU^) . 

The  next  order  of  complexity  in  the  flat  surface  is  the  surface 
composed  of  uniform  layers  (see  Fig.  16).  Here  the  cw  scattering 
properties  are  the  same  as  for  a  half  space,  i.e.,  the  beam  is  reflected 
coherently  as  if  emanating  from  the  image  point;  the  reflection  coef¬ 
ficient  must  be  computed  for  the  layered  structure  (see  below).  The 
brightness  temperature,  however,  particularly  when  the  layers  are  taken 
to  be  at  different,  (but  uniform)  temperature  must  be  computed  by 
summing  the  contributions  from  each  layer.  This  can  most  easily  be 
accomplished  by  the  use  of  a  reciprocity  principal.  This  asserts  that 
if  a  plane  wave  is  incident  on  the  surface  from  direction  0^  and  if  a 
fraction  f  of  the  total  incident  power  is  absorbed  in  a  certain  layer 
at  a  Kinetic  or  "thermometer"  temperature  T,  then  the  contribution  of 
that  layer  to  the  total  brightness  temperature  of  the  radiation  leaving 
the  surface  in  the  direction  e.  is  fT.  To  implement  this  computation, 
and  to  determine  the  required  reflection  coefficients,  an  iterative 
procedure  introduced  by  Richmcnd[5]  has  been  found  convenient.  Consider 
the  location  of  the  N  layers  to  be  specified  by  their  upper  surfaces, 
zQ  =  0  (the  bottom  "layer"  extending  to  z  =  -®),  2^,  z^  ••••  z^.  Each 
layer  has  thickness  d(i  and  uniform  Kinetic  temperature  T  ,  with  the 
boundary  between  the  surface  and  the  atmosphere  being  at 
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N 

ZN  =  |  dn‘ 

Consider  now  a  plane  wave  incident  from  the  direction  ,  i.e.,  along  * 
+n.. .  It  will  have  a  reflection  coefficient  R(e.. )  and  will  set  up  a 
combination  of  downward  and  upward  waves  in  each  layer.  Considering 
first  the  "horizontal"  or  "perpendicular"  or  !ITE"  polarization  the 
incident  wave  would  be  given  by 


(16)  E.nc  =  Tx  exp( jky  sin  e.  +  jkz  cos 

where  is  an  incident  electric  field  strength.  In  each  layer  the 
appropriate  fields  may  then  be  written  (e.g.,  for  the  layer  whose  upper 
surface  is  z(i) 


(17) 


(Ane 


-y  z  jky  sin  e. 
e  "  )  e  1 


-y  z  jky  sin  e. 
_  n  \  ^  i 

e  )  e 


where  ZQ  =  (,J0/e0)  2ien »  are  the  relative  electric  and  magnetic  per¬ 
mittivities  of  the  nth  layer,  and  the  complex  propagation  constant  is 
given  by 


(18)  yn  =  jk  (yp  en  -  sin2ei)Jj. 


By  enforcing  the  boundary  conditions,  one  finds[5]  that  the  coefficients 

St 

for  the  n+1  layer  can  be  found  from  those  of  the  :ayer  be^w  ;t, 
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(19) 


where 

P„  =  Cs)(l  +  P„)  sxp((vri  -  yn+1)zn) 
Q„  =  (*s)0  -  Pn)  exp((-Yn  -  Yn+1  )Z) 
(20)  Rn  =  (J$)(l  -  Pn)  exp((vn  +  Yn+i)z) 

Sn  -  (H)(l  +  pn)  exp({-  yn  +  Yn+1)z) 

pn  =  <Vl  Yn)/(|1n  Yn+1>' 


An  iterative  solution  for  the  relative  field  strength  is  easily  found 
by  setting  Aq  =  1 ,  Bq  -  0  (no  upward  wave  in  the  lower  halfspace  i  <  0; 
if  the  bottom  most  layer  is  a  perfect  conductor,  however,  one  should  put 
Bq  =  -  1).  Since  the  incident  electric  field  strength  is  A^+1  and  the 
reflected  field  strength  is  B^  ,  the  power  reflection  coefficient  is 


(21) 


R(ei) 


BN+1/AN+1 


The  next  problem  is  the  determination  of  f  ,  the  fraction  of  the  incident 
power  density  dissipated  or  absorbed  in  the  nth  layer.  This  is  seen 
to  be,  for  0  <  r.  <  N, 


(22) 


f-  =  (P 


Pn-1>/Pinc 


where 
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■ Re  2=z 

r 

Pinc  =  VrVl*/^ 


and  one  must  define  P_-|  to  be  P_-|  =  0.  If  now  the  sky  brightness 
is  given  by  (e^)  then  the  total  brightness  temperature  of  the  surface 
is  given  by 


,2  N 


(23)  Tb(0i)  -  Ts(ei)  R(ei)  +IfnTn 


Since  all  the  energy  not  reflected  is  absorbed,  we  must  also  have 


K  ■  i  -  «  • 


This  produces  the  standard  result  when  all  layers  are  at  the  same 
temperature.  For  the  case  of  "parallel",  "vertical"  or  "TM"  polarization, 
the  results  are  identical  except  that  now  the  fields  are  found  from  an 
assumed  incident  magnetic  field 

(24)  Hinc  =  \  Tx  exP(jk*  sin  0i  +  Jkz  cos  ei^ 
with  the  layer  fields  given  by 


x  ,  V  "V  3^  sin  ei 

H*  =  A  e  n  +  B  e  n  )  e  1 

n  n  n 


Fy  _  .  #  7  r,z  -Ynz  jky  sin  e. 

n  Yn  Zo/jken)  (An  e  *'  ‘  Bn  e  °  )  e 
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(26) 


pn  '  <En+1  Yn)/(En  Vl> 

Pinc  =  Vl  AN+1*  Zo 
Pn  '  Re  K  H„X*j  • 

The  expressions  for  R( ) ,  1^(6.)  etc.  are  identical  to  the  TE  case. 
Attempts  to  find  closed  form  expressions  for  the  brightness  lead  to 
formulas  of  great  complexity,  even  for  the  simple  case  of  two  layers. 
However,  by  computing  the  fields  layer  by  layer  upward  from  the  infinite 
half  space  representing  the  deepest  layer,  it  is  a  relatively  simple 
matter  to  find  the  actual  brightness  temperature  for  any  given  dis¬ 
tribution  of  layers.  One  interesting  case  which  has  received  some 
attention  in  the  li terature[l 2]  is  the  brightness  temperature  of  a 
layer  of  ice  forming  over  a  water  surface.  Here  oscillations  in 
brightness  of  over  100°K  were  observed  due  to  the  oscillations  in 
reflection  coefficient  of  the  layered  structure. 


C.  The  Slightly  Rough  Surface 

The  slightly  rough  surface  is  one  for  which  the  surface  height 
t(x,y)  above  the  mean  surface  satisfies  the  condition  c(x,y)/x  «  1; 
as/ax  «  1;  at/ay  «1.  That  is,  the  surface  height  is  every  where 
significantly  less  than  a  wavelength  and  the  surface  slopes  are  not  too 
large.  If  these  conditions  are  satisfied,  the  scattering  can  be  found 
by  a  perturbation  method  described  by  S.O.  Rice. [13]  The  first  order 
bistatic  coefficients  for  such  a  surface  have  been  worked  out  in  detail 
by  Bi.rrick;[l]  the  second  order  coefficients  have  been  obtained  by 
Valenzuela[6]  for  the  backs catte ring  case. 

The  general  feature  of  the  scattering  from  such  a  surface  is 
that  there  is  a  specular  (and  also  coherent)  reflection  process  identical 
to  that  occurring  from  a  perfectly  flat  surface  (see  Fig.  14a)  but 
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p 

diminished  in  amplitude  by  a  factor  of  order  exp(-  2(kh  cos  )  )  where  h 
is  the  rms  surface  roughness,  and  thus  provides  one  criterion  for 
judging  whether  the  surface  may  be  considered  "slightly  rough". 

The  bi static  coefficients  are  given  by [1 ] 

a°j  (ei,  es,  $s)  =  k4  h2  cos2ei  cos2es|aij|2  I 
=  TT2  k4  cos2ei.  cos2es|al.J.  |2  W 

k  -  2n/X 

2  7 

h  =  mean  square  surface  roughness  =  <c  > 

°r  I - 

I  ~  2tt  J  r  (r)  JQ(rk  J^2  +  5y2)  dr 
o 

I  =  (*2/h2)  W  (k  J  5X2  +  Sy2) 

E,  =  sin  e.  -  sin  e  cos  6  (2  sin  0.  for  backscatter) 

X  1  S  S  1 

C  =  sin  0  sin  4  (0  for  backscatter)  . 

y  s  s 

Mere  p(r)  is  the  correlation  function  for  the  surface  height 
p(r)  =  <  c(x,y)  c(x'y')  >  /h2 

(28) 

r  = 

and  W,  the  Hankel  transform  of  p,  is  the  roughness  spectral  density. 

The  matrix  elements  a.,  which  connect  the  polarization  states  of  the 

*  J 

incident  and  scattered  waves  are  given  in  the  general  case  in  Reference  1. 
for  the  simpler  case  yc  =  1  (i.e.,  a  non  magnetic  surface)  they  are 


(27) 

where 


(cr  -  1)  cos  4S 


*hh 


|  2  \  2 
(cos  6^  +  ec  -  sin  e.)  (cos  ©s  +  4  ec  ”  S1’n  6S) 


-sin  <t»s  (ec  -  1 )  \  ec 


sin2  0. 


vh 


!  2  2 
(cos  0.  +  4  ec  -  sin  0^ )  (ec  cos  0S  +  Jec  -  sin  0g) 


(29) 


sin  a 


(tc  -  11  J 


e„  -  Sin  0. 
c  1 


hv 


2  I  2 

(ec  cos  0i  +  jec  -  sin  0.)  (cos  0g  +  \  eq  -  sin  0g) 


(e_  -  l)[(er  sin  0.  sin  0C  -  cos  a 


s  >1  c 


e  -  sin  0. 


U  c 


-sin2  oc)] 


vv 


(ec  cos  0.  +  J  ec  -  sin2  0i)(e,  cos  e$  +  J  ec  -  sin2  0$) 


A  very  large  number  of  special  cases  of  these  formulas,  including 
those  for  perfect  conductors  (ec  ®),  for  circular  and  tilted  linear 
polarizations,  for  particular  types  of  correlation  functions  etc.  are 
given  in  Reference  1,  together  with  many  numerical  examples.  Here  we 
mention  only  two  of  these  special  cases.  First,  for  circular  polari¬ 
zation  states  L  (left  circular)  and  R  (right  circular)  the  matrix 
elements  become 


“LL  *  C(“hh  ‘  \v>  *  j  !“hv  +  “vh13 
RR 

(30) 

“LR  -  h  [(“hh  +  “vv’  1  j.Kv  '  “vh'3  ' 

PI 
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For  linear  polarization  states  p.  =  n.  sin  e.  -  v.  cos  3.  and 

->  lllll 

p  =  n  sin  3  -  v  cos  3  we  have  for  the  "direct"  polarized  return 

j  j  j  j  j 

6S  -  8j 

(3,a)  “siai  =  “w  cos2  6i  +  “hh  si"2  8i 

and  for  the  "cross "-polarized  return  3S  =  3^  +  tt/2 

(31t>)  \  6,  „  ,,  -  '  ’  “hh>  Sl"  6i  C0S  6i  ' 

The  scattering  in  any  direction  for  the  slightly  rough  surface 

is  proportional  to  the  power  spectral  density  of  the  surface  height 

appropriate  to  t.he  wave  numbers  ,  k?  .  These  correspond  to  the 

x  y 

Bragg  condi t  on  for  scattering  from  a  periodic  lattice  and  for  this 
reason,  particularly  in  the  literature  on  sea  return,  such  scattering 
is  referred  to  as  Bragg  scatter. 

There  have  been  two  useful  applications  of  the  slightly  rough 

surface.  The  first[7]  is  to  surfaces  such  as  asphalt  or  concrete 

runways  or  roadways,  w!..ch  clearly  satisfy  the  conditions  of  the 

theory  at  microwave  frequencies.  Here  if  the  surface  is  known  to  be 

2 

slightly  rough,  then  one  may  estimate  the  roughness  h  from  the  hack- 
scattering  coefficient;  the  dielectric  constant  ec  can  be  estimated 
from  the  difference  between  the  ’.'ackscattering  coefficients  for  vertical 
and  horizontal  polarization  (see  Figs.  37  and  38  of  Chapter  V).  The 
second  application  is  to  the  sea  surface, [3 ,8, 9]  where  the  capillary 
waves  at  the  lower  microwave  frequencies , [9]  and  the  entire  surface  at 
H.F.  and  lower  frequencies  satisfy  the  conditions  of  the  model. 

The  radiometric  properties  of  the  slightly  rough  surface  are 
essentially  the  same  as  those  of  a  flat  surface,  due  to  the  piedominant 
effect  of  the  coherent  reflection  process.  If  the  corrections  due  to 
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the  small  scale  roughness  is  required  one  must  compute  the  brightness 
in  two  parts.  The  first  is  the  "specular"  contribution  identical  to 
that  for  a  flat  surface  as  given  by  Eq.  (13),  but  with  the  Fresnel 
coefficients  |Rj  multiplied  by  an  additional  factor  F  to  account  for 
the  loss  of  energy  due  to  the  diffuse  scatter.  This  factor  is  of 
order  exp ( -2( kh  cos  e.)  )  but  must  be  computed  by  subtracting  the  total 
power  diffusely  scattered  into  the  upper  hemisphere,  i.e.,  the  factor 
is,  for  vertica1  polarization 


(32) 


Fy  =  1  -  (4ir  cos  ei)‘1 


[° vv ( 0 i  ’  V 

hemisphere 


°,h  (9i  es  +s> 


dn_ 


with  a  corresponding  expression  (interchange  v  and  h)  for  horizontal 
polarization.  To  this  must  be  added  the  "diffuse"  contribution, 


(33) 


Ald  =  0  "  Fv)  rg  +  (4tt  cose..}"1 


+  o  , )  T  ,  (e  *  ) 
vh'  skyv  s^V 


in  accordance  with  Eq.  (14)  of  Chapter  Ii.  Thus,  for  example,  a  flat 
surface  with  normal  n  in  the  z-di  recti  or,  will  have  a  brightness 


yB 


(6i>  Tsky  <%> 


*  T.  0  -  F„ 


)  *  ^  (.,) 


where 


The  spectrum  of  the  power  scattered  by  a  slightly  rough  surface 
has  been  derived  by  Barrick[14]  and  Valenzuela. [15]  For  the  back-scattering 
case,  the  spectrum  for  the  "natural  line  width",  due  to  motion  of  the 


69 


surface,  may  be  expressed[l  4]  in  terms  of  a  radar  cross-section  per  unit 
area  per  unit  bandwidth,  o  (f)  such  that  the  standard  cross-section  is 
given  by 


o 


0 


* 

O0  (f)  df. 


Here 


(34)  °0(f)  =  ^  cos^  9i 


WT(-2k  sin  9.  0,  2ir(f  -  f  ) 


where  f  is  the  transmitter  frequency  and  the  joint  spatial,  frequency 
spectrum  W-  is  given  by 


(35) 


(p»  q »  u)  ■  (h2/*2) 


oo 

. 


-00 


00 


-00 


°f  - jpAx  -  jqAy  -  juAt 

Rj  (ax, Ay, At)  e 

-00 

d( ax)  d(Ay)  d(At) 


where 


(36)  h2  R1  (ax.  Ay,  At)  =  <t(x,y,t)  c(x+AX,  y+Ay,  t+At)> 

and  c(x,y,v)  is  the  surface  height  as  a  function  of  x,y  and  t.  The 
exact  form  of  the  correlation  function  R-j  depends  on  the  dynamic 
processes  controlling  the  structure  of  the  surface.  If  one  assumes, 
for  example,  a  gaussian  height  and  time  correlation, 

R]  =  exp(-Ax2/£2  -  Ay2/£2  -  At2/!2) 

where  £  and  T  ure  the  correlation  length  and  correlation  time  respec¬ 
tively,  then  aQ(t)  becomes [14] 


70 


where  a  ( e - )  is  the  conventional  back  scattering  coefficient.  The 
0  1 

reason  for  referring  to  the  spectrum  due  to  surface  motion  as  the 
natural  line  width  is  now  apparent.  The  bandwidth  (in  this  case)  does 
not  depend  on  the  transmitter  frequency,  but  only  on  the  correlation 
Lime  T.  For  example,  for  the  capillary  scattering  of  microwaves  by 
the  ocean,  where  T  ^  1/5  sec,  the  bandwidth  is  of  order  2.5  Hz.  For 
the  HF  range  where  the  scatter  is  due  to  the  entire  structure  of  the 
ocean  surface  T  is  considerably  longer  (of  the  order  of  several  seconds) 
and  the  bandwidth  correspondingly  narrower,  viz  a  few  tenths  of  a  Hertz. 

D .  The  Rough  Continuous  Surface  -  Physical 

Optics  Models 

When  the  surface  is  of  the  gently  undulating  type,  of  continuous 
structure  with  radius  of  curvature  everywhere  several  times  larger  than 
an  electromagnetic  wavelength,  the  scattering  may  be  computed  by  physical 
optics  methods.  With  the  possible  excepticii  of  tK'  sand  dunes  in 
certain  barren  desert  areas,  r.o  known  natural  surface  satisfies  the 
conditions  of  this  model  at  microwave  frequencies,  aUhough  the  surface 
of  the  sea  might  be  considered  an  appropriate  one  at  optical  frequencies. 
Nevertheless  an  enormous  amount  of  theoretical  work  has  been  devoted  to 
this  model.  A  recent  review  is  g*ven  in  Reference  [15].  Although  the 
formulation  of  such  theories  has  taken  many  different  forms,  it  can  be 
shown[15]  that,  at  least  insofar  the  average  scattered  power  is 
concerned,  the  scattering  mechanism  is  equivalent  to  geometrical  re¬ 
flection  from  the  specular  points  of  the  surface  (see  Fig.  14b).  Thus  the 
cross-section  aQ  can  be  found  oy  determining  the  average  number  and 
radii  of  curvature  of  the  specular  points  from  the  statistical  properties 
of  the  surface,  or  even  more  directly  from  the  statistical  distribution 
of  surface  slopes,  since  at  each  specular  point  the  slope  must  be 


appropriate  for  specular  reflection  into  the  specified  scattering 
direction.  The  bistatic  cross-section  is  found  to  be 


(38)  (ei  ,es,<t>s)  =  tt  sec  0. 


Ris(a) 


p(e,<j>) 


it  sec  0. 


P'ix.ty) 


where 


p(e,<j>)  sine  de  d<j>  =  probability  that  surface  normal  n  lies 

in  dn  =  sine  de  d<j> 
c(x,y)  ■-  i^face  height 

p(c  ,c  )  dc  dc  =  joint  probability  that  surface  slopes  r  , 

a  Jr  A  j  A 

5y  lie  in  d;x  dCy 

n  ■  (tz  "  Vx  '  ^y  ^y^/( 1  +  +  4^ 

R^s(a)  =  Fresnel  coefficient  for  states  i,s. 

a  =  local  angle  of  incidence  defined  by  9^,  e$. 

In  evaluating  the  probabilities  above,  the  angles  e,<j>  for  the  surface 
normal,  or  5,  5  for  the  surface  slopes  are  to  be  chosen  such  that 
n  corresponds  to  a  direction  which  reflects  specularly  into  the  desired 
scattering  direction  65,^.  The  appropriate  relations  between  n  and  ng 
and  the  local  angle  of  incidence  a,  are  given  in  Eqs.  (l)-(4).  Similarly 
the  appropriate  Fresnel  coefficient  for  scattering  from  a  state  p^ 
into  a  polarization  state  p$  is  given  by  Eq.  (11). 

From  the  scattering  cross-section  of  the  physical  optics  surface, 
one  can  estimate  the  dielectric  constant  and  the  distribution  of 
surface  slopes.  It  is  customary  to  assume  particular  forms  for  p(e,4>) 
or  p( ;,?..);  a  large  number  of  special  cases  (gaussian,  exponential, 
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etc.)  are  given  in  References  1  and  15.  For  example,  for  backscattering, 

where  a  =  0,  c,  =  tan  e.,  c,  =  0,  the  direct  polarized  return  is 
x  i  y 


(39) 


O  _  O  _  0  _ 


11 


-  a 


hh 


-  a 


vv 


=  a 


RL 


=  [3  sec4  e.|/S2] 


Rx(0)  2  exp(-  /iT  tan  e^/S) 


for  the  exponential,  and 


(40) 


n 


[sec4  e./S2] 


Rx(0) 


2  exp[-tan2  e^S2] 


2 

for  the  gaussian  slope  distribution  models.  Here  S  is  the  mean  square 
surface  slope.  These  two  cases  are  plotted  in  Fig.  17. 


The  radiometric  properties  of  the  physical  optics  model  have 
been  worked  out  in  detail  by  Stogryn,[16]  who  has  given  a  number  of 
examples  of  application  of  the  model  to  the  sea  surface. [17]  Since  the 
bi static  coefficients  are  known,  for  any  given  slope  probability 
distribution,  it  is  only  necessary  to  substitute  the  expression  Eq.  (38) 
into  the  general  result  (Eq.  (15)  of  Chapter  II.).  Thus,  if  the 
effects  of  the  atmosphere  between  the  surface  and  ground  are  ignored, 
one  has 
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where  for  each  point  dftg  in  the  integration,  a,  e  and  <j>  must  be  found 

from  Eqs.  (l)-(4);  here  state  j  is  the  desired  or  receiver  polarization, 

and  k  is  the  orthogonal  polarization  state,  with  R..,  R..  from  Eq.  (11). 

J  J  JK 

The  evaluation  of  these  formulas  is  straight  forward  but  computationally 
tedious.  Examples  are  given  in  Reference  17.  The  results,  especially 
for  vertical  polarization,  are  very  similar  to  the  brightness  temperatures 
for  a  flat  surface.  The  model  appears  to  predict  the  brightness  tempera- 
turesof  the  ocean  at  microwave  frequencies  fairly  well  as  long  as  no 
significant  part  of  the  ocean  surface  is  foam  covered.  (See  Fig.  171b.) 


The  spectrum  of  the  signal  scattered  by  the  physical  optics 
surface,  aQ(f),  depends  on  the  motion  of  the  specular  points  as  the 
surface  changes  shape,  and  is  given  by [1 4]  for  backscattering, 


(42) 


C„(f) 


it  sec 


9i 
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sec  e- 
4tt 


(f  -  f0)  c 

tan  6i ’  °*  2  cos  6.  f 


where  here  p{?  ,  c  ,  5  )  is  the  joint  probability  density  function  for 
x  y  t 

the  surface  slopes  c  >  C  and  the  surface  vertical  velocity  c.  (ct  =  3c/3t); 

x  y  l  x 

that  is  p(c  ,  l  ,  c.)  dt  dt  dct  is  the  joint  probability  that  the  slopes 

A  II  L  A  y  w 

lie  in  dc  dc  ,  and  the  vertical  velocity  r.  of  the  surface  lies  in  dct. 

x  y  l  l. 


E.  The  Composite  Surface 

The  very  restricted  conditions  under  which  one  may  use  the 
slightly  rough  surface  and  the  physical  optics  models  have  led  to 
the  development[15,18]  of  the  composite  model,  in  which  the  total 
surface  height  c(x,y)  is  considered  to  be  the  sum  of  two  independent 
random  processes  u(x,y)  and  s(x,y),  i.e., 

dx,y)  =  u(xy)  +  s(x,y). 


74 


Here  u(x,y),  called  the  large  scale  component,  satisfies  the  conditions 
for  the  physical  optics  model  (radius  of  curvature  »  x).  Superimposed 
on  this  large  scale  structure  is  a  small  scale  structure  s(x,y)  which, 
relative  tc  the  local  surface  specified  by  u(x,y),  satisfies  the 
conditions  of  the  slightly  rough  surface  model,  i.e.,  slopes  ?s/3x, 
as/9y  and  normalized  height  s/x  all  «  1.  In  the  case  of  the  sea  at 
microwave  frequencies,  the  "sea"  and  "swell"  components  of  the  wave 
spectrum  would  constitute  the  large  scale  prccess,  and  the  capillary 
waves  or  ripples  would  constitute  the  small  scale  process.  A 
satisfactory  solution  for  the  scattering  from  such  a  surface  is 
quite  difficult,  but  to  a  first  approximation,  it  is  usual ly [1 3  assumed 
that  the  total  cross  section  for  the  composite  model  is  the  sum  of  two 
parts , 

(43)  o°(ei,  es,  $s)  =  cr°p(e.j ,  es,  <i>s)  +  <ajj(e.,  es,  *s)> 

where 

a°  =  scattering  cross  section  for  u(x,y)  process  alone, 

sp 

computed  trom  Eq.  (38) 

«jj>  =  scattering  cross-section  for  s(x,y)  process  computed 
from  Eq.  (27)  in  a  co-ordinate  system  referred  to 
local  normal  of  u(xy)  process,  and  then  averaged  over 
the  distribution  of  normals  specified  by  the  large 
scale  process. 

The  process  of  averaging  over  the  slope  distribution  must  be 
carried  out  numerically.  For  this  purpose,  one  chooses  an  incident 
angle  and  a  normal  n.  The  incident  wave  must  now  be  resolved  into 
parallel  and  perpendicular  components  in  the  local  surface  co-ordinate 
system  defined  by  the  local  normal,  and  the  vector  perpendicular  to 
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the  local  plane  of  incidence,  (i.e.,  the  unit  vectors  n,  ,  t.  in 
Fig.  15).  Now  the  local  angle  of  incidence  a  and  the  local  polarization 
state  (determined  by  the  angle  in  Fig.  15  for  linear  incident  polari¬ 
zation  states)  may  be  found  from  Eqs.  (1-4)  and  (6).  The  local 
scattering  angles  e^,  4/  are  then  found  from  cose^  =  ns  ■  n,  and  tan<}/  = 
-ns-s^/(ns  •  tj ).  These  angles  ej,  4/  can  now  be  substituted  into 
Eqs.  (27)  and  (29)  to  compute  the  contribution  of  the  slightly  rough 
component  to  scattering  from  n.  intc  n  direction,  when  the  slightly 
rough  surface  is  tilted  at  an  angle  specified  by  the  given  n.  The 
average  over  the  slope  distribution  is  then  carried  out  according  to 

(44)  <oij(0i’  V  4>s)>  =  /  o0(6.  0^  P(n)  dfi 

where 

p(n)dn  =  probability  that  the  normal  n  of  the  u(xy) 
process  lies  in  solid  angle  dn. 

p 

If  the  mean  square  slope  S  of  the  u(xy)  process  is  not  too 
large,  the  diffuse  contribution  for  the  direct  polarized  components 
<0vv>’  are  near]y  ^ie  same  as  the  unaveraged  values,  except 

for  horizontal  polarization  near  grazing  incidence,  where  the 
unaveraged  cross-section  decreases  much  more  rapidly  than  the  averaged 

The  principal  difference  for  the  composite  model  is  the  existence 
of  a  cross -polarized  term  <0^  >  in  the  backscattering  direction,  as  a 
consequence  of  the  fact  that  the  incident  polarization  is  no  longer 
perpendicular  to  the  local  plane  of  incidence.  Thus,  for  backscattering, 
for  example 
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(45) 


<0hv(81)> 


4tt  k 


avv(a)  "  ahh(a) 


sin2r,. 


2 

COS  T)j> 


W(-2k  sine. ,  0). 


Although  the  average  is  difficult  to  evaluate,  one  can  illustrate  the 

order  of  magnitude  of  the  cross -polarized  term  by  assuming  the  large 

scale  process  has  a  gaussian  slope  distribution  with  mean  square 
2 

slope  S  .  In  this  case 


(46)  <ahv^®i^''  —  2-rr  S2  k4  cos4  0. 


cosec  ei 


a  (e.)  -  a,  .( 0  • ) 
vv'  v  hh'  v 


W(-  2k  sine^ 


0 


The  only  significant  application  of  the  composite  model  has 

been  to  scattering  of  microwaves  by  the  ocean. [8, 9]  In  this  case 

certain  properties  of  the  ocean  surface  height  spectrum  permit  an 

explanation  of  a  number  of  features  of  the  back-scattering  behavior, 

in  particular,  (i)  the  weak  dependence  on  wavelength  (because  the  wave 

-4 

number  dependence  of  the  surface  height  spectrum  W(k)  ^  k  approximately 
"cancels"  the  k4  dependence  in  Eq.  (29));  (ii)  the  weak  dependence  on 
wind  velocity  at  angles  away  from  normal  (because  the  ocean  surface 
spectrum  W(k)  tends  to  become  saturated  at  moderate  wind  velocities 
for  those  larger  wave  numbers  which  control  the  microwave  scattering); 

(iii)  the  difference  between  the  horizontal  and  vertical  direct 
polarized  returns  (because  ahh  and  ayy  differ  in  the  off  normal  direction); 

(iv)  the  upwind/downwind/crosswind  ratios  (due  to  the  asymmettric  dis¬ 
tribution  of  the  slopes  of  the  large  scale  component);  (v)  the  power 
spectral  density  of  the  return  signal.  Other  features,  such  as  the 
return  from  near  grazing  angles,  and  certain  features  of  the  cross- 
polarized  return,  are,  however,  not  well  explained  by  the  composite 
model . 
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F. 


The  Lambert  Model 


Because  of  the  relative  simplicity  of  the  model  surfaces  described 
in  the  preceding  section,  they  have  attracted  the  attention  of 
theoreticians  to  an  extent  that  is  quite  out  of  proportion  to  their 
possible  application  to  terrestrial  surfaces.  Few  natural  surfaces 
can  be  conceptualized  as  continuous  Gaussian  randan  processes,  and  most 
must  be  treated  by  quite  different  kinds  of  models.  One  useful  group 
of  models  are  those  based  on  empirical  scattering  laws,  of  which  the 
best  known  is  the  "Lambert"  law,  with  bistatic  cross-section  given  by 

(47)  o-.  (e. ,  es,  $s)  =  r  cose^  cosbs 

o 

(e^)  =  r  cos  (backscatter) 

where  r  is  a  constant.  For  the  Lambert  surface,  the  scattered  energy 
is  considered  to  be  completely  depolarized  by  the  scattering  process. 

For  a  completely  reflecting  Lambert  surface,  ( i . e . ,  a  surface  with 
aTbedo  of  unity)  the  constant  r  has  a  value  of  r  =  2,  since  the  total 
scattered  power  must  be  summed  over  both  polarizations.  For  terrain 
surfaces  the  parameter  r  can  be  estimated  from  the  back-scattering 
cross-section. 

The  type  of  surface  likely  to  exhibit  Lambert  behavior  at  micro- 
wave  frequencies  is  one  which  is  highly  re-entrant  with  large  openings 
and  many  reflecting  planes,  for  example  a  surface  composed  of  large 
blocks  of  pumice  (see  Fig.  34).  Because  most  of  the  surface  structure 
is  large  in  terms  of  wavelength,  the  typical  Lambert  surface  should 
have  scattering  behavior  almost  independent  of  wavelength.  The 
emissive  properties  of  the  Lambert  surface  are  independent  of  viewing 
angle,  i.e.,  the  surface  acts  like  a  grey  body.  Backscattering  and 
brightness  temperature  for  a  typical  Lambert  surface  is  shown  in 
Figs.  50  and  61 . 
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There  are  a  number  of  other  empirical  scattering  laws  which 
have  been  proposed,  several  of  the  more  useful  of  which  are  illustrated 
in  Table  I.  This  tabie  also  gives  the  albedo  and  emissivity  for  each 
law.  Other  examples  of  empirical  scattering  laws,  mainly  of  theoretical 
interest,  are  given  in  Ref.  1. 


G.  Scattering  by  Individual  Elements 

(Lommell -See Tiger  Scattering! 

The  second  characteristic  type  of  diffusely  scattering  surface 
is  that  composed  of  a  number  of  individual  scatterers  disposed  in 
either  a  surface  or  a.  volume  distribution.  The  theoretical  approach 
to  this  type  of  scattering  is  via  the  scattering  functions  for  the 
individual  particles,  which  are  then  used  to  compute  the  total 
scattering  from  the  distribution  of  particles.  Various  levels  of  com¬ 
plexity  in  the  resulting  formulas  then  depend  on  the  extent  to  which 
multiple  scattering  effects  are  considered  and  the  probability 
distribution  assimed  for  the  particle  orientations.  Two  rather  detailed 
models  of  this  type  have  been  proposed  to  represent  the  scattering  from 
vegetated  surfaces  at  microwave  frequencies.  The  first  is  appropriate 
to  scattering  from  such  linear  vegetation  as  wheat,  oats,  sucian 
grass  etc,  where  the  individual  elements  may  be  approximated  by  a 
surface  distribution  of  thin  cylinders,  which  can  be  assumed  to  be 
much  less  than  a  wavelength  in  diameter.  In  this  case  the  bistatic 
scattering  cross  section  is  given  by  formulas  of  the  type[l ,/] 


(48)  aiJ(e1.es,*s) 


(AN)(AKZ) 


V1 


T  o 

28tt[3(ci/K)  +  3(cose..  +  coses) 

2  2 

+  (sin  e.  +  sin  e  -  2  sine,  sine  cos<j>  )] 
^  s  l  s  s 


79 


where  ec  =  relative  complex  dielectric  constant  of  the  cylinder  material 
A  =  cross-sectional  area  of  each  cylinder 
N  =  number  of  cylinders  per  unit  surface  area 
K  =  2tt/a. 
a  from  Eq.  50 

The  polarization  functions  f.  .  are  factors  of  order  unity  which  depend 

■  J 

on  the  scattering  and  incidence  angles,  and  the  parameter 
t  2/(i  +  Reuc)). 

For  material  of  large  dielectric  constant,  t  is  small,  and  the  polari¬ 
zation  factors  are  approximately 

fhn  *  <3  -  2  5i"\> 

2  2  2 
fyh  =  1  +  2  sin  e..  +  2  cos  0^  sin  $s 

(49)  fhv  =  1+2  sin2es  +  2  cos2es  sin2^. 

2  2  2 

f.„.  "3  •  2  cos  a  cos  0„  sin  if,, 
vv  i  s  s 

1 2(  si  no .  sino  -  cose,  coso  cos<f>  )  (sino.  sinep). 

1  S  1  S  'S  I  5 

The  parameters  (for  horizontal  incident  polarization)  and  ay 
(for  vertical  incident  polarization)  are  the  attenuation  factors  for 
the  incident  field  as  it  propagates  through  the  cylinders,  and  are 
given  by 

(50)  ah  =  (3K/16)  AN  (Ime  )  secoi  (1  +  3t2) 

av  =  (3K/16) AN  (Imec)  Seco.  [1  +  3t2  +  sin2©.  (1  -  t2)]. 
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Thus  1/a  ,  1/a^  give  estimates  of  the  depth  of  penetration  (optical 

depth)  of  microwave  radiation  into  such  vegetated  surfaces.  Because 

these  formulas  contain  a  number  of  parameters,  (A,  N,  cc  etc.)  which 

are  rather  difficult  to  estimate  for  actual  vegetated  surfaces,  one 

cannot  consider  them  to  give  quantitative  agreement  with  measured 

results.  They  can,  however,  be  useful  in  estimating  penetration 

depths,  and  surveying  the  effects  of  changes  in  such  parameters  as 

vegetation  moisture  content,  biomass  density  etc.,  on  the  radar  return. 

Figure  18  shows  typical  behavior  of  the  back. cattering  cross-sections 

for  this  model,  which  may  be  compared  with,  for  example,  Figs.  36,  43 

and  107.  A  number  of  calculations  based  on  these  formulas,  together 

with  higher  order  terms  in  the  expressions  for  f..,  are  given  in 

*  J 

references  1  and  7. 


The  second  type  of  vegetated  surface  for  which  fairly  detailed 
scattering  models  have  been  worked  out  is  that  for  which  the  leaves 
can  be  modelled  by  a  volume  distribution  of  their  circular  discs. 

Such  commercial  crops  as  soybeans,  cotton  etc.,  as  well  as  many 
deciduous  trees  and  shrubs  can  be  modelled  in  this  way.  Reference  20 
gives  a  detailed  derivation  of  the  bistatic  scattering  cross-sections 
for  these  volume  distributions  of  individual  scatterers,  viz 


(51)  0^(9,  ,6s.ts)  - 


cose.j  ccses 
cose^  +  coses 


where 


asJ  is  the  scattering  cross -sect; on  of  an  individual  particle, 
(i.e.,  leaf)  averaged  over  the  probability  distribution  of 
the  spatial  orientation  of  the  particle. 
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agc  is  the  total  scattering  cross  section  of  the  particle  for 
incident  angle  8^  and  incident  polarization  state  specified 
by  index  { i ) ,  averaged  over  all  spatial  orientations. 

is  the  absorption  cross  section  of  the  particle  for 

a 

incident  angle  8.  and  incident  polarization  state  specified 
by  the  index  {i),  averaged  over  all  spatial  orientations. 

Thus  a  detailed  study  of  the  scattering  process  of  the  individual 
particle,  and  a  knowledge  of  the  probability  distribution  of  the 
spatial  orientations  of  the  particles,  permits  one  to  estimate  the 
parameter  r  (see  Table  I)  used  in  the  empirical  form  of  the  Lonmel- 
Seeliger  Law.  Detailed  formulas  based  on  low  frequency  and  high 
frequency  approximations  to  scattering  by  thin  discs  are  given  in 
Ref.  20,  which  also  contain  a  number  of  calculations  for  the  back- 
scattering  cross-sections  for  the  model.  Figure  19  shows  typical 
calculations  for  backscattering  for  this  model,  which  can  be  compared 
with,  for  example,  Figs.  35  and  105.  Again,  because  of  the  many 
parameters  in  the  model,  one  should  not  expect  quantitative  agreement. 
However,  using  reasonable  values  for  leaf  thickness,  area  and  dielectric 
constant,  such  models  predict  qualitatively  the  dependency  of  the 
scattering  cross-sections  and  attenuation  lengths  on  frequency,  moisture 
constant  etc.  Figure  20  shows  estimates  of  the  attenuation  length 
(optical  depth)  of  vegetated  layers,  based  on  this  model.  The  models 
are  less  successful  in  predicting  the  depolarized  scattering,  presumably 
because  multiple  scattering  effects,  which  tend  to  increase  the  cross- 
polarized  power,  are  neglected. 
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"Slightly  rough"  surface  h  «  A.  Scattered  beamwidth  determined 
by  transmitter  beamwidth;  -  predominantly  coherent. 


Undulating  "rough"  surface  h  >>  A,  slopes  <<  1.  Scattered  beam- 
width  =  r.m.s.  surface  slope;  specular  (not  coherent). 


The  Lambert  type  of  "rough"  surface  The  Lommel-Seeliger  type  of 
with  large  reentrant  cavities.  "rough"  surface  with  many 

independent  scatterers. 


Fig.  14 
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Fig.  15.  Geometry  of  the  scattering  problem 
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Fig.  16.  Geometry  of  layered  surface. 
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10  LOO  tr„.  on*  10  LOG  0-h0h  MB)  t0  L(10  —  o  (M)<  (0  UQ0  -•  (dB) 


exbobchtul  *ror 


0  10  ?0  W  40  SO  *0  TO  «0  *0  0  10  ?0  50  40  SO  40  TO  *0  *0 

Ang«  Frtyr)  V*ff<oi ,  OoQ'm %  An q»  from  V**  »*CDi ,  Mgrm 


(*)«r«e.O  (4)  €r  •  8.0 

(b) 

The  small  scale  component  has  qaussian  surface  height  correlation 
coefficient  with  ^(r)  =  exo(-r2/«2)  (*  is  smal 1  scale  correlation 
length,  h  is  the  mean  square  height  of  the  small  scale  process. 
Solid  curve  kP.=l .  kh=0.2;  dotted  curve  k?-l,  klWJ.l.  The  complex 
dielectric  constant  55(1- j .55)  corresponds  roughly  to  sea  water 
at  X-band.  The  large  scale  component  has  rms  slope  =  s -j . 

Fig.  17.  Backscattering  cross-section  for  composite  model. 
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Fig.  20.  Optical  depth  vs.  frequency 
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V. 


EXPERIMENTAL  DATA 


A.  Introduction 

In  this  section,  a  large  number  of  data,*  both  radar  and  radio- 
metric,  are  presented  in  a  manner  chosen  to  illustrate  the  relations 
between  the  sensor  responses  (aQ,  y,  and  T^),  and  those  parameters  of 
the  surface  and  the  sensor  which  control  them. 

Among  the  parameters  which  can  affect  the  magnitude  of  the  sensor 
response  are: 

1 .  Surface  roughness 

2.  Complex  dielectric  constant 

3.  Frequency 

4.  Polarization 

5.  Angle  of  incidence 

For  the  most  part,  we  have  chosen  to  illustrate  the  effects  of  surface 
roughness,  dielectric  constant,  frequency,  ana  polarization,  on  the 
sensor  responses  by  present  ng  the  data  as  a  function  of  angle  of 
incidence  (from  normal).  Some  of  the  large  angle  data,  however,  are 
plotced  as  a  function  of  the  grazing  angle.  Because  of  the  interrela¬ 
tions  between  the  parameters  above  (e.g.,  polarization  and  dielectric 
constant  effects  for  smooth  surfaces)  there  is  a  certain  amount  of 
redundancy  in  the  data  presented. 

The  data  are,  of  course,  by  no  means  complete,  brt  it  is  believed 
they  represent  most  of  the  categories  of  surface  for  which  data  are 
available.  While  a  large  part  of  the  terrain  data  represent  only  the 
measurement  of  a  single  investigator,  there  are  a  number  of  surfaces 
(particularly  the  sea  and  the  planets)  which  have  been  studied  inde¬ 
pendently  by  several  investigators.  In  general  it  may  be  assumed  that 

*Most  of  the  backscattering  data  in  this  section  display  the  parameter 
y  =  o0  secej  rather  than  o0. 


the  calibration  procedures  of  any  one  investigator  will  be  consistent 
from  run  tc  run,  but  that  there  may  be  systematic  differences  between 
the  absolute  levels  measured  with  different  experimental  configurat.ons . 

The  most  evicent  case  where  measurements  of  ostensibly  the  same  surface 
give  significantly  different  absolute  levels  occurs  for  tne  sea  surface. 

In  this  case,  however,  because  the  state  of  the  surface  can 
change  rapidly,  and  is  often  specified  indirectly  (by  using  the  wind 
velocity  rather  than  the  surface  roughness  as  the  descriptor)  the 
observed  discrepancies  may  be  due  to  the  fact  that  different  surfaces 
were  measured,  rather  than  to  systematic  calibration  errors.  A  further 
note  of  caution  should  be  inserted  about  cross-polarized  return.  There 
is  a  tendency,  especially  over  smooth  surfaces  near  normal  incidence, 
for  the  measured  depolarized  return  to  equal  the  polarization  isolation 
of  the  antenna  used  in  the  experiment.  In  these  cases,  the  measured 
return  must  be  considered  only  as  an  upper  limit  to  the  actual  de¬ 
polarized  power. 

The  data  to  be  presented  ar^  organized  in  the  following  manner. 

First,  the  general  characteristics  of  scattering  from  terrain  surfaces 
are  illustrated  via  the  bis:atic  scattering  pattern.  Although  complt*'1 
hemispherical  patterns  have  not  been  measured  at  microwave  frequencies, 
certain  cuts  of  the  bistatic  cross  section  (per  unit  surface  area) 
aQ( 0i ,  efi>  <f>) ,  have  been  measured[l]  at  10  GHz  for  a  number  of  well- 
defined  surfaces,  including  smooth  and  rough  sand,  and  smoothed  loam, 
dry  grass,  and  soybeans.  The  measurements  cover  a  wide  range  of 
incidence,  scattering,  and  azimuth  angles,  for  vertical,  horizontal, 
and  crossed  linear  polarizations.  Selected  curves[l,2]  are  presented 
in  Section  B. 

Ground-based  measurements  of  the  normalized  radar  backscattering 
c^oss  section  (per  unit  "projected"  area),  y{Q^),  for  a  number  of 
geological  and  vegetation-covered  surfaces,  have  been  made  by  this 
facility[3,4,5,6]  in  recent  years,  at  1.8  GHz,  10  GHz,  15  GHz,  and  35  GHz. 
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A  number  of  these  results  are  presented  in  Section  C  to  illustrate 
the  relation  between  the  backscatterirg  cross  section  and  the  parameters 
of  the  surface  and  the  sensor.  It  should  be  noted  that  the  absolute 
levels  of  our  earlier  data[3]  are  consistently  lower  than  our  more  recent 
data[5,6]  for  similar  surfaces.  This  difference  may  be  due  to  a  otic 

error  in  the  earlier  calibration  technique. 

Measurements  of  the  brightness  temperature,  T^(e^),  of  a  number  of 
terrestrial  surfaces  of  geological  and  agricultural  interest  have  been 
obtained, [7, 8]  over  the  past  several  years,  with  the  Ohio  State  University 
truck-mounted  radar- radiometer  facility,  at  10  GHz  and  35  GHz.  Curves, 
dra.  n  from  these  data  together  with  a  number  of  measurements  due  to 
other  investigators,  are  presented  in  Section  D  to  illustrate  the 
dependence  of  the  radiometric  temperature  upon  the  surface  roughness, 
dielectric  constant,  frequency,  and  polarization. 

In  Section  E,  representative  measurements  of  radar  bi static  and 
backscattering  cross  sections  and  radiometer  brightness  temperatures, 
are  provided  for  several  special  surface  categories  including  geo¬ 
logical  surfaces,  vegetated  surfaces,  developed  land  areas  (cities), 
and  sea  surfaces.  The  data  compiled  in  these  curves  were  gathered 
from  many  sources  and  should  provide  a  convenient  survey  of  microwave 
radar  responses. 

B.  Bistatic  Scattering  Data 

In  this  section,  a  number  of  data  curves  are  presented,  in 
Figs.  22-30,  to  illustrate  the  general  features  of  the  bistatic  pattern 
of  tyoical  surfaces,  and  to  show  the  relation  of  surface  roughness, 
dielectric  constant,  and  antenna  polarization  tc  the  bistatic 
scattering  patterns.  Other  bi static  data  are  given  in  Figs.  75-89. 

A  complete  description  of  the  expe?'i  mental  processes,  and  a  number  of 
additional  data  curves  are  given  in  Reference  1.  In  order  to  separate 
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more  clearly  the  several  sets  of  data  in  each  figure,  the  measured 
points  have  been  arbitrarily  connected  by  lines,  but  it  should  not 
be  assumed  that  measurements  at  intermediate  points  would  lie  on 
the  indicated  curves.  Photographs  are  shown  in  Fig.  21,  of  the  actual 
measured  surfaces  which  include:  sand  which  had  been  smoothed  with 
a  board  so  that  surface  inequalities  were  less  than  1/4  wavelength 
("smooth  sand");  sand  which  had  been  moulded  into  hillocks  with  a 
roughness  scale  of  several  wavelengths  ("rough  sand");  loam  which  had 
been  raked  smooth  so  that  it  exhibited  irregularities  of  about  1/4 
wavelength  ("loam");  dry  grass  several  inches  long  and  partly  matted 
down;  loam  with  a  stubble  covering;  and  full  grown  soybean  plants. 

1 .  Effects  of  surface  roughness 

Figure  22,  which  exhibits  the  bi static  pattern  in  the  plane  of 
incidence  ($  =  0)  for  fixed  values  of  incidence  angle  illustrate  the 
characteristic  behavior  of  a  smooth  surface,  i.e.,  one  for  which  the 
roughness  is  less  than  1/4  wavelength.  A  pronounced  forward  lobe,  in 
approximately  the  specular  direction  is  visible  in  each  case  with  the 
smoother  surface  (smooth  sand)  having  the  larger  return.  Fig.  23 
illustrates  the  characteristic  difference  between  smooth  and  rough 
surfaces;  for  the  latter  the  forward  lobe  is  much  less  prominent 
and  the  scattering  more  nearly  isotropic.  Some  trace  of  the  forward 
lobe  is  still  visible  for  the  "dry  grass"  surface,  showing  penetration 
through  the  grass  and  reflection  from  the  earth  beneath.  Deep 
surfaces  of  many  independent  scattererc.  which  should  give  an 
essentially  isotropic  scattering  pattern,  arb  illustrated 
by  the  soybean  measurements  of  Figs.  24  and  28.  A  second  type  of  rough 
surface  (the  "rough  sand"  surface  composed  of  a  distribution  of  large 
hemispheres)  also  exhibits  a  nearly  isotropic  scattering  pattern  ^s 
can  be  seen  in  Figs.  79(a)  and  79(b).  Here  the  scattering  is  nearly 
isotropic  for  the  same  reason  that  the  scattering  from  a  sphere  is 
isotropic,  although  multiple  scattering  effects  are  clearly  of  some 
importance  for  this  surface. 
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Scattering  behavior  in  the  azimuth  cone  is  illustrated  in  Fig.  24, 
(with  8^  =  6^  fixed;  <j>  measured  from  forward  direction).  Again  a 
pronounced  forward  lobe  is  evident  for  the  smooth  surfaces.  It  is  inter¬ 
esting  that  in  most  cases  the  minimum  scattering  occurs  when  <j>  =  90°, 
with  relatively  larger  returns  in  the  backscattering  direction. 


2.  Effects  of  complex  dielectric  constant 

The  best  method  to  illustrate  the  effects  of  dielectric  constant 
on  the  bi static  return  is  to  examine  the  cross-sections  for  smooth 
surfaces  (for  example,  see  Fig.  22).  For  those  situations  where  the 
scattering  exhibits  a  distinct  lobe  in  the  specular  direction,  it  is 
possible  to  calculate  an  equivalent  "Fresnel"  reflection  coefficient 
from  the  measured  cross  section  by  the  expression[l] 


where  oQ  is  the  value  measured  et  the  specular  angle  (e  =  0. ,  <J>  -  0) , 
because  this  is  analogous  to  the  case  of  a  perfectly  smooth  surface. 

(i  ire  I  ^  tt/(B  cos  9. )»  where  is  an  approximation  to  the  one¬ 

way  power  pattern  of  identical  transmitting  and  receiving  antennas). 

The  reflection  coefficients  for  smooth  sand  and  loam  at  the  two 
antenna  polarizations  were  calculated  in  this  manner  and  plotted  in 
Fig.  25.  The  solid  curves  are  calculated  values  of  |R|  assuming 
perfectly  smooth,  lossless  surfaces  with  relative  dielectric  constants 
2  and  3,  and  the  points  clustered  around  these  curves  are  the  experi¬ 
mentally  measured  values.  The  contrast  between  the  vertically 
polarized  and  horizontally  polarized  returns  clearly  illustrates  that 
tha  Brewster  effect  is  present,  and  suggests  that  the  dry  sand  and 
dry  loam  surfaces  have  an  effective  dielectric  constant  of  about  3. 
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Figure  26  shows  results  of  reflection  coefficient  measurements 
made  by  Sherwood  and  Ginzton,[9]  indicating  that  for  the  purpose  of 
revealing  the  Brewster  angle  of  a  surface,  (and  therefore,  the  relative 
dielectric  constant),  a  number  of  closely  spaced  measurements  about 
the  estimated  Brewster  angle  is  required  for  satisfactory  results. 


3.  Effects  of  polarization 

Effects  of  polarization  are  illustrated  by  Figs.  27-28  (azimuthal) 
and  29-30  (plane  of  incidence).  (See  also  Fig.  22).  As  might  be  expected 
from  the  behavior  of  the  Fresnel  coefficients,  the  horizontal  com¬ 
ponent  is  significantly  larger  than  the  vertical  component  for  smooth 
surfaces  in  the  specular  direction,  particularly  near  the  Brewster 
angle.  However,  at  angles  near  <j>  =  90°  the  situation  is  reversed,  with 
vertical  predominating.  For  vegetated  surfaces  (e.g.,  Fig.  30)  there  is 
little  difference  between  the  two  polarizations.  Of  some  interest  are 
the  cross  polarized  returns,  which  for  smooth  surfaces,  are  small 
(relative  to  direct  polarization)  in  the  forward  lobe,  as  would  be 
expected,  but  are  comparable  to  the  direct  returns  at  angles  near 
4)  =  90°  (Fig.  29).  For  vegetated  surfaces,  the  cross-polarized  return 
tends  to  follow  the  direct  return  (Fig.  30)  but  at  a  somewhat  lower 
power  level  (see  also  Fig.  92).  (Note  in  Fig.  30,  the  last  few  data 
points  were  made  as  snow  began  to  fall.) 


C.  Backseat tering  Data 

Curves  are  presented  here  to  illustrate  the  dependence  of  the 
radar  backscattering  cross  section,  y,  upon  such  parameters  as  surface 
roughness,  dielectric  constant,  frequency,  and  polarization.  These 
data  have  been  drawn  from  measurements  by  several  investigators[3,4,5,6] 
at  this  facility.  A  complete  description  of  the  Ohio  State  University 
truck-mounted  radar-radiometer  facility  and  its  operation,  is  given 
in  Reference  10. 
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i .  Effects  of  surface  roughness 


Figure  31  illustrates  the  characteristic  behavior  of  tne  radar 
backscattering  cross  section,  y,  for  "smooth"  surfaces,  ',e,,  the 
return  decreases  monotonically  as  the  angle  of  incident  increases 
and  the  magnitude  of  the  horizontally  polarized  return  is  less  than 
that  of  the  vertically  polarized  as  predicted  by  the  perturbation 
model.  For  such  slightly  rough  surfaces,  the  magnitude  of  the  return 
provides  an  estimate  of  the  rms  roughness,  and  the  difference  between 
the  vertical  and  horizontal  return  provides  an  estimate  of  the  dielectric 
constant.  The  glacially  polished  limestone  surface  is  level,  with  small 
scale  roughness  <  1  mm. 

One  form  of  transition  from  a  "smooth"  to  a  "rough"  surface  is 
exhibited  by  the  radar  backscattering  cross  sections  of  rubble  type 
surfaces,  at  X-band,  shown  in  Fig.  32.  Here  the  surfaces  were  similar 
in  character,  i.e.,  sand,  pea  gravel,  and  crushed  limestone  (with  a 
diameter  of  3-5  cm  for  the  individual  rocks;,  but  differed  in  roughness, 
with  the  sand  satisfying  the  Rayleigh  criteria  for  a  "smooth"  surface  at 
3  cm  wavelength,  while  the  crushed  stone  was  "rough".  Again  the  backscat¬ 
tering  increases  with  roughness,  as  does  the  Depolarized  return  relative  to 
the  direct,  until  the  Rayleigh  criterion  is  exceeded.  For  the  crushed 
stone,  the  significant  multiple  scattering  produces  a  cross-polarized 
return  almost  equal  to  the  direct.  Figure  33  illustrates  the  limiting 
form  for  this  kind  of  roughness,  that  is,  the  discontinuously  rough 
surface  having  large  re-entrant  cavities.  Such  surfaces  exhibit  a  back- 
scattering  cross  section  almost  independent  of  frequency,  and  behave 
more  or  less  like  a  Lambert  Law  surface  (  y(e.)  a  cos  e.).  A  photograph 
of  the  actual  pumice  surface,  at  Mono  Crater,  is  shown  in  Fig.  34.  The 
range  of  diameters  of  the  individual  boulders  is  30  to  100  cm. 

The  other  type  of  discontinuously  rough  surface  (the  Lornmel- 
Seeliger  type  of  scattering  model;  y  =  constant)  is  illustrated  by 
vegetation-covered  terrain.  Figures  35  and  36  exhibit  the  backscattering 
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cross  sections  at  X-band,  for  vertical  polarization,  for  two  distinct 
types  of  vegetation:  those  in  which  the  leaf  structure  is  similar  to 
a  flat  disc  (soybeans),  and  those  in  which  the  leaves  are  essentially 
thin  cylinders  (oats).  For  vegetated  surfaces,  the  magnitude  of  y  is 
seen  to  be  almost  independent  of  angle  of  incidence,  as  one  would 
expect  from  the  model.  It  is  also  not  strongly  dependent  on  crop  type, 
although  wheat  tends  generally  to  have  a  lower  y,  and  sugar  beet  a 
significantly  higher  y  than  the  crops  illustrated  here, 

2.  Effects  of  complex  dielectric  constant 

The  effect  of  the  dielectric  constant,  for  two  slightly  rough 
surfaces  (the  limestone  and  asphalt  surfaces  of  Fig.  31),  is  illu¬ 
strated  in  Figs.  37  and  38,  showing  the  difference  between  the  vertical 
and  horizontal  backscatter.  The  solid  curves  are  computed  values  (using 
Eqs.  27,29  of  Chapter  IV)  of  the  difference  between  y„y  and  yhh,  independent 
of  the  surface  roughness  parameter.  Thus,  by  measuring  the  backscattering 
for  vertical  and  horizontal  polarizations,  as  a  function  of  angle  of 
incidence,  it  is  possible  t.o  estimate  the  dielectric  constant  for  a 
slightly  rough  surface. 

For  surfaces  of  the  rubble  type  witn  irregularities  comparable 
to  the  Rayleigh  criteria,  the  effect  of  roughness  and  dielectric 
constant  on  the  radar  backscatter  is  clarified  if  the  data  at  a  number 
of  angles  and  wavelengths  is  plotted  in  terms  of  a  normalized  roughness 
parameter,  5  =  (D  cos  e.)/x,  where  D  is  the  diameter  of  the  individual 
particles.  For  example,  if  the  return  from  gravel  at  3  wavelengths 
in  Fig.  39,  and  that  of  similar  measurements  on  crushec  rock 
(D  :>  3-5  cm)  and  sand  are  plotted  versus  the  normalized  roughness 
parameter,  5,  the  result  is  Fig.  40.  It  is  seen  for  £  <  1/4,  the  return 
falls  off  roughly  as  1/X‘,  whereas  for  5  >  1/4  the  return  is  a  constant, 
more  or  less  independent  of  roughness  (corresponding  to  "optical" 
scattering).  The  general  level  of  the  curve  would  move  up  or  down 
for  materials  of  different  dielectric  constant.  Thus,  by  measuring 
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the  backscattering  from  a  surface  of  this  type  at  a  number  of  angles 
and  wavelengths,  it  should  be  possible  to  estimate  the  roughness  by 
locating  the  break-point  (5  ±0.3)  and  to  estimate  the  dielectric  con¬ 
stant  by  measuring  the  optical  level. 

The  dependence  of  the  backscattering  cross  section  upon  the 
dielectric  constant  is  further  illustrated  by  data[3]  on  seasonal 
changes  of  grass  at  three  frequencies  in  Figs.  41-43.  At  X-band 
(Fig.  43)  the  variation  in  the  magnitude  of  y ,  which  may  be  as  large 
as  8  dB  (regardless  of  the  variation  in  height  of  the  vegetation),  is 
presumably  due  to  changes  in  water  content,  with  the  resulting  change 
in  the  effective  dielectric  constant  of  the  grass  blades.  Other  vege¬ 
tated  surfaces,  however,  (e.g.,  soybeans  Fig.  35,  oats  Fig.  36,  grass 
at  ku  and  Ka  band  Figs.  41  and  43)  show  little  change  in  y  over  most  of 
the  growing  season. 

Figures  44  and  45  illustrate  the  effects  of  rain  on  smooth 
(asphalt  road)  and  rough  (grass)  surfaces,  at  K  -band.  Water  on  a 

a 

relatively  smooth  surface  tends  to  lower  the  return  since  the  surface 
appears  to  be  smoother,  and  this  effect  overrides  the  tendency  of 
the  water  to  increase  the  effective  dielectric  constant.  On  vege¬ 
tation,  however,  the  return  may  increase  because  the  water  drops  on 
the  leaves  have  effectively  increased  both  the  leaf  thickness  and 
its  effective  dielectric  constant. 

3 .  Effects  of  frequency 

Figure  46  exhibits  the  frequency  dependence  of  radar  backscat¬ 
tering  for  a  typical  smooth  (i.e.,  "slightly  rough")  asphalt  surface. 

As  the  frequency  decreases,  the  magnitude  of  y  should  decrease  at 

2 

least  as  rapidly  as  f  .  The  particular  surface  shown  has  a  rms 
roughness  of  approximately  0.3  mm  and  a  relative  dielectric  constant 
of  er  =  4.3  +  jO.l . 
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Figures  47-49  illustrate  the  frequency  dependence  of  the  radar 
backscattering  for  three  rubble-type  surfaces  (sand,  gravel,  and 
crushed  stone).  The  gravel  surface  (Fig.  48)  shows  the  kind  of  in¬ 
crease  in  y  with  frequency  expected  for  a  slightly  rough  surface. 
However,  it  is  interesting  to  note  that  for  the  case  of  crushed 
stone  surface  the  largest  return  was  obtained  at  X-band,  i.e.,  at 
the  lowest  frequency.  It  is  possible  that  either  resonance  scattering, 
or  a  significant  decrease  in  dielectric  constant  with  frequency, 
could  cause  this  behavior. 

The  characteristic  behavior  of  the  frequency  dependence  of  v 
for  the  two  types  of  discontinuously  rough  surfaces,  large  blocks 
(pumice)  and  vegetation  (oats),  is  shown  in  Figs.  50  and  51.  As 
previously  noted,  the  large  blocks  of  pumice  exhibit  a  radar  return 
almost  independent  of  angle  and  frequency  (i.e.,  "optical"  behavior). 

For  the  oats  surface,  however,  a  significant  difference  in  the  magni¬ 
tude  of  y  exists  between  the  1.8  GHz  and  10  GHz  return.  This 
frequency  behavior  would  be  expected  on  the  basis  of  theoretical 
models,  since  at  1.8  GHz  the  vegetation  layer  may  be  less  than  one 
optic?1  Jepth  thick,  and  the  scattering  from  an  individual  cylinder 
decreases  rapidly  as  wavelength  increases. 

4.  Effects  of  polarization 

Figure  52  illustrates  the  polarization  dependence  of  the  radar 
backscattering  cross  section,  y>  for  a  "smooth"  surface  (glacially 
po!ished  limestone  located  near  Marblehead,  Ohio);  i.e.,  the  hori¬ 
zontally-polarized  backscatter  (y^)  1S  ^ess  than  that  cf  the  vertically 
polarized  (yyy)»  and  the  cross-polarized  components  (y^  and  y^y)  are 
at  least  10  dB  less  than  the  .:i rect-oolarized  components. 

The  polarization  dependence  (for  both  linear  and  circular 
polarization)  of  the  radar  backscattering,  at  X-band,  for  the  three 
rubble-like-surfaces  described  previously,  is  shown  in  Figs.  53-55. 
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It  can  be  seen  that  the  vertically  polarized  return  is  larger  than  the 
horizontal  for  the  quas^-smooth  sand  surface,  and  the  cross  polarized 
return  is  at  least  10  to  15  dB  down.  However,  as  before,  the  depolar¬ 
ized  component  increases  as  the  surface  becomes  rough  in  terms  of 
wavelength. 

Figures  33  and  56  illustrate  the  polarization  behavior  of  the 
radar  backscatter  for  discontinuously  rough  surfaces  (pumice  and  green 
oats),  i.e.,  y  almost  independent  of  polarization  for  the  direct 
returns,  and  cross-polarized  returns  (for  the  oats  surface  in  particular) 
which  are  only  slightly  lower  than  the  direct  returns.  This  is  in 
contrast  with  the  relatively  small  cross-polarized  returns  for  a  slightly 
rough  surface  (e.g.,  Fig.  52),  but  is  characteristic  of  surfaces  in  which 
multiple  scattering  plays  a  significant  role. 

D.  Brightness  Temperature  Data 

Measured  microwave  brightness  temperatures  are  presented  in  this 
section  to  show  their  dependence  on  the  parameters  of  the  surface 
and  the  sensor.  The  measurements  are  primarily  those  obtained  by 
this  facility, [7, 8]  although  a  number  of  other  data, [11]  are  included. 

In  the  case  of  the  data  of  References  7  and  8,  the  integral  equation 
relating  brightness  temperature  and  antenna  temperature  has  been 
inverted  by  a  bootstrap  method,  so  that  the  data  plotted  are  estimates 
of  brightness  temperature.  In  most  other  cases,  raw  antenna  temp¬ 
erature  is  plotted. 

1 .  Effects  of  surface  roughness 

The  measured  brightness  temperatures  for  smootn  asphalt,  at 
10  and  35  GHz,  are  shown  in  Figs.  57  and  58,  respectively.  The  effect 
of  the  near  zero  reflection  co*efficient  (and  thus,  emissivity  very  near 
unity)  at  the  Brewster  angle  for  vertical  polarizati.  ,  is  prominent 
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for  both  frequencies.  A  fairly  good  estimate  of  the  dielectric 
constant  of  the  surface  may  be  obtained  from  the  characteristic  shape 
of  the  curves. 

Figure  59[11]  illustrates  the  transition  from  "smooth"  to  ’’rough" 
for  very  homoqeneous,  rubble-type  surfaces,  composed  of  different 
size  fractions  of  limestone,  at  13.5  GHz.  The  gravel  surface  exhibits 
the  charactei istic  of  a  slightly  rough  surface.  The  cobble  surface 
shows  some  of  the  characteristics  of  a  "smooth"  surface  but  the  Brewster 
angle  effect  is  decreased.  The  cobbles  included  size  fractions  passing 
a  10  cm  by  15  cm  screen.  The  boulder  surface  (diameters  of  several  feet) 
is  characteristically  rough;  i.e.,  the  brightness  temperature  is 
approximately  independent  of  angle  of  incidence,  i.e.,  approximately 
Lambertian . 

Figures  60-62  present  a  set  of  X-band  measurements [7]  made  at 
Mono  Craters,  California  on  a  sequence  of  lavas  of  essentially  the  same 
chemical  composition  and  structural  roughness  but  different  physical 
properties,  ranging  from  light  pumice  to  obsidian  (see  photographs  in 
Fig.  63).  The  brightness  temperatures  again  exhibit  the  characteristic 
behavior  for  a  Lambert  Law  surface,  i.e.,  almost  independent  of  angle  of 
incidence.  Here  the  denser  obsidian  has  a  higher  dielectric  constant, 
a  higher  reflection  coefficient,  and  thus,  a  lower  brightness  temperature. 
Since  the  radar  return  for  the  surfaces  also  show  the  characteristic  of 
a  Lambert  surface  (e.g.,  see  Fig.  33),  the  brightness  temperature  can 
be  related  directly  to  the  radar  return,  and  estimated  from  it. 

Measurements  of  X-band  brightness  temperatures  of  cultivated 
farmland,  in  particular  plots  of  alfalfa,  wheat,  and  oats,,  are  shown 
in  Figs.  64-66.  The  vegetated  surfaces  behave  as  rough,  diffusely 
scattering  surfaces,  as  evidenced  by  the  lack  of  angular  iependence 
for  the  brightness  temperatures. 
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2.  Effects  of  complex  dielectric  constant 

Figure  67  illustrates  the  effect  of  the  dielectric  constant 
on  measured  X-band  brightness  temperatures  (vertical  polarizaticn)  for 
asphalt,  glacially  polished  limestone,  and  a  smooth  coal  bed  exposed 
by  strip  mining  operations.  The  three  curves  represent  brightness 
temperatures  computed  from  the  Fresnel  reflection  coefficients, 
using  values  of  dielectric  constant  appropriate  to  the  actual  surfaces. 
It  is  clear  that  not  only  is  there  good  agreement  between  measured 
and  computed  temperatures,  but  that  such  measurements  offer  a  means  for 
remotely  determining  the  value  of  dieTectric  constant  for  exposed, 
smooth  surfaces. 

The  dependence  of  the  br  ghtness  temperature  of  soil  on  the 
soil  moisture  content  (and  hence  on  dielectric  constant)  is  exhibited 
by  13.4  GHz  data, [11]  in  Fig.  68,  which  were  obtained  on  Harper  Lake, 
a  playa  deposit  in  the  Mojave  Desert,  California.  It  is  noted  from 
these  data  that,  as  the  moisture  content  increases,  the  observed 
temperatures  decrease,  and  the  difference  between  the  vertical  and 
horizontal  temperatures  increase.  This  is  in  agreement  with  the  fact 
that  the  effective  dielectric  constant,  and  thus  the  reflection 
coefficient  of  the  surface  increase^  '!th  the  moisture  content. 

Figure  69  illustrates  the  effect  of  dielectric  constant  on  the 
measured  X-band  brightness  temperatures  for  three  pumice  surfaces 
(see  also  Figs.  60-62).  As  noted  previously, [12]  the  differences  in 
temperature  for  these  surfaces,  ere  directly  related  to  differences 
in  dielectric  constant,  and  indirectly  to  differences  in  density. [12] 

The  brightness  temperatures,  at  X-band,  tor  irrigated  and  non- 
irrigated  soybeans,  are  shown  in  Fig.  70.  For  these  surfaces,  at 
least  part  of  the  difference  in  brightness  temperature  must  be 
ascribed  to  a  difference  in  moisture  content  of  the  two  categories 
of  vegetation. 


3.  Effects  of  frequency 


The  frequency  dependence  of  for  a  "smooth"  surface,  asphalt 
at  X-baid  and  K  -band,  is  shown  in  Fig.  71.  The  Brewster  angle  effect 

a 

is  quite  prominent  at  both  frequencies,  indicating  a  surface  rough¬ 
ness  less  than  the  Ray^igh  criteria.  This  is  in  agreement  with  the 
known  rms  roughness  =  0.3  m.m.  The  surface  is  cooler  at  the  lower 
frequency  a;  one  would  expect  from  the  larger  dielectric  constant  at 
X-band. 

The  transition  from  smooth  to  rough  is  illustrated  in  Fig.  72, 
which  shows  the  brightness  temperatures  for  a  slightly  rough 
lapilli  surface  (photograph  in  Fig.  73),  at  10  GHz  and  35  GHz.  The 
brightness  temperature  at  10  GHz  is  similar  to  that  of  a  smooth  surface 
(for  instance,  Fig.  57).  At  35  GHz,  however,  the  return  shows  a 
character  intermediate  between  a  "smooth"  surface  and  a  "rough"  surface 
(c.f.,  Fig.  61). 

In  Fig.  74,  brightness  temperatures  are  shown  for  a  rough 
surface  (dark  pumice),  at  10  GHz  and  35  GHz.  In  contrast  to  the  lapilli 
surface  (Fig.  72),  the  brightness  temperatures  of  the  pumice  surface 
are  approximately  independent  of  angle  of  incidence  at  both  frequencies, 
as  expected  for  a  Lambert  surface;  furthermore,  if  the  brightness 
temperatures  are  scaled  according  to  the  actual  surface  temperatures, 
the  corresponding  emission  ties  are  practically  identical,  i.e., 
independent  of  frequency. 

4.  Effects  of  polarization 

The  dependence  of  the  microwave  brightness  temperature  on  the 
antenna  polarization  can  be  noted  from  data  in  the  previous  sections. 

An  excellent  example  of  the  polarization  behavior  of  a  "smooth" 
surface  (asphalt)  was  shown  in  Fig.  57,  i.e.,  a  prominent  Brewster 
angle  effect  for  the  vertical  component  and  a  continuous  decrease 
with  nadir  angle  for  the  horizontal  component. 
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For  rough,  diffusely  scattering  surfaces,  however,  the  microwave 
brightness  temperatures  are  almost  independent  of  polarization  This 
lack  of  polarization  dependence  is  clearly  evidenced  by  the  twa  typical 
rough  surfaces:  dark  pumice  (Fig.  61)  and  oats  (Fig.  66). 

E.  Special  Surface  Categories 

1 .  Bistatic  scattering  data 

In  this  section,  X-band  measurements [1]  of  aQ  are  presented  for  the 
following  land  surfaces:  bare  soil  (smooth  and  rough  sand,  and  loam), 
and  vegetation  (stubble,  dry  grass,  and  soybeans).  A  description  and 
photographs  (Fig.  21)  of  the  actual  surfaces  measured,  were  presented  in 
Section  B.  The  measurements  cover  a  wide  range  of  incidence,  scattering, 
and  azimuth  angles,  for  vertical,  horizontal,  and  crossed  linear  polari- 
zati on . 

C-band  measurements,  by  Pidgeon,[13]  of  0Q  fer  the  sea  are  also 
included.  These  measurements  cover  transmitter  depression  angles 
between  0.2°  and  3°  below  the  horizontal,  and  receiver  depression 
angles  between  10°  and  90°,  for  vertical  and  crossed  linear  polarizations. 

a)  Bare  soil  surfaces. 

The  bistatic  scattering  patterns  for  smooth  sand  are  shown 
in  Figs.  75  (plane  of  incidence)  and  76  (azimuth).  Similar 
measurements  for  smoothed  loam  are  shown  in  Figs.  77  and  78.  These 
measurements  again  illustrate  the  pronounced,  polarization  depen¬ 
dent  forward  lobe  for  smooth  surfaces. 

Figure  79  shows  the  bi static  pattern  in  the  plane  of 
incidence  for  rough  sand.  For  this  surface  the  scattering  is 
more  nearly  isotropic,  corresponding  to  the  behavior  expected 
from  a  distribution  of  hemispheres  on  a  plane. 
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b)  Vegetated  surfaces 


The  scattering  behavior  in  the  plane  of  incidence  for  three 
vegetated  surfaces  (stubble,  dry  grass,  and  soybeans)  is  shown 
in  Figs.  80-82.  For  the  soybean  surface,  the  azimuth  scattering 
is  shown  in  Fig.  83,  for  =  60°  and  80°.  For  the  vege¬ 

tated  surfaces,  there  is  little  difference  between  the  two  direct- 
polarized  returns,  and  the  scattering  pattern  tends  to  be  nearly 
isotropic. 

c)  Sea  surface  (Data  from  Ref.  13) 

Figures  84-86  show  bistatic  cross  sections  of  the  sea  at 
C-band  for  vertically-polarized  transmitter  and  receiver,  for 
three  values  of  sea  state.  The  receiving  antenna  is  located  in 
the  plane  of  incidence.  Figures  87-89  present  the  bi static  cros^. 
sections  for  a  vertically-polarized  transmitter  and  a  horizontally- 

polarized  receiver.  The  results  indicate  values  of  a  at  sea 

c 

state  3  are  approximately  10  dB  larger  than  corresponding  values 
for  sea  state  1.  The  cross-polarized  return  is  10  to  15  dB  less 
than  the  polarized  return  for  transmitter  depression  angles  less 
than  1°,  but  only  5  to  8  dB  less  at  depression  angles  near  3°. 

F.  Backscattering  Data 

Measurements  of  y  (or  c  )  are  presented  in  this  section,  for 
the  following  surface  categories:  geological,  vegetated,  developed 
land  areas  (cities),  and  sea.  A  limited  number  of  near-grazing  angle 
measurements  are  included  separately. 

1 .  Geological  surfaces 

The  backscattering  cross  sections  for  several  slightly  rough 
surfaces,  including  smooth  concrete,  concrete,  smooth  asphalt,  rough 
asphalt,  and  gravel,  are  shown  in  Fig.  90,  at  K  -,  K  -,  and  X-band. 
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It  is  evident  that  the  dependence  upon  angle  decreases  and  the  absolute 
level  increases  as  the  roughness  increases  in  terms  of  wavelength.  The 
mean  square  roughnesses  were  determined  for  four  of  the  actual  surfaces: 

"smooth  concrete" 

"concrete" 

"smooth  asphalt" 

'rough  asphalt" 

Figures  91-94  present  measurements  of  the  backscattering  cross 
section  for  a  number  of  terrestrial  surfaces  of  geological  interest. 

In  particular,  airborne  data  for  "normal"  desert  surfaces,  at  X-band 
and  horizontal  polarization,  are  presented  in  Fig.  91.  These  data 
were  obtained  by  the  Naval  Research  Laboratory[14 ,1 5]  and  by  the 
Goodyear  Aircraft  Corporation. [16]  Measurements  of  aQ  for  the  "Arizona 
Desert"  surface,  at  428  MHz,  1.228  GHz,  4.455  GHz,  and  8.91  GHz, 
for  horizontal,  vertical,  and  crossed  linear  polarizations,  may  be 
found  in  Reference  15. 

In  Fig.  92,  airborne  measurements  of  y,  for  a  "Pisgah  Lava 
Flow"  surface,  are  compared  with  ground-based  data[5]  for  a  similar 
surface.  It  is  noted  that  agreement  between  airborne  and  ground- 
based  measurements  is  quite  good. 

Figures  93  and  94  exhibit  the  backscattering  cross  sections  for 
Bristol  Dry  Lake  (Amboy,  California)  and  Lavic  Dry  Lake  (Pisgah 
Crater  in  the  Mojave  Desert,  California),  respectively.  A  photograph 
of  the  Lavic  Dry  Lake  playa  (shown  in  Fig.  95)  shows  the  surface  to  be 
almost  perfectly  flat,  having  a  hard,  dry  crust. 

Backscattering  cross  sections,  y(e..),  for  the  Moon, [18]  Venus[19] 
and  Earth, [17,20]  are  shown  in  Fig.  96.  The  Moon  (X  =  3.8  cm,  23  cm, 

68  cm)  and  Venus  (x  =  12.5  cm)  radar  transmissions  were  circularly 


Z2  =  2.3  x  10'4  cm2 
Z2  =  2.6  x  10‘4  cm2 
Z2  =  1.6  x  10"3  cm2 

~ 2  9 

r  =  2.7  x  10  J  cm'-. 
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polarized.  The  vertically-polarized,  30  cm  data  from  the  surface  of 
New  Mexico  (White  Sands  Missile  Range)  were  obtained  by  a  rocket  radar 
which  was  designed  to  transmit  vertical  polarization,  and  receive  both 
the  vertical  and  cross-polarized  components.  The  horizontally- 
polarized,  2.3  cm  results  vie  re  obtained  by  the  Ryan  13.3  GHz  Scatter- 
ometer. 

Figures  97-102  present  a  portion  of  the  radar  backscattering 
cross  section  (a  or  y)  data  for  the  lunar  surface.  The  wavelength 
dependence  of  the  circularly  polarized  backscatter  is  exhibited  in 
Fig.  97,  for  wavelengths  ranging  from  8.6  mm  to  6  m.  These  values  of 
oQ  are  obtained  by  normalizing  the  relative  power  measurements  to  a 
total  cross  section  of  0.065  times  the  geometrical  lunar  cross 
sect ion. [21] 

The  polarization  dependence  of  the  lunar  backscatter, [22]  at  23  cm 
and  68  cm,  is  shown  in  Figs.  98  and  99,  respectively.  Measurements  of 
the  circularly-depolarized  components  indicate  that  the  depolarization 
is  increasing  with  decreasing  wavelength.  It  should  be  noted  that  the 
"linear  depolarized"  measurements  at  23  cm  (see  Fig.  98),  include  the 
return  from  a  complete  range  ring  of  constant  delay,  and  hence  represent 
a  statistical  mixture  of  "vertical"  and  "horizontal"  components  relative 
to  the  local  plane  of  incidence  co-ordinates.  As  mentioned  previously, 
near  normal  incidence,  the  separation  between  direct  and  crossed  polarized 
return  tends  to  equal  the  polarization  separation  of  the  antenna  used  in 
the  experiment.  Thus  the  rise  in  the  depolarized  return  near  normal 
incidence  is  more  likely  to  be  an  artifact  of  the  experiment  rather 
than  a  property  of  the  surface. 

The  radar  cross  section  of  the  lunar  surface,  at  13.3  GHz, [23]  for 
Surveyors  V,  VI,  and  VII,  is  shown  in  Figs.  100,  101,  and  102,  respectively. 
These  data  (a  =  2.3  cm)  are  in  good  agreement  with  the  Earth-based  3.6 
cm  results[21]  (solid  curve). 
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2.  Vegetated  surfaces 


BacKscattering  cross  sections  for  a  number  of  vegetated  surfaces, 
are  presented  in  Figs.  103-104.  A  large  quantity  of  additional  data 
may  be  found  ir.  References  3  and  6.  In  Fig.  103,  airborne  X-band 
radar  measurements  of  y  for  wooded  terrain,  obtained  by  NRL[14]  and 
Goodyear, [16]  are  compared.  These  data  are  in  fairly  good  agreement. 

It  is  noted  that  the  radar  return  at  X-band  is  approximately  inde¬ 
pendent  of  tree  type  and  height. 

Figure  104  exhibits  the  backscattering  cross  sections,  at  X-band, 
for  some  typical  Arizona  farmland.  The  irrigated  farmland  ar.n 
included  the  following  crops:  onions  (less  than  one  foot  high), 
alfalfa  and  wheat  (2  feet),  and  maize  (2  to  3  feet). 

The  backscattering  cross  sections  (at  S-,  X-,  and  K  -bands)  for 
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soybeans  and  oats  are  shown  in  Figs.  105  and  106,  respectively.  These 
surfaces  are  representative  of  the  two  typical  structures  for 
agricultural  crops:  leafy  (soybeans)  and  cylindrical  (oats).  The 
curves  presented  are  the  upper  and  lower  limits  of  the  measured  data 
for  the  number  of  indicated  dates.  Only  the  vertically-polarized 
direct  component  has  been  plotted  since  the  return  is  essentially 
independent  of  polarization. 

Figure  107  presents  the  backscattering  cross  section  for  a  wheat 
surface[25]  at  297  Wz,  5.87  GF’,  and  9.375  GHz  for  horizontal  polari¬ 
zation.  A  point  of  interest  is  that  oats  and  wheat  have  a  rather 
similar  structure,  yet  the  X-band  radar  return  from  the  wheat  (Fig. 

107)  is  nearly  an  order  of  magnitude  less  than  from  oats  (Fig.  106). 

A  similar  result  was  noted  by  Peake, [12]  from  measurements  at  this 
facility. 
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3 .  Developed  land  areas  (cities) 

Figures  108  and  109  provide  airborne  measurements  of  the  back- 
scattering  cross  section,  y,  for  cities.  Measurements  by  NRL[15]  at 
428  MHz  and  Westinghouse[26]  at  220  MHz,  are  shown  in  Fig.  108,  for 
vertical  polarization.  In  Fig.  109,  X-band  values  of  y  for  cities, 
including  Chicago,  Phoenix,  and  New  Jersey  residential,  are  plotted 
from  measurements  by  NRL. [14,15] 

4.  Sea  surfaces 

The  radar  backscattering  cross  section  for  the  sea,  aQ(ei. ) 

(where  e.  is  measured  from  normal),  is  presented  in  Figs.  110-129, 
for  frequencies  ranging  from  428  MHz  to  48.7  GHz.  The  order  of  pre¬ 
sentation  is  such  that  the  :•  recent  data  are  shown  first. 

Figures  110-117  present  the  newer  NRL  data[27]  which  show  little 
dependence  of  aQ  on  wind  velocity  at  angles  between  2C°  and  70°.  In 
particular,  the  difference  in  return  for  "smooth"  and  "rough"  sea 
conditions  is  less  than  10  dB  at  all  frequencies.  X-band  backscatter 
measurements [27]  on  sea  surfaces  with  oil -slicks  are  shown  in  Figs.  118 
and  119.  together  with  data  for  four  ’'oughest  sea  conditions  in  July 
1968.  These  results  indicate  very  little  scattering  at  angles  more 
than  20°  or  30°  off  normal  incidence  for  sea  surfaces  with  oil- 
slicks,  in  agreement  with  the  concept  that  the  "diffuse"  part  of  sea 
surface  scattering  (i.e.,  the  backscattering  between  20°  and  70°) 
is  due  to  Bragg  scatter  by  the  capillary  waves.  These  are  suppressed 
by  the  oil  slick. 

The  radar  data[28,29,30]  which  show  a  more  pronounced  dependence 
of  oQ  on  wind  velocity  are  shown  in  Figs.  120-129.  At  the  present 
time,  the  discrepancy  between  these  results  remains  unresolved.  It 
seems  likely,  however,  that  a  part  of  the  difference  is  due  to  different 
experimental  geometries,  and  to  the  fact  that  some  measurements  were 
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taken  over  fulty  developed  seas  in  the  open  ocean,  whereas,  others  were 
taken  under  conditions  of  limited  fetch  (e.g..  Fig.  121).  In  any  case 
it  is  clear  that  there  is  a  significant  difference  between  sea  return 
at  very  low  wind  velocities  (say  less  than  5  knots)  and  that  at  the 
higher  wind  velocities. 

5.  Near-grazing  angle  data 

Figure  130  presents  small  grazing  angle  measurements  for  trees[31] 
(X-  and  S-band)  and  cotton  seedl ings[16]  (X-band).  The  most  significant 
aspect  of  the  data  is  the  increase  in  y  as  the  grazing  angle  decreases 
from  1  to  3  degrees.  This  tendency  of  the  backscatlering  near 
grazing  to  increase  as  grazing  incidence  is  approached  has  been 
observed  for  many  years,  but  has  received  no  general  explanantion. 

Ground  based  X-band  measurements  of  aQ(ei.),  by  Li nel  1 , [32]  for  cultivated 
and  wooded  terrain,  are  shown  in  Figs.  131  and  132,  with  the  percentile 
value  as  the  parameter  indicated  on  the  curves.  No  difference  in  the 
results  between  horizontal  and  vertical  polarizations,  was  observed 
for  either  terrain.  The  curves  for  the  cultivated  terrain  (Fig.  131) 
exhibit  a  characteristic  upturn  in  y  for  angles  of  less  than  1.5  degrees. 

The  variation  of  a  with  season  is  shown  for  cultivated  terrain 
o 

in  Figs.  133-135  and  for  a  forested  area  in  Fig,  116.  In  these 
figures,  the  grazing  angle  is  fixed  and  aQ  is  plotted  against,  the 
percentage  value.  The  point  of  interest  is  that  a  strong  variation 
of  aQ  with  season  exists  for  the  cultivated  terrain,  but  not  for  the 
forested  area.  It  was  noted  that  the  forest  consisted  primarily  of 
pines  and  firs,  which  remain  green  throughout  the  year. 

Measured  distributions  of  average  values  of  a  for  different 
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types  of  terrain,  at  grazing  angles  of  1°,  5°,  and  30°,  are  shown 
in  Figs.  137-139.  These  measurements  were  obtained  by  Ericson,[33] 
of  the  Research  Institute  of  National  Defense,  Stockholm,  Sweden, 
using  an  airborne  X-band  radar  with  a  horizontally-polarized  antenna. 
Table  I  provides  a  description  of  the  different  types  of  terrain,, 
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TABLE  I  Types  of  terrain.  (Ref.  33) 


The  Map  Shows 

Code 

Designation 

predominantly 

remainder 

0 

cultivated  land  1. 

fields, farms .gardens 
bushes, single  trees, 
etc. 

groves, woods  up  to  25% 

1 

cultivated  land  2. 

as  above 

as  above  but  up  to  50% 

2 

Swamos 

swamps, sparsely  cove¬ 
red  with  trees  and/or 
houses 

wood  up  to  25% 

3 

Open  wood 

a)  open  wood 

b)  wood 

fields, pastures, farms,  up  to 
swamps, open  wood  etc  50% 

4 

Rocks 

moss-covered  or  na¬ 
ked  rocks  or  moun¬ 
tains, possibly  co¬ 
vered  with  bushes 
and  single  trees 

wood  up  to  25% 

r~  - 

5 

Wood 

dense  wood 

as  in  3b)  above,  but  up 
to  25% 

6 

Archipelago 

terrain 

sparsely  covered 
smaller  islands 
and  holmes,  naked 
rocks 

(water) 

7 

Villa  districts 

dispersed  smaller 
houses  surrounded 
by  gardens 

8 

Small  houses 
districts 

concentrated  lower 
buildings  up  to  3 
floors  (suburbans, 
provincial  towns) 

9 

1 - 

High  houses 
districts 

Very  concentrated 
higher  buildings 
from  4  floor  and  up 
(to  about  12  floors) 
(The  center  of 
Stockholm) 

1 

Wood:  every  kind  of  tree  included  in  the  designation. 
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G.  Brightness  Temperature  Data 


Measurements  of  the  microwave  brightness  temperature,  which 
are  presented  in  this  section,  are  taken  mostly  from  investigations 
by  The  Ohio  State  Universi ty[7 ,8]  and  Space-General  Corporation. [34,35] 
The  Ohio  State  data  have  been  corrected  by  a  boot-strap  inversion  of 
the  integral  equation;  other  data  are  in  raw  form.  Data  are  given  for 
a  number  of  surfaces  in  the  following  categories:  geological,  vege¬ 
tated,  and  sea.  Finally,  a  comparison  is  made  of  the  radiometric 
temperatures  for  a  variety  of  natural  surfaces. 

1 .  Geological  surfaces 

Figures  140-142  are  measured  brightness  temperatures  for  asphalt, 
glacially  polished  limestone  (Marblehead,  Ohio),  and  a  smooth  coal  bed 
exposed  by  strip  mining  operations  (Cadiz,  Ohio).  These  surfaces 
exhibit  the  characteristic  of  a  specular  surface  (see  Section  D). 

Radiometric  deviation  envelopes  and  mean  value  curves[34]  of  barren 
mud  (Palo  Alto,  California),  without  vegetation,  are  presented  in 
Figs.  143  and  144.  (The  mean  value  curve  is  the  simple  average  of  all 
the  data  recorded  at  a  particular  angle).  These  curves  also  exhibit 
the  characteristic  shape  for  "smooth"  surfaces. 

Figures  145  and  146  are  the  deviation  envelopes  and  mean  value 
curves[34]  of  the  radiometric  temperature  for  La  Jolla  Beach  sand  at 
13.5  GHz  and  37  GHz,  respectively.  The  sand  shows  a  definite  polari¬ 
zation  dependence  at  both  13.5  GHz  and  37  GHz.  Examination  of  the 
general  shape  of  the  mean  sand  curves  shows  that  the  sand  does  not 
appear  as  a  specular  surface  or  a  totally  diffuse  surface,  but  has 
some  of  the  characteristics  of  each. 

The  brightness  temperatures  of  lapilli  (Mono  Crater,  California) 
and  bare  soil  (Ohio  State  University  Farms)  are  shown  in  Fgs.  147  and 
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148,  respectively.  As  was  the  case  for  the  La  Jolla  Beach  sand,  these 
surfaces  also  show  a  character  intermediate  between  a  smooth  surface 
and  a  diffuse  surface. 

The  radiometric  temperatures[35]  for  two  other  soil  surfaces,  silt 
loam  and  sandy  loam,  are  presented  in  Figs.  149  and  130.  These  data 
show  only  the  horizontally-polarized  component  of  the  brightness. 

Fiyures  151-153  exhibit  measurements [7]  made  at  Mono  Craters, 
California,  on  a  sequence  of  lavas  of  essentially  the  same  chemical 
composition  but  different  physical  properties,  ranging  from  a  very 
light  pumice  to  obsidian.  The  brightness  temperatures  are  almost 
independent  of  angle  of  incidence,  which  is  characteristic  of  a  Lambert 
Law  surface,  but  depend  on  the  density  (i.e.,  the  effective  dielectric 
constant)  of  the  material. 

Figures  154  and  155  are  radiometric  measurements [34]  collected 
near  Pisgah  Craters,  California,  at  13.5  GHz  and  37  GHz,  for  horizontal 
and  vertical  polarizations.  The  response  of  a  series  of  thick  playa 
sediments,  a  thick  lava  flow,  and  a  1.5  foot  layer  of  playa  sediments 
covering  the  lava  flow,  are  shown  in  each  figure.  In  all  cases,  the 
radiometric  temperature  of  the  sediments  is  seen  to  be  substantially 
different  than  that  for  the  exposed  lava. 

2.  Vegetated  surfaces 

Figures  156-164  provide  measurements,  by  the  Ohio  State 
University,  of  X-band  brightness  temperatures  for  a  number  of  native 
agricultural  crops  including  oats,  wheat,  sorghum,  and  soybeans.  The 
point  of  interest  is  that  the  oats,  wheat,  sorghum,  and  alfalfa 
surfaces  have  brightness  temperatures  ranging  from  260°K  to  290°K, 
while  the  soybean  brightness  temperatures  are  somewhat  cooler,  ranging 
from  240°K  to  270°K.  The  lower  brightness  temperatures  for  the  soybean 
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surface  may  partially  result  from  irrigation,  which  decreases  the  soil 
temperature  (for  example,  compare  soil  temperatures  for  alfalfa  in 
Fig.  160  and  soybeans  in  Fig.  161. 

The  deviation  envelopes  and  mean  values  curves[34]  for  the  bright¬ 
ness  temperature  of  marsh  foliage  (salicornia)  at  13.5  GHz  and  37  GHz 
are  presented  in  Figs.  165  and  1G6.  The  curves  at  both  frequencies 
have  the  characteristics  of  a  completely  diffuse  surface,  i.e.,  a 
brightness  temperature  almost  independent  of  angle  of  incidence  and 
polarization. 

Figures  167-170  are  measurements,  by  Space-General  Corporation, [35] 
of  the  radiometric  temperatures  for  marsh  vegetation,  pasture  grass, 
flowers,  and  celery  crops.  It  is  interesting  to  note  that  the 
brightness  temperatures  at  94  GHz  for  the  grass,  flower,  and  celery 
surfaces,  are  lower  than  at  the  longer  wavelengths. 

3.  Sea  surface 

Figure  171(a)  shows  measured  brightness  temperatures,  by  Nordberg,[36] 
of  a  "smooth"  and  "rough"  portion  of  the  Salton  Sea.  The  "smooth"  data 
agrees  quite  well  with  the  theoretical  curve  calculated  by  Stogryn[37]  for 
a  smooth  sea.  The  "rough"  data,  however,  are  significantly  higher  than 
what  would  be  predicted  for  a  rough  sea;  it  is  presumed  that  the  higher 
brightness  temperature  is  due  to  the  higher  emissivity  of  foam  patches 
over  the  rough  sea.  Figure  171(b)  shows  tower  based  measurements  by 
Hollinger[39]  at  1.4,  8.4  and  19.3  GHz,  under  conditions  where  the  wind 
velocity  and  sea  state  could  be  carefully  monitored.  These  measure¬ 
ments  exhibit  the  behavior  predicted  by  Stogryn,  and  confirm  the 
concept  that  the  brightness  temperature  is  controlled  by  the  slope 
distribution  as  long  as  foam  is  absent. 


Figures  172  and  173  are  data[34]  for  Palo  Alto  Sea  water  at  13.5 
GHz  and  37  GHz,  respectively.  Comparisons  between  13.5  GHz  and  37  GHz 
show  the  37  GHz  temperatures  are  from  30°K  to  70°K  hotter  than  those 
from  13.5  GHz. 

4.  Comparison  of  measured  data 

Figures  174-176  present  some  general  comparisons  of  measured 
radiometric  temperatures  of  a  variety  of  surfaces  at  13.5  GHz, 

37  GHz,  and  35  GHz,  for  horizontal  and  vertical  polarizations.  The 
data[34]  in  Figs.  174  and  175,  are  composites  of  the  mean  value  curves 
presented  in  the  previous  section.  These  curves  illustrate  the  effects 
of  surface  roughness  and  moisture  content  (mud  and  inundated  mud).  The 
35  GHz  data  in  Fig.  176  are  older  data[36]  which  are  seen  to  agree 
favorably  with  the  more  recent  data  of  Figs.  174-175. 
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Fig.  24.  Bistatic  scattering  behavior,  in  the  azimuth  cone, 
for  smooth  sand,  loam,  and  soybean  foliage. 
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(b) 

Fig.  25.  Experimental  and  theoretical  reflection 

coefficients;  solid  curves  are  theoretical 
values  for  smooth  surfaces  -  Reference  1. 
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Reflection  Coefficient  |R| 


Fig.  26.  Reflection  coefficients  measured  by  Sherwood  and  Gi.nzton 
for  a  tidal  flat  with  some  organic  material;  solid  curves 
are  theoretical  for  e  =  10  -  Reference  9. 


Fig.  27.  Effects  of  polarization  on  the  azimuthal  bistatic 
scattering  pattern  for  smooth  sand.  Dlstdtlc 
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Fig.  31.  Radar  backscattering  cross  sections  for  three 
smooth  surfaces  (limestone,  asphalt,  and  con¬ 
crete)  at  X-band. 


Fig.  Radar  backscattering  cross  sections  for  three 

rubble  type  surfaces  (sand,  pea  gravel,  and 
crushed  stone)  at  X-band.  C.P.  (Circular 
Polarization) 
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Fig.  33.  Radar  backscattering  cross  section  for  large  blocks 
of  pumice  at  Mono  Crater,  California,  (p  -•  density) 
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Fig.  34. 

Pumice  blocks  at  Mono  Crater.  The  photograph  shows 
the  large  angular  blocks  of  pumice  that  underlie 

*  * 

the  rugged  topography  near  the  summit  of  one  of  the 

Mono  Craters.  Note  the  large  size  of  the  vesicles. 

The  scale  is  in  inches. 
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BACKSCATTERING  CROSS  SECTION,  y  <dB) 


ANGLE  OF  INCIDENCE  (  DEGREES  ) 


Fig.  35.  Radar  backscatter  for  green  soybeans  at  X-band. 


ANGLE  OF  INCIDENCE  (DEGREES) 

Fig.  36.  Radar  backscatter  fcr  green  oats  at  X-band  - 

Reference  6.  (Stepped  data  indicates  variability 
of  measurements  made  on  same  day) 
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Fig.  37.  Measured  and  computed  difference  between  vertical  and 

horizontal  backscatter  for  •'smooth"  limestone  at  X-band. 


Fig.  33.  Measured  and  computed  difference  between  vertical  and 
horizontal  backscatter  for  "smooth"  asphalt  at  X-band. 
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HORIZONTAL  POLARIZATION 


Fig.  3S.  Angular  dependence  of  the  radar  backscatter  for  slightly  rough 
gravel;  average  diameter  of  rounded  gravel  =  1  cm.  - 
Reference  4. 


Radar  return  vs.  normalized  roughness  for  crushed  stone,,  gravel 
and  sand  surfaces  at  X-,  K  - ,  and  K  -bands  -  Reference  4. 
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Seasonal  changes  of  grass  at  X-band  -  Reference  3 
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Fig.  43.  Seasonal  changes  of  grass  at  K  -band  -  Reference  3 
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ANGLF  OF  INCIDENCE  (DEGREES) 

Fig.  44.  Effects  of  rain  on  a  smooth  asphalt  road 
at  K  -band  -  Reference  3. 
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ANGLE  OF  INCIDENCE  (DEGREES) 

Fig.  45.  Effects  of  rain  on  two-inch  grass  at  K  -band  - 
Reference  3.  a 
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BACKSCATTERING  CROSS  SECTION  , 


ANGLE  OF  INCIDENCE  (DEGREES) 

Fig.  46.  Frequency  dependence  of  radar  backscattering  for  a 
smooth  asphalt  surface.  (Roughness  -  0.3nm) 
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Fig.  47.  Frequency  dependence  of  radar  backscattering  for  a 
smooth  sand  surface  -  Reference  4. 
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BACKSCATTERING  CROSS  SECTION,  y  (dB) 
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Fig.  48.  Frequency  dependence  of  radar  backscattering  for  a 
slightly  rough  gravel  surface:  average  diameter  of 
rounded  gravel  =  1  cm  -  Reference  4. 
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BACKSCATTERING  CROSS  SECTION, 


Fig.  4g.  Frequency  dependence  of  radar  backs cattering  for  a 
rough  crushed  stone  surface;  average  diameter-  of 
scone  «  3  to  5  cm  -  Reference  4. 
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Fig.  50.  Frequency  dependence  of  radar  backscattering  for 
rough  grey  purtice  surface. 
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Fig.  51  Frequency  dependence  of  radar  backscattering 
for  a  veoe  Fated  surface  (oats). 
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olarization  dependence  of  the  radar  backscattering  for 
smooth"  limestone  surface  at  X-band  -  Reference  5. 


BACKSCATTERING  CROSS  SECTION  ,  v  (  dB )  BACKSCATTERING  CROSS  SECTION 


Fig.  53.  Polarization  dependence  of  the  radar  backscattering 
for  a  smooth  sand  surface  at  X-band  -  reference  4. 


Fig.  54.  Polarization  dependence  of  the  radar  backscattering 
for  a  slightly  rough  gravel  surface  at  X-band  - 
Reference  4. 
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58.  Brightness  temperature  of  smooth  asphalt  at  o5  GHz 


RADIOMETRIC  TEMPERATURE 


ANGLE  OF  INCIDENCE  (DEGREES) 


Fig.  59.  Brightness  temperatures  for  different  size 
fractions  of  limestone  at  13.5  GHz.  - 
Reference  11. 
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BRIGHTNESS  TEMPERATURE  (#K)  BRIGHTNESS  TEMPERATURE  C*K) 


Fig.  60.  Brightness  temperature  for  dense,  broker.,  grey 
pumice  at  10  GHz.  (p  =  1.9;  e  =  2.7) 


Fig.  61.  Brightness  temperature  for  light  weight,  dark 
pumice  at  10  GHz.  (p  =  0.45) 
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BRIGHTNESS  TEMPERATURE  (°K> 


Fig.  62.  Brightness  temperature  for  large  blocks  of  obsidian  at 
10  GHz.  (p  =  2.3;  e  «  4.7) 
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Fig.  63.  Photographs  of  volcanic  materials  at  Mono  Craters, 
California  -  Reference  7. 
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BRIGHTNESS  TEMPERATURE 
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Fig.  67.  Computed  and  measured  brightness  temperatures  for  coal 
(near  Cadiz,  Ohio),  limestone  (near  Marblehead,  Ohio), 
and  asphalt  surfaces  at  10  GHz,  vertical  polarization. 
(Data  normalized  to  3G0°K  at  the  Brewster  Angle). 
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Fig.  68.  Brightness  temperatures  of  playa  sediments, 
with  variable  moisture,  at  13.4  GHz  - 
Reference  11. 
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Fig.  69.  Brightness  temperature  of  three  pumices  (Mono  Crater, 
California)  at  X-band. 


BRIGHTNESS  TEMPERATURE 


Fig.  70.  Brightness  temperatures  for  irrigated  and  non-irrigated 
soybeans  at  X-band. 
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BRIGHTNESS  TEMPERATURE 


Fig.  72.  Brightness  temperature  of  1  api Hi  (volcanic  ash)  at 
10  GHz  and  35  GHz  (Mono  Crater,  California). 


r 


Fig.  73.  Close-up  of  lapilli  surface.  Note  the  concentration  of  coarse  pumice 
lapiili  in  botches  and  bands,  and  the  sparse  ground  vegetation  and 
scattered  pine  needles.  Scale  is  in  inches. 
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Fig.  74.  Brightness  temperature  of  "rough"  dark  pumice  at  10  GHz 
and  35  GHz  (Mono  Crater,  California). 
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Fig.  76.  Bistatic  scattering  cross  section,  in  the  azimuth 
cone,  for  smooth  sand  at  X-band. 
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Fig.  78.  (Continued) 
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Fig.  79  Bistatic  scattering  cross  section,  in  the  plane 
of  incidence,  for  rough  sand  at  X-band. 
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F'g.  79.  (Continued) 


Fig.  80 


Bistatic  scattering  cross  section,  in  the  plane  of 
incidence,  for  stubble  at  X-band.  (<p  =  0  or  180° 
except  wnere  marked) 
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Fig.  81.  Bistatic  scattering  cross  section,  in  the  plane  of 
incidence,  for  dry  grass  at  X-band. 
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Fig.  84.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  3)  at  C-band,  vertical  polarization  - 
Reference  13. 
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Fig.  85.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  1)  at  C-band,  vertical  polarization  - 
Reference  13.  (Depression  angle  =  grazing  angle  = 
90°  -  e. .) 
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Fig.  36.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  2)  at  C-band,  vertical  polarization  - 
Reference  13. 
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Fig.  87.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  3)  at  C-band,  vertically-polarized  trans¬ 
mitter  and  horizontally-polarized  receiver  - 
Reference  13. 
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Fig.  88.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  1)  at  C-band,  vertically-polarized  trans¬ 
mitter  and  horizontally-polarized  receiver  - 
Reference  13. 
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Fig.  &.  Bistatic  scattering  cross  section  of  the  sea  (sea 
state  2)  at  C-band,  vertically-polarized  trans¬ 
mitter  and  horizontally-polarized  receiver  - 
Reference  13. 
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Fig.  91.  Radar  backs cattering  cross  section  for  normal 
desert  at  X-band. 
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Fig.  92.  Radar  backscattering  cross  section  for  Pisgah 
Lava  Flow  at  X-band  and  Ku-band. 
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Fig.  93.  Radar  backscattering  cross  section  for  Bristol 
(Dry  Salt)  Lake  at  X-band  -  Reference  16. 
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Fig.  94.  Radar  backscattering  cross  section  for  Lavic 
Dry  Lake  playa  at  X-band  and  ku~bard. 
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NOT  REFRODUCIBLE 


Fig.  95.  Lavic  Lake  Playa.  The  surface  of  the  playa  is  almost  perfectly 
flat  except  for  small  scale  roughness  generated  by  contraction 
(dessication)  cracks.  These  cracks  break  the  playa  surface  into 
polygons  of  three  distinct  sizes:  A  -  measured  in  terms  of  10’s 
of  feet;  B  -  measured  in  terms  of  a  few  feet;  and  C  -  measured 
in  terms  of  a  few  inches.  In  tne  intermediate  distance  is  the 
radar  truck  being  used  to  make  calibration  measurements  and  in 
the  far  distance  to  the  left  are  the  black  basalt  lava  flows 
from  Sunshine  Crater.  The  view  is  to  the  northwest. 


191 


20 


ANGLE  OF  INCIDENCE  (DEGREES) 


Fig.  96.  Radar  backscattering  cross  sections  for  the 
lunar  and  planetary  surfaces. 
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Fig.  97.  Frequency  dependence  of  the  radar  backscattering 
cross  section  for  the  Moon  -  Reference  21. 


192 


10  LOG  er0  (dB) 


I 


IOi 


Fig.  98.  Polarization  dependence  of  the  radar  backscattering 
cross  section  for  the  Moon  at  23  cm  -  Reference  22. 


Fig.  99.  Polarization  dependence  of  the  radar  backscattering 
cross  section  for  the  Moon  at  68  cm  -  Reference  22. 
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Fig.  100.  Radar  cross  section  of  the  lunar  surface  vs. 

angle  of  incidence  for  Surveyor  V  (13.3  GHz). 


Fig.  101.  Radar  cross  section  of  the  lunar  surface  vs. 

angle  of  incidence  for  Surveyor  VI  (13.3  GHz). 


Fig.  102.  Radar  cross  section  for  the  lunar  surface  vs. 

angle  of  incidence  for  Surveyor  VII  (13.3  GHz). 
(Data  f i  om  Ref.  23) 
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Fig.  103.  Radar  backscattering  cross  section  for  wooded  terrain  at  X-band. 
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Fig.  105.  Radar  backscattering  cross  section  for  24"-36 
green  soybeans  at  S-,  X-,  and  K  -bands. 

Shaded  Area  is  envelope  of  measured  curves. 
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Fig.  107.  Radar  backscattering  cross  section  for  wheat  at  297 
MHz,  5.87  GHz,  and  9.375  GHz  -  Reference  25. 
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Fig.  110.  Radar  backscattering  cross  section  for  the  sea  at 
428  Wz,  vertical  polarization  -  Reference  27. 
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Fig.  112.  Radar  backscattering  cross  section  for  the  sea  at 
1.228  GHz,  vertical  polarization  -  Reference  27. 


Fig.  113.  Radar  backscattering  cross  section  for  the  sea  at 
1.228  GHz,  horizontal  polarization  -  Reference  27. 
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Fig.  114.  Radar  ba  ^scattering  cross  section  for  the  sea  at 
4.455  GHz,  vertical  polarization  -  Reference  27. 


Fig.  115. 


Radar  backscattering  cross  section  for  the  sea  at 
4.455  GHz,  horizontal  polarization  -  Reference  27. 
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Fiu.  116.  Radar  backscattering  cross  section  for  the  sea  at 
8.91  GHz,  vertical  polarization  -  Reference  27. 


Fig.  117.  Radar  backscattering  cross  section  for  the  sea  at 
8.91  GHz,  horizontal  polarization  -  Reference  27. 
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g.  118.  Radar  backscattering  cross  section  for  the  sea  with 

slicks  at  8.91  GHz,  vertical  polarization  -  Reference  27. 
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Fig.  119.  Radar  backscattering  cross  section  for  the  sea  with 

slicks  at  8.91  GHz,  horizontal  polarization  -  Reference  27. 
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Fig.  121 


Radar  backscattering  cross  section  for  the  sea 
at  8.8  GHz,  (average  of  vertical  and  horizontal 
polarization)  -  Reference  29. 
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Fig.  122.  Radar  backscattering  cross  section  for  the  sea  at 
9.4  GHz,  vertical  polarization  -  Reference  24. 
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Fig.  123.  Radar  backscattering  cross  section  fo.'  the  sea  at 
24  GHz,  vertical  polarization  -  Reference  24. 
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Fig.  124.  Radar  backscatteri ng  cross  section  for  the  sea  at 
35  GHz,  vertical  polarization  -  Reference  24. 
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Fig.  125.  Radar  backscattering  cross  section  for  the  sea  at 
9.6  GHz,  vertical  polarization  -  Reference  30. 
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Fig.  126.  Radar  backscattering  cross  section  for  the  sea  at 
24  GHz*  vertical  polarization  -  Reference  30. 
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Fig.  127.  Radar  backscattering  cross  section  for  the  sea  at 
35  GHz,  vertical  polarization  -  Reference  30. 
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Fig.  128.  Radar  backscattering  cross  section  for  the  sea  at 
48.7  GHz,  vertical  polarization  -  Reference  30. 
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hg.  131.  Near-grazing  radar  backscatcer  for  cultivated 
terrain  at  10  GHz  -  Reference  32.  Parameter 
is  per  cent  of  returns  exceeding  given  ordinate. 


GRAZING  ANGLE  (DEGREES' 

Fig.  132.  Near-grazing  radar  backscatter  for  wooded  ter¬ 
rain  at  10  GHz  -  Reference  32.  Parameter  is 
per  cent  of  returns  exceeding  given  ordinate. 
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Fig.  133.  Seasonal  variation  in  the  radar  backscatter 
distribution  for  cultivated  terrain,  at  10 
GHz,  for  a  1.25°  depression  angle  - 
Reference  32. 
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Fig.  134.  Seasonal  variation  in  the  radar  backscatter 
distribution  for  cultivated  terrain,  at  10 
GHz,  for  a  2.5°  depression  angle  - 
Reference  32, 
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Fig.  135.  Seasonal  variation  in  the  radar  backscatter 
distribution  for  cultivated  terrain,  at  10 
GHz,  for  a  5*  depression  angle  -  Reference  32. 
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Fig.  136.  Seasonal  variation  in  the  radar  backscatter 
distribution  for  wooded  terrain,  at  10  GHz, 
for  a  0.7°  depression  angle  -  Reference  32. 
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Fig.  137.  Distribution  of  average  values  of  o0  for  different  types 
of  terrain,  at  X-band,  horizontal  polarization,  for  a  1° 
depression  angle  -  Reference  33. 
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Fig.  I38.  Distribution  of  average  values  of  n0  for  different  types 
of  terrain,  at  X-band,  horizontal  polarization,  for  a  5° 
depression  angle  -  Reference  33. 
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Fig.  139.  Distribution  of  average  values  of  a0  for  different  types 
of  terrain,  at  X-band,  horizontal  polarization,  for  a  30° 
depression  angle  -  Reference  33. 
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(b) 

Deviation  envelope  and  mean  value  curve  for  the 
radiometric  temperature  of  Palo  Alto  Marine  mud 
at  13.5  GHz  -  Reference  34. 
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Fig.  144.  Deviation  envelope  and  mean  value  curve  for  the 
radiometer  temperature  of  Palo  Alto  Marine  mud 
at  37  GHz  -  Reference  34. 
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Deviation  envelope  and  nean  value  curve  for  the 
radiometric  temperature  of  LaJolla  Beach  sand  at 
13.5  GHz  -  Reference  34. 
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iviation  envelope  and  mean  value  curve  for  the 
idiometric  temperature  of  LaJolla  Beach  sand  at 
r  GHz  -  Reference  34. 


BRIGHTNESS  TEMPERATURE 


Fig.  H7.  Brightness  temperature  for  lapilli  (volcanic  ash)  at 
10  GHz  and  35  GHz. 
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ig.  U9.  Radiometric  temperature  for  silt  loam  soil  (Davis, 
California)  at  13.4  GHz  and  37  GHz  -  Reference  35. 
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Fig.  150.  Radiometric  temperature  for  sandy  loam  soil  with 
10%  vegetation  cover  (Rocky  Mtn.  Arsenal)  at  13.4 
GHz  and  37  GHz  -  Reference  35. 
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.  1 51 .  Brightness  temperature  for  dense,  broken  grey  pumice 
at  10  GHz  and  35  GHz. 
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Fig.  152.  Brightness  temperature  for  light  weight,  dark  pumice 
(block)  at  10  GHz  and  35  GHz. 
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Fig.  154.  Radiometric  measurements  of  playa  sediments 
lava  flow,  and  lava  flow  beneath  playa  sedi 
ments  at  13.5  GHz  -  Reference  34. 
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Fig.  i5  7.  Brightness  temperature  of  ripe  oats  at  10  GHz. 
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BRIGHTNESS  TEMPERATURE  { 


Fig.  158.  Brightness  temperature  of  green  wheat  at  1C  GHz. 
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Fig.  159.  Brightness  temperature  of  green  sorghum,  in  tassel, 
at  10  GHz. 
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rig.  162.  Brightness  temperature  of  irrigated,  green  soybeans  at 
10  GHz.  (4  September) 


ANGLE  OF  INCIDENCE  (DEGREES) 

Fig.  163.  Brightness  temperature  of  irrigated  soybeans  at  10  GHz. 
(13  September) 
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Fig.  165.  Deviation  envelope  and  mean  value  curve 
for  the  radiometric  temperature  of  Palo 
Alto  marsh  foliage  at  13.4  6Hz  - 
Reference  34. 
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Fig.  1 66  Deviation  envelope  and  mean  value  curve 
for  the  radiometric  temperature  of  Palo 
Alto  marsh  foliage  at  37  GHz  - 
Reference  34. 


Fig,  168.  Radiometric  temperature  for  pasture  grass  (Thousand  Oaks) 
at  15.4  GHz  and  94  GHz  -  Reference  35. 


Fig.  169.  Radiometric  temperatu^  for  flower  stocks  (Oxnard  at  13.4, 
37,  and  94  GHz  -  Reference  35. 
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RADIOMETRIC  TEMPERATURE 


Fig,  170.  Radiometric  temperature  for  celery  crops  (Oxnard)  at 
13.4,  37,  and  94  GHz  -  Reference  35. 
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Fig.  171(a).  Brightness  temperature  of  the  Sal  ton  Sea  at 
1.55  cm  -  Reference  36. 
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Fig,  172.  Deviation  envelope  and  mean  value  curve 
for  the  radiometric  temperature  of  sea 
water  (Palo  Alto)  at  13.5  GHz- 
Reference  34, 
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Fig.  174.  Comparison  of  mean  radiometric  temperatures  at  13.5  GHz  -  Reference  34. 
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Comparison  of  mean  radiometric  temperatures  at  37  GHz  -  Reference  34. 
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Comparison  of  measured  radiometric  temperatures 
at  35  GHz,  horizontal  polarization  -  Reference  38. 
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